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1.  INTRODUCTION 


The  main  problem  studied  in  this  thesis  is  that  of  computing  an  optimal  routing  strategy  in 
a  dynamic  network  where  the  node  demands  and  link  capacities  are  deterministic  functions  of 
time,  and  where  the  commodity  being  routed  is  represented  by  continuous  (rather  than 
discrete)  variables.  More  specifically,  we  consider  a  network  with  node  set  N  and  a  single  des¬ 
tination  node  d  €N .  where  the  link  flows  and  node  demands  are  real- valued  measurable  func¬ 
tions  on  some  time  interval  [0  J  ],  and  where  storage  is  allowed  at  the  nodes.  The  flow  utj  (t  ) 
on  link  (i  .j  )  at  time  t  represents  the  rate  at  which  the  commodity  moves  from  node  i  to  node 
j  at  time  t .  and  the  demand  rj(t)  for  node  j  at  time  t  represents  the  rate  at  which  commo¬ 
dity  is  generated  at  node  j  at  time  t .  Thus  the  storage  x}  (t  )  at  node  j  satisfies  the  differential 
equation 


Xj  (t )  =  ri  (f  )  +  £  Ujj  (r  )  -  £  Ujj  ( t  ). 

i  i 

The  link  flows  and  the  storage  at  the  nodes  are  constrained  by  (possibly  time-varying)  link 
and  storage  capacities,  and  the  main  problem  is  to  find  a  dynamic  flow  assignment  (routing 
strategy)  u  =  (,ui}  (t  ):  i  .j  €JV .  t  €[0.T ])  which  minimizes  the  cost  function 

r 

/(«)*/  £  xj(t)dt. 

0  jiN-d 

J  (u  )  is  commonly  referred  to  as  the  total  delay  due  to  u  . 

In  this  thesis,  two  approaches  are  used  to  compute  an  optimal  routing  strategy.  The  first 
approach  (valid  when  the  node  demands  and  link  capacities  are  constant  in  time)  is  to 
transform  the  problem  via  "flow  relaxation"  (Chapter  2)  into  a  finite  dimensional  convex  pro¬ 
gramming  problem  which  can  be  solved  using  gradient  descent  methods.  This  approach  leads  to 
a  decentralized  iterative  algorithm  (Section  2.5). 
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The  second  approach  is  the  one  emphasized  in  this  thesis  and  is  valid  even  when  the  link 
capacities  are  time-varying.  This  approach  is  combinatorial  in  nature,  and  makes  use  of  the 
max-flow  min-cut  theorem  for  static  networks  [Ford  and  Fulkerson.  1962].  The  key  to  this 
approach  is  explained  below. 

If  u  is  a  dynamic  flow  assignment  and  r€[0.r  ],  let  vT(u  )  denote  the  total  amount  of  com¬ 
modity  reaching  the  sink  d  up  to  time  t.  We  define  a  T-maximum  dynamic  flow  to  be  one 
which  maximizes  v  r(u  )  over  all  dynamic  flows.  When  the  link  capacities  and  node  demands 
are  constant  or  piecewise-constant,  the  problem  of  finding  a  r-maximum  flow  is  equivalent  to  a 
static  maximum  flow  problem.  The  key  to  the  second  approach  is  the  discovery  that  a  dynamic 
flow  has  minimum  total  delay  if  and  only  if  it  is  r-maximum  for  all  r€[0.r  ].  In  view  of  this 
fact,  the  max-flow  min-cut  theorem  can  be  used  to  divide  the  problem  of  finding  a  minimum- 
delay  dynamic  flow  into  two  subproblems.  This  is  the  main  idea  behind  the  algorithms  of  Sec¬ 
tions  2.6  and  3.5.  which  compute  a  minimum  delay  dynamic  flow  by  solving  a  series  of  static 
maximum  flow  problems. 

In  the  case  of  constant  capacities,  these  algorithms  require  only  O  (  I  N  1 4)  arithmetic 
operations,  which  is  very  efficient  considering  that  the  fastest  algorithm  for  computing  a  static 
maximum  flow  in  dense  networks  requires  O  ( I N  1 3)  arithmetic  operations. 

In  order  to  apply  the  above  idea  to  dynamic  networks  with  continuously-varying  capaci¬ 
ties.  we  must  generalize  the  max-flow  min-cut  theorem  so  that  it  can  be  applied  to  such  net¬ 
works.  This  is  done  in  Chapter  4,  where  we  define  a  continuous  network  (called  a  product- 
algebra  network)  whose  flows  and  capacities  are  additive  set  functions.  We  prove  a  max-flow 
min-cut  theorem  for  product-algebra  networks  (which  generalizes  the  classical  max-flow  min- 
cut  theorem)  and  establish  the  existence  of  a  flow  in  such  networks  which  simultaneously  max¬ 
imizes  the  flow  into  an  increasing  family  of  node  sets. 

The  results  of  Chapter  4  are  applied  in  Chapter  5  to  establish  a  version  of  the  max-flow 


m in-cut  theorem  for  dynamic  networks  with  continuously-varying  capacities,  and  to  establish 
the  existence  of  a  dynamic  flow  in  such  networks  which  is  r-maximum  for  all  r€[OJ’].  It  is 
hoped  that  these  results  will  contribute  to  the  understanding  of  dynamic  networks  and  may 
lead  to  new  methods  for  computing  optimal  dynamic  flows. 

In  Chapter  6,  we  depart  from  the  study  of  dynamic  networks,  and  we  generalize  the 
theory  of  continuous  networks  presented  in  Chapter  4.  We  no  longer  require  that  the  capacities 
be  additive  set  functions.  Instead,  the  flows  are  identified  with  linear  functionals  and  are  con¬ 
strained  by  a  sublinear  functional.  Using  the  Hahn-Banach  Theorem,  a  generalization  of  the 
max-flow  min-cut  theorem  is  proved  under  a  certain  condition  on  the  constraining  sublinear 
functional.  As  an  application,  we  consider  a  class  of  finite  networks  which  is  related  to  the 
polymatroid  networks  of  [Lawler  and  Martel.  1982], 


2.  MINIMUM  DELAY  DYNAMIC  ROUTING:  CONSTANT  CAPACITIES 


2.1.  Introduction 

We  consider  a  minimum  delay  dynamic  routing  problem  using  a  network  model  intro¬ 
duced  by  Segall  (1977).  where  the  flow  on  each  link  is  a  deterministic  real-valued  function  of 
time  and  where  storage  is  permitted  at  the  nodes.  This  chapter  is  organized  as  follows.  The 
dynamic  routing  problem  is  formulated  in  Section  2.2  following  the  Segall  model.  In  Section 
2.3  the  concept  of  flow  relaxation  is  introduced  and  an  algorithm  for  its  efficient  computation  is 
described.  In  Section  2.4  the  problem  of  flnding  an  optimal  time-dependent  control  is 
transformed  via  flow  relaxation  into  a  vector  optimization  problem  with  convex  cost  and  sim¬ 
ple  linear  equality  and  inequality  constraints.  As  a  corollary,  new  characterizations  are 
obtained  for  arbitrary  optimal  controls. 

Three  algorithms  are  provided  in  Sections  2.5  -  2.7  for  solving  the  vector  optimization 
problem.  The  algorithm  in  Section  2.5  is  a  decentralized  iterative  algorithm,  while  the  algo¬ 
rithm  in  Section  2.6.  called  OPTFLO.  is  a  combinatorial  algorithm.  OPTFO  reduces  the  problem 
to  a  series  of  max-flow  problems  and  computes  the  exact  optimal  vector  in  at  most  O  ( \N  1 4) 
computations.  The  third  algorithm,  given  in  Section  2.7,  reduces  the  vector  optimization  prob¬ 
lem  to  a  search  for  "bottlenecks"  in  the  network.  This  third  algorithm  is  not  yet  computation¬ 
ally  competitive  with  the  other  two.  but  is  included  for  conceptual  purposes. 

Since  the  state  evolution  in  our  model  is  deterministic  (given  the  control  and  the  initial 
state),  optimal  controls  can  be  implemented  in  closed  —loop  or  open  —loop  form.  Closed-loop 
implementation  requires  that  an  optimal  control  value  be  found  for  each  possible  state  (i.e.,  a 
feedback  policy  is  needed).  This  problem  is  studied  in  the  series  of  papers  [Segall.  1977.  1979], 
[Moss  and  Segall.  1978,  1982],  [Moss.  1977],  [Jodorkovsky  and  Segall.  1979.  1981].  In  these 
papers,  max-flow  algorithms  were  used  to  compute  closed-loop  control  laws  for  networks  with 
less  than  10  nodes.  However,  it  appears  that  both  the  computational  requirement 
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and  the  storage  requirement  for  the  feedback  policy  itself  grow  exponentially  with  the  size  of 
the  network,  and  therefore  the  method  is  impractical  for  large  networks.  Since  computation  of 
closed-loop  policies  can  be  made  off-line,  the  storage  requirement  is  the  most  critical  limitation. 

Thus  we  are  led  in  this  chapter  to  seek  open-loop  solutions  whereby  the  optimal  control  is 
computed  (usually  in  real  time)  as  a  function  of  the  given  initial  state.  Since  the  optimal  con¬ 
trol  is  sought  for  only  a  single  initial  state  it  is  hoped  that  the  computation  time  (which  is  a 
critical  parameter  for  open-loop  implementation)  will  be  much  smaller  than  that  required  for 
closed-loop  solutions.  We  believe  that  the  algorithms  presented  in  this  chapter  bear  this  out. 

An  interesting  algorithm  for  solving  the  open-loop  problem  is  given  by  Shats  and  Segall 
(1980).  Our  algorithm  OPTFLO  is  more  efficient,  primarily  for  two  reasons.  First,  although 
both  the  algorithm  of  Shats  and  Segall  and  OPTFLO  are  recursive,  only  OPTFLO  divides  the 
original  problem  into  independent  subproblems  during  each  iteration.  Second.  OPTFLO  takes 
advantage  of  existing  special-purpose  algorithms  for  solving  max-flow  problems;  whereas  the 
algorithm  of  Shats  and  Segall  requires  the  solutions  of  linear  programming  problems  which  are 
solved  by  less  efficient  general  linear  programming  methods. 

Yet  another  approach  to  the  open-loop  dynamic  routing  problem  is  to  begin  by  considering 
a  discrete-time  formulation.  Then  if  the  network  is  time  expanded  [Ford  and  Fulkerson.  1962] 
into  a  new  network  which  contains  a  duplicate  of  the  original  network  for  each  unit  of  time, 
then  the  discrete-time  version  of  the  problem  we  consider  in  this  chapter  is  equivalent  to  a  cer¬ 
tain  static  minimum-cost  flow  problem.  However,  the  well-known  technique  of  "building  up" 
optimum  solutions  for  such  static  problems  [Ford  and  Fulkerson]  is  much  less  efficient  for  this 
particular  minimum-cost  problem  than  is  the  discrete-time  version  of  our  algorithm  OPTFLO. 
Indeed,  the  technique  of  Ford  and  Fulkerson  requires  the  solution  of  at  least  one  max-flow 
problem  per  time  step;  whereas  OPTFLO  often  requires  many  fewer  since  it  exploits  the  fact 
that  an  optimal  control  will  often  be  constant  over  many  consecutive  time  intervals. 
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2.2  Problem  Formulation 

A  single  destination  network  N  is  a  quadruple  (N ,  L  ,C  ,d)  where  (N .  L)  is  a  finite 
directed  graph  with  a  set  N  of  nodes  and  a  set  L  of  links  with  L  C.N  xN  ,d  is  a  distinguished 
node  in  N  called  the  destination,  and  C  =  (Ct :  l  €L  )  is  a  capacity  assignment  vector  such  that 
C,  ^ 0  for  each  link  Z .  The  notation  I A  I  will  be  used  to  denote  the  number  of  elements  in  a 
set  A  .  A  link  leading  from  node  i  to  node  j  will  often  be  denoted  (i  ,j ).  Define  E  (i )  to  be  the 
collection  of  nodes  k  such  that  (i  Jc  )  is  a  link,  and  7  (i  )  to  be  collection  of  nodes  k  such  that 
( k  ,i  )  is  a  link.  For  each  node  i  in  N  .  xt  (t  )  is  a  real  value  which  denotes  the  amount  of  traffic 
at  node  i  at  time  t  (measured  in  bits,  packets,  vehicles,  or  messages,  for  example). 

A  demand  for  the  network  is  a  pair  ( x  (0).  r  )  where  x,  (0)  for  i  €N  denotes  the  (nonnega¬ 
tive)  initial  amount  of  traffic  at  node  i ,  and  r,  for  i  denotes  the  (nonnegative)  rate  at 
which  traffic  enters  node  i  from  outside  the  network.  Given  a  control  u  =  (u,  (r  ):  l  67. .  t  ^0). 
where  u,  (t  )  denotes  the  instantaneous  flow  on  link  Z  at  time  t .  the  state  equation  of  the  net¬ 
work  is 


xi(t)  =  ri  +  £  £  Ujy(r).  i^d. 

We  emphasize  that  r,  is  not  time-dependent.  By  convention,  xd  (t  )=rd  =0  for  all  t .  The  state 
equation  can  be  written  in  vectc  /  notation  as 

x  (*  )  =  r  +  Bu  U  )  (2.1) 

where  B  is  the  node-arc  incidence  matrix  of  the  graph  (N ,  L  )  with  the  row  corresponding  to  d 
deleted.  It  is  convenient  to  define  Ci}  =0  for  pairs  of  nodes  i  .j  such  that  a  link  (i  ,j  )  does 
not  exist. 


Throughout  this  chapter,  when  sets  of  nodes  are  used  as  subscripts,  the  summation 


convention  applies.  Thus. 


rA~'Lr,. 

iiA  itA 

and 

Cab  =  £  £  Cij . 
i  €  .4  ;  €  B 

and  the  state  equations  can  be  written  as 

xA  (t)  =  rA  +  uN-A  A  (f  )  —  uA  j,_A  ( t  ),  A  CN  —d. 

Note  that  in  this  last  equation  we  have  abused  notation  by  writing  d  when  we  really  mean  the 
set  \d  }. 

A  state  trajectory  ( x  (r  ):  t  ^0)  is  well-defined  by  the  state  equation  (2.1)  for  any  demand 
(jc  (0),  r  )  and  any  control  u  in  the  set 

U  >  {u :  u  is  a  measurable  function  on  R  +  and  O^u  } 

where  vector  inequalities  such  as  are  to  be  interpreted  coordinate-wise.  A  control  u 

in  U  is  termed  admissible  for  the  demand  (x  (0),  r )  if  the  corresponding  state  trajectory 
satisfies  the  constraint  x  (f  )^0  for  all  t .  For  such  a  control  the  number 

oo 

/  (u  )  =  f  x.v  ( t  )dt 
o 

is  the  total  delay  in  the  network  incurred  by  all  traffic.  Given  a  network  N  and  demand 
(x  (0).  r  ).  we  consider  the  problem: 


minimize  {/(u  ):  u  €U,  u  is  admissible} 


(P) 


and  we  let  J*  denote  the  minimum  delay.  Since  the  input  rates  ri  are  not  time  varying,  if  the 
total  delay  is  finite  for  some  control,  then  it  is  possible  to  empty  each  of  the  nodes  in  finite 
time.  Thus,  the  problem  we  consider  might  be  termed  an  optimal  evacuation  problem. 


23.  Flow  Relaxation 

Suppose  u  is  a  control  in  U  which  is  not  necessarily  admissible  for  a  given  demand 
(x  (0).  r  ).  Then  a  control  u  is  defined  to  be  a  relaxation  of  u  for  (x  (0).  r  )  if  u  is  admissible 
and  if.  letting  x  denote  the  corresponding  (nonnegative)  state  trajectory,  the  following  condi¬ 
tions  hold  for  each  t  ^0: 

0  ^  u  (t )  ^  u  (r ).  (2.2) 

Ujj  (r  )  =  Uij  (t  ).  if  x,  (t  )  >  0.  (2.3) 


Thus  a  relaxed  control  differs  from  the  original  control  only  in  that  nodes  with  zero  traffic  may 
have  reduced  flow  on  outgoing  links,  which  may  be  necessary  to  keep  x  nonnegative.  It  will  be 
shown  that  any  control  in  U  has  at  least  one  relaxation  for  any  given  demand,  but  first  we 
show  how  a  relaxation  of  a  piecewise-constant  control  can  be  computed  using  a  series  of  "flow 
relaxations."  On  a  first  reading,  the  reader  may  now  wish  to  skip  to  the  next  section. 

Given  a  nonnegative  flow  vector  /  for  a  directed  graph  (N  ,L  ).  a  node  j  is  said  to  be  p 
links  upstream  from  node  i  (for  /  )  if  for  some  path  (i.e.,  sequence  of  nodes) 

/  =*0.*! . np  ~  i  •  f  nknk+l  >  OforO <p.  If  i  -  j  and  p  is  at  least  one  (equivalently 

two),  then  the  path  is  called  a  loop  for  /  .  The  flow  /  is  called  loop-free  if  it  has  no  loops. 
The  algorithm  implied  by  the  next  proposition  will  require  that  the  input  flow  be  loop-free.  An 
O  ( I N  I  I L  I )  algorithm  for  removing  loops  can  be  found  in  the  appendix  of  the  paper  by  Hajek 


*  ! 
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and  Ogier  ( 1984). 

Given  a  nonnegative  flow  vector  /  =  (/; :  Z  €Z, )  on  the  network  (N  .  L  ,C  .  d  )  with  input 
rate  vector  r  .  and  given  a  subset  O  CN  .  a  flow  /  is  called  a  flow  relaxation  of  /  over  O  (for 
r  )  if 

0  <  /  <  /  .  (2.4) 

f  ij  ~  f  >j .  >f  *  €ZV  — O  or  Bf ,  +  r,  >0.  (2.5) 

Bf,+ri>  0.  if  i  €0.  (2.6) 

The  next  proposition  gives  an  iterative  procedure  in  which  a  flow  relaxation  of  /  over  0  is 
computed  essentially  by  the  network  itself.  In  each  step  of  the  procedure,  every  node  i  in  O 
such  that  the  net  flow  into  i  is  negative  reduces  its  flow  on  outgoing  links  so  that  the  net  flow 
into  i  becomes  zero.  This  step  must  be  repeated  since  the  flow  on  links  going  into  nodes  in  O 
may  have  been  reduced  by  other  nodes  in  O .  If  the  flow  /  is  loop-free,  then  a  relaxation  is 
obtained  after  a  finite  number  of  steps.  The  relaxed  flow  may  not  be  unique  since  any  priority 
among  outgoing  links  may  be  used  in  choosing  how  the  outgoing  flow  is  reduced. 

Proposition  2.1 

Given  a  nonnegative  flow  vector  f  and  a  subset  O  €  N ,  let  Tf  denote  the  set  of  flow  vec¬ 
tors  /'  such  that  if 

i  €0  and  £/*i  -2 'J ij  +  r>  <  °> 

*  ) 

then 

0  *  f'i)  <  f  ij  and  £/* *  =  +  rt  *0, 

*  ) 

and  f\j  ~ftj  for  all  other  i .  Then  Tf  is  nonempty  for  every  /  ,  and  if  (/  "  :  n  ^0)  is  a 


sequence  of  flows  such  that  /  °—f  and  f  n+l€Tf  n  .  then  /  n  converges  to  some  flow  relaxation 
/  of  /  over  O .  If  /  is  loop-free  then  /  n  —f  for  all  n  ^  I N  I  and  so  /  can  be  computed  in 
O  ( I N  I  I L  I  )  computations. 

Proof 

Tf  is  clearly  nonempty  for  every  /  .  Since  0^fn+1^fn  ./  "  converges  to  a  limit  / 
which  satisfies  (2.4).  If  i  €N— O  .  then  f  ”j  +1= / .  and  so  f ,,  =f  t] ,  establishing  part  of  (2.5). 
If  Bf  f+r,  ^0  for  some  n  and  some  node  i .  then 

Bf?+l  +  n  =  Zf &+1  - Zf u+1  +  n  <  LfZ-LfZ+'  +  n  <  o. 

k  j  k  j 

and  therefore  Bf,  +r,  ^0.  Thus  if  Bf,+r,> 0,  then  Bff+r,  >0  for  each  n:  therefore. 
/  /) +1=/  [)  and  thus  fij  ~fij  for  all  i  ■  This  establishes  the  rest  of  (2.5).  If  i  €0  .  then 

Bfr  +  r,  >  L/*+1  -  Zfij  -0. 

j  J 

which  implies  condition  (2.6).  Therefore  /  is  a  relaxation  of  f  over  O . 

The  last  statement  of  the  proposition  will  be  proved  by  contraposition,  and  so  we  suppose 
that  f  I) +l  <  f  [)  for  some  link  ( i  .j  )  and  some  rt  ^  i  N  I .  Then  i  must  be  in  O  and 

Zfk.  -ZfZ  <  o  <  Z  fir1  -  Zf!)  +  n 

‘  j  <  j 

and  so  f *,  <  /Z1,-1  for  some  node  k  such  that  fa  >  0.  Continuing  by  induction,  we  deduce 
that  for  any  m  with  1  < m  (in  particular,  m  =  \N  I  )  there  exists  a  node  which  is  m  links 
upstream  from  node  i  for  flow  /  .  This  implies  that  /  is  not  loop-free.  ■ 

The  next  proposition  shows  how  a  relaxation  of  any  constant  control  (and  therefore  of  any 
piecewise-constant  control  as  stated  in  the  corollary)  can  be  constructed  from  a  series  of  flow 
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relaxations.  Together  with  Proposition  2.1.  this  gives  a  procedure  for  computing  a  relaxation  of 
any  piecew ise-constant  control.  The  procedure  is  particularly  efficient  for  constant  controls 
which  are  loop-free,  which  is  a  case  of  special  interest  in  later  sections. 

Proposition  2.2 

Let  u  be  a  constant  control  in  U  such  that  u  (r  )=/  for  all  t  ^0.  Then  a  relaxation  u  of  u 
for  demand  (x  (0).  r  )  can  be  constructed  as  follows: 

u  (t )  =  /  (k  )  t  €[f* .  tk  +1) 

where  1 0=0.  and  the  constants  /  (0),  t  ,.  /  (l).  1 2 . are  defined  recursively  by  (letting  x 

denote  the  state  trajectory  for  u  and  letting  /  (—1)=/  ) 

/  (k  )  is  a  flow  relaxation  of  /  (i  — 1)  over  Ok  .  where  Ok  =  {i  €W :  x*  ( tk )  =  0}. 

tk  +1  =  min  It  >tk :  x*  ( t  )=0  for  some  i  €A/  — Ok  }. 

Let  K  ~  min  {k :  tk  +1=+oo(.  Then  K  ^\N—d  I  and  if  /  is  loop-free,  then  u  can  be  computed 
in  at  most  O  ((A-  +1)  I N  I  I  L  I  )  computations. 

Proof 

Using  induction  on  k  .  the  definition  of  Ok ,  and  property  (2.6)  of  flow  relaxations,  we  see 
that  (u(t ),  x  (t ))  is  well  defined  and  x(f)>0  for  t  in  [r*  ,  tk  +1)  for  each  k  ^0.  Now  if 
x,  (tt  )=0  for  some  node  i  and  some  k^  1  then  ( Bf  (k  — l)+r  ),  ^0  and  therefore,  since  f  (k  )  is 
a  flow  relaxation  of  /  (£—1),  (Bf  (k  )+r  =0.  Thus  the  sequence  of  sets  Ok  is  strictly 

increasing  for  l^k  .  and  so  K  ^  \  N—  d  I.  Therefore  (u  (t ),  x  (t ))  is  well  defined  for  all 
t  >0.  Since  /  (k)^f  for  all  k  .  if  /  is  loop-free  then  so  is  /  ( k  ).  Therefore  the  last  state¬ 
ment  is  a  consequence  of  the  last  assertion  of  Proposition  2.1.  ■ 
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Corollary  2.1  (Relaxation  of  piecewise-constant  controls.) 

Let  u  be  a  piecewise-constant  control  in  U  with  the  structure  v, 

u(t)  =  f  (k).  t  €[j*  ,  ri+1) 

where  s  0=0  and  st  increases  to  infinity.  Then  a  relaxation  u  of  u  for  demand  (x  (0).  r  )  can  be 
constructed  as  follows: 

u  (f )  =  uk  (r  —  sk  ),  t  €[r* .  sk  +1) 

where  the  segments  ( x  (r  ):  t  €(sk  ,  sk  +1])  and  controls  uk  are  defined  recursively  by 

uk  is  a  relaxation  of  the  constant  control  which  is  identically  equal  to  /  (k  ),  for  demand 
( x  ).  r  ). 

Proposition  13  (Relaxation  of  general  controls.) 

(a)  A  relaxation  u  exists  for  any  control  u  in  U  and  any  pair  (x  (0),  r  ). 

(b)  Suppose  that  u  is  admissible  for  (x  (0),  r  )  and  that  u  is  a  relaxation  of  u  for  some  other 
demand  (x  (0),  r  )  such  that  x  (0)^x  (0)  and  r  .  Then  x  (t  )^x  (r  )  for  all  t  ^0. 

Proof 

Given  a  control  u  in  U  there  exists  a  sequence  of  piecewise-constant  controls  un  such  that 
u—un  converges  to  zero  in  L  +.  RL  ).  By  Corollary  2.1  there  exists  a  relaxation  of  each 
control  in  this  sequence,  thus  producing  a  sequence  of  admissible  controls.  Now  the  collection 
of  all  admisible  controls  is  a  closed,  bounded  subset  of  the  space  Z,°°(l?+.  FL)  of  bounded 
measurable  functions  from  J?+  to  JRL .  Hence,  by  the  Banach- A laoglu  theorem  [Rudin,  1973] 
there  exists  a  subsequence  of  this  sequence  of  admissible  controls  which  converges  in  the  weak* 
topology  of  Lco(It+,  EL)  to  a  control  u.  Since  L°°  is  the  dual  of  Ll  [Rudin,  1974]  the  weak* 
convergence  implies  that  the  corresponding  state  trajectories  converge  uniformly  on  bounded  t 


sets.  It  is  then  easily  shown  that  u  is  a  relaxation  of  u  if  u  is  redefined  on  a  subset  of  it  +  with 
zero  measure. 

We  now  turn  to  the  proof  of  (b).  By  definition,  the  state  trajectory  xt  for  any  node  i 
is  absolutely  continuous  with 


x,  (t  )  =  r,  +  Ink,  )• 

t  j 

Therefore,  by  conditions  (2.2)  and  (2.3).  if  xt  (r  )>0  then 

x,  (c )  ^  r,  +  (r )  —  (r )  =  x,  (r ). 

*  j 

This  and  the  facts  x  (t  )^0  and  x  (r  )^0  for  all  t  imply  (b).  ■ 

2.4.  Reducing  the  Problem  to  Finding  an  Optimal  Initial  Flow 

Fix  an  initial  state  x  (0)  and  input  rate  vector  r  throughout  this  section.  For  any  set  of 
nodes  A  and  t  ^0  define 

z  (A  X  )  =  (0)  -  t  (CA  jt—A  ~  rA  ). 

Note  that  z(A  X  )=0  if  A  is  the  empty  set. 

Lemma  2.1 

For  any  admissible  control  u  with  corresponding  state  trajectory  x  . 
xA(t)Zz(Ax) 


for  all  t  ^0  and  all  A  CN  —d  .  In  particular. 
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xN  (t )  ^  max  { z  (A  X):  A  C  N  —d  }  (2.7) 

for  all  t. 

Proof 

By  the  state  equations  and  control  constraints.  xA  (t  )^rA  — CA  _y-A  ,  and  the  lemma  is 
immediate.  ■ 

An  immediate  consequence  of  the  lemma  is  that  the  minimum  cost  J*  for  problem  P  is 
bounded  below  by 

CO 

J  max{z  (A  x):  A  CN  —d  \dt. 

o 

Next,  an  upper  bound  to  J*  will  be  identified. 

Given  a  flow  vector  /  =(/z :  l  €L  )  such  that  0^ /  .  let  uf  be  the  control  in  U  such 

that  uf  {t)  is  equal  to  /  for  all  t.  Suppose  that  uf  is  a  relaxation  of  for  the  pair 
(x  (0).  r  )  and  that  x  is  the  corresponding  state  trajectory.  Then  x,  (f  )^Bf ,  +r,  for  a.e.  t  such 
that  XjCt  )>  0.  and  so 

Xi  (t  )  <  [xz  (0)  +  t  (Bf  i  +  rt  )]+. 

The  vector  /  is  said  to  be  feasible  if  0^  /  ,  Bf  +r  ^0,  and  Bf ,  +r,  <0  for  nodes  i  such 

that  x,  (0)>0.  If  /  is  feasible  then  by  the  above  inequality  the  cost  J  (uf  )  is  no  greater  than 
Jo(f  ).  where 

oo 

Uf  )  =  /  I  [Xi  (o)  +  t  {Bf,  +  rt  )\+dt 

0  i$N 

_  V2X,  (0)2 
I  :  (0)»*0  Bfi+ri 


Therefore,  an  upper  bound  for  the  minimum  cost  /"for  problem  P  is  the  minimum  cost  J’0  for 
the  problem  P0  defined  by 

minimize  {y 0C /  ):  /  is  feasible  ).  (P0) 

If  no  feasible  vector  /  exists  we  set  J‘0  =+oo. 

Throughout  this  section  we  assume  that  there  exists  a  feasible  vector  /  .  Since  /0  is  con¬ 
tinuous  and  the  set  of  feasible  /  such  that  /0(/  )^7  is  nonempty  and  compact  for  sufficiently 
large  J  .  there  exists  a  solution  /  to  problem  P0. 

Theorem  2.1 

Suppose  that  /  is  a  solution  to  problem  P0.  Then  any  relaxation  uf  of  uf  (for  the 
demand  (x(0).  r ))  is  an  optimal  control  for  the  original  problem  P.  The  minimum  cost 
J*  =J  ( u  *  )  satisfies  J*  =J'0  ■ 

The  proof  of  Theorem  2.1  will  follow  two  lemmas. 

Remark 

Theorem  2.1  shows  that  problem  P  can  be  solved  in  two  steps.  The  first  step  is  to  solve  the 
vector  optimization  problem  P0  to  obtain  /  and  the  second  step  is  to  compute  a  relaxation  of 
uf  .  However,  since  a  relaxation  of  uf  can  be  computed  using  a  simple  decentralized  algo¬ 
rithm.  for  most  practical  purposes,  knowing  /  is  just  as  good  as  knowing  an  optimal  control  u 
in  advance.  In  fact,  because  of  the  many  degrees  of  freedom  available  in  determining  a  flow 
relaxation,  knowledge  of  /  may  even  be  preferable. 

Lemma  1.2  (Optimality  conditions  for  problem  P0.) 

A  feasible  vector  /  is  a  solution  to  problem  P0  if  and  only  if  there  exists  a  vector 
t=(t,  :  t  )  with  rt  ^0  for  i  such  that 
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rd  =0. 

T,  =  — Xj  (0)/(B/  j  +r;  ).  Bfi+ri< 0. 


and  for  every  link  (i  .j  ). 


f,j=c'j ■  Ti>'rj- 
fij-0.  t,<t  j. 


Proof 


Suppose  /  is  feasible.  Since  /y  appears  in  at  most  two  terms  in  the  sum  defining  J0(f  )• 
we  find  that 


377  2  7  ‘ 


where 


—x,  (0)/(Bf  j  +rj  ).  Bf ,  +r,  <0, 

0.  (0)=0. 


Then  the  Kuhn-Tucker  stationarity  condition  for  problem  P0  becomes 


Mir  2  —  Vir,2  +  otj  —  a,  +  Xy  —  Mij  =  0  for  all  links  (i  .j  ) 


where  a,  X.  and  /x  are  Lagrange  multipliers  corresponding  to  the  respective  constraints 
Bf  +r  ^0.  /  — C  ^0.  and  —  /  ^0  which  must  satisfy  the  complementary  slackness  conditions 


«i  ^0.  a;  =0,  unless  Bf  i  +rt  =0, 

Xy  ^0.  Xy  =0,  unless  /  y  =0. 

Mu  >°.  Mi;*0,  unless /y=Cy. 


Define  the  new  variables  r,  by 


O  I 
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(2a,  )*  x,  (0)=0. 
Tj.  X,  (0)5*0. 


Then  the  conditions  (2.8)  and  (2.9)  are  equivalent  to  those  of  the  lemma.  The  Kuhn-Tucker 
conditions  (2.8)  and  (2.9)  are  sufficient  as  well  as  necessary  for  optimality  since  the  cost  J0(f  ) 
is  convex  and  /  ranges  over  a  convex  set.  ■ 

Remark 

Note  that  if  Bf  >  +r,  <0.  then  r,  is  the  time  node  i  would  first  empty  if  the  flows  were 
fixed  to  be  /  for  all  time  (ignoring  state  constraints). 

Lemma  2.3 

Let  /  be  an  optimal  solution  to  problem  P0  (with  finite  cost)  and  let  T  and  /  satisfy  the 
optimality  conditions  of  Lemma  2.2.  Fix  <x>0.  Then  z  (A  .  cr)  is  maximized  over  A  C.N  —d  by 
a  set  R  C.N  —d  if  and  only  if  the  three  following  conditions  hold: 

/  r  „v  -r  -Cr  jv  -s  and  f  n-r  a  =0< 

(b)  r,  ^cr  for  i€N—R  withxi(0)>0, 

(c)  for  i  €R  with  Xj  (0) > 0. 


Proof 

For  A  C.N  —d  . 

*a  =  Z  x,  =  £  t,  (-5/ ,  ~r,  ) 

i  €  A  i  €  A 

and 


°^rA  ~Ca  A -A  )  =  <r(rA  /a  ft  —A  +fff-A  A  )  +  <K/  A  Jt-A  ~Ca  Jf-A  ~f  N -A  A  ) 

~  ~  T,  —  r, )  +  <r(f  a  a -a  ~~Ca  a -a  ~f \-a  a  )• 

<  €.4 
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Therefore. 


z{A  ,cr)  =  £  (r,  -cr)(-B/,-r,)  +  cr{f  A^_A-CA^.A-f  y_A  A  ).  (2.10) 

Now  define  S  =  {i  €  N :  f  j  ^o-}.  Then  the  optimality  conditions  for  problem  P0  imply  that 

fsy-s  =Csy-s  indf.v-ss  =  0-  (2.11) 

By  equations  (2.10)  and  (2.11). 

z  ( A  .cr)  -  z  (S  .cr)  =  £  (r,  -c r)(-5/ ,  -r,  )  ..v 

i€A:  T,  «r  K  '  ' 

"*■  X  (cr— r,  )(—Bf  i  —r, )  +  c r(  /  4  v  -.4  “ C.a  _y~A~f  x  -.4 .4  )• 

i€  V  —A  — d  :  r,  ><r 

Now  for  all  i  £N—d  . 


>0. 

^0. 


x*  >0. 

Xi  =0. 


Thus  all  terms  on  the  right-hand  side  of  (2.12)  are  nonpositive  and  are  equal  to  zero  if  and 
only  if  A  is  equal  to  a  set  R  which  satisfies  conditions  (a)  -  (c).  ■ 

Proof  of  Theorem  2.1 

Suppose  that  /  is  a  solution  to  problem  P0.  We  know  that 


OO 

f  max(z  (A  X  ):  A  C.N  —d  \dt  <7*  ^J(uf) 

0 

oo 

*  /  Z  [*,  (0)  +  t  (Bf ,  +r,  ))+dt  =/„(/)  =  / o- 

0  iCN 


If  the  two  integrals  in  this  string  of  inequalities  are  shown  to  be  equal,  then  it  will  follow  that 


J*  =  J  {uf  )  =  J  o  and  the  theorem  will  be  proved.  Hence,  if  for  fixed  t  >0  we  let  R  denote  the 
set  which  maximizes  z  (A  .t  )  over  all  subsets  A  of  N  —  d  ,  it  suffices  to  prove  that 


z(R  . r )  =  X  (0)  +  t  (Bf ,  +r,  )]+• 

i€,V 

For  any  node  i  with  x,  (0)>0,  the  optimality  conditions  imply  that  r,  =  —  x,  (Q)/(Bf ,  +r,  ).  and 
so 


x,  (0)  +  t  (.Bf  ,+r,  )  =  Xi  (OXl—  t  /t,  ). 

It  follows  from  this  and  conditions  (b)  and  (c)  of  Lemma  2.3  that  for  nodes  i  with  x,  (0)>0. 

x,  (0)  +  t  (Bf ,  +r, )  >  0,  i€R. 

x,(0)  +  t(Bf,+r,  )  <  0.  i€N-R. 

These  relations,  along  with  condition  (a)  of  Lemma  2.3,  imply  the  three  equalities 

z  (R  X  )  =  xR  (0)  —  t  (/ g  ji-g  —fx-g  jg  —rR  ) 

=  £  *,  (0)  +  t  ( Bf ,  +r,  ) 

i€  J? 

=  X  (0)  +  t  {Bf  i  +r-j  )]+. 

i  €  y 

and  the  proof  is  complete.  ■ 

Corollary  2.2 

For  any  admissible  control  u  with  corresponding  trajectory  x  ,  u  is  optimal  if  and  only  if 
equality  holds  in  (2.7)  for  the  particular  optimal  control  uf  found  in  Theorem  2.1. 
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Proof 

In  view  of  Lemma  2.1,  it  suffices  to  note  that  equality  holds  in  (2.7)  for  the  particular 
optimal  control  uf  found  in  Theorem  2.1.  ■ 

Corollary  2.3 

For  any  admissible  control  u  with  corresponding  trajectory  x  ,  u  is  optimal  if  and  only  if, 
for  each  (r>0,  if  z  ( A  ,cr)  is  maximized  over  A  C.N  —d  by  a  set  R  C.N  —d  then 

(a)  xjv-j?  (°r)  =  0.  (b)  The  net  flow  out  of  R  under  u  (t )  is  at  capacity  for  a.e.  t  with 
0<r  ^ a . 


Proof 

By  the  definition  of  R  and  Corollary  2.2,  u  is  optimal  if  and  only  if  Xy((r)  =  z(R  ,cr). 
Now  by  the  state  equation  and  the  definition  of  z  , 


x_v(c)  z  (R  ,o0  =  xN_R  (<r)  +  J  [Cgjr-g  —  uR  ^ (t  )  +  uN #  (r  )]dr . 


Since  the  integrand  is  nonnegative,  the  right  side  of  this  equation  is  zero  if  and  only  if  condi¬ 
tions  (a)  and  (b)  hold.  ■ 


Remarks 


(1)  Using  the  theory  of  necessary  conditions  for  optimal  control  subject  to  state  con¬ 
straints  a  necessary  and  sufficient  condition  for  the  optimality  of  a  control  u  is  given  by  Moss 
and  Segall  (1982).  For  the  problem  P  the  conditions  include  the  existence  of  costate  functions 
X  =  [X,  it ):  i  6  A/  ]  which  satisfy  certain  backwards  differential  conditions.  As  noted  by  Moss 
and  Segall.  the  costate  functions  are  often  far  from  unique.  The  difficulty  of  applying  the 
optimality  conditions  appears  to  be  the  computation  of  X.  However,  it  is  not  difficult  to  verify 
that  X  defined  by 


are  suitable  costate  variables  for  the  optimal  control  uf  over  the  interval  [0.  cK  ].  Thus,  one 
consequence  of  our  algorithms  given  below  for  solving  problem  P0  is  that  a  particular  costate 
function  for  problem  P  can  be  readily  computed. 

(2)  In  this  chapter,  we  consider  only  a  single-commodity  routing  problem.  A  more  gen¬ 
eral  formulation  allows  multiple  commodities  with  different  destination  nodes  for  different 
commodities  [Segall,  1977].  Unfortunately,  for  this  formulation,  there  need  not  exist  an 
optimal  control  which  is  a  relaxation  of  a  constant  control  (see  examples  in  [Moss.  1977]). 
Thus,  our  algorithms  based  on  flow  relaxation  do  not  directly  generalize  to  the  multiple- 
destination  problem.  However,  it  may  be  possible  to  compute  some  choice  of  costate  variables 
in  an  efficient  way  (see  first  remark  above)  even  for  the  multiple-destination  problem,  thereby 
achieving  the  same  effect.  This  possibility  warrants  further  investigation. 

2.5.  A  Decentralized  Algorithm  for  Computing  an  Optimal  Control 

A  decentralized  algorithm  is  given  in  this  section  for  solving  problem  P0.  By  Theorem  2.1, 
this  algorithm  yields  an  optimal  control  for  the  routing  problem  P.  Since  we  will  no  longer  con¬ 
sider  state  trajectories,  the  notation  x  -  (x, :  i  €N  )  will  be  used  in  place  of  x  (0).  The  con¬ 
straints  0^  /  for  the  problem  P0  are  rather  simple,  but  the  other  constraints 

Bf  +r  ^  0  (2.13) 

can  cause  complication.  However,  if  we  assume  that 

Xj  >0  for  all  i  ,  (2.14) 

then  the  constraints  (2.13)  will  not  be  active  since  the  term  corresponding  to  node  i  in  the  cost 


22 


J0(f  )  tends  to  positive  infinity  as  Bf ,  +  r,  tends  to  zero.  Thus,  for  many  iterative  descent 
algorithms,  the  constraint  (2.13)  can  be  ignored  under  the  assumption  (2.14).  as  long  as  an  ini¬ 
tial  guess  for  /  satisfies  (2.13).  The  problem  P0  then  becomes  a  "simply  constrained  problem" 
in  the  terminology  of  Bertsekas.  Various  methods  (both  old  and  new)  for  such  problems  are 
discussed  by  Bertsekas  (1982). 

We  shall  describe  one  of  the  most  simple  methods  for  solving  problem  P0  under  the 
assumption  (2.14),  namely  gradient  descent  with  bending  and  constant  step  size  multiplier 
[McCormick.  1969].  As  noted  below,  the  solution  to  problem  P0  in  general  can  then  be  obtained 
by  a  straightforward  penalization  method. 

INITIALIZATION.  Find  a  feasible  vector  /  0  and  let  J0(f  °)=J  0.  Choose  M  to  be  larger  than 
the  maximum  eigenvalue  of  the  positive  semidefinite  Hessian  matrix  of  second  derivative 
V2/<i(  /  )  for  all  feasible  /  such  that  /0(/  )^7 o-  Set  T)=l/M . 

RECURSION  ( /  n  given  ).  Step  RJ .  Set  t£=Q  and  compute 
r,"  =  x,  /{Bf  "  +r, )  for  i*d. 

Step  R2.  Compute 

/V1  =  /.;  -  ^[(t;)2  -  (t/*)2] 

and  set 

ozfrs'sci,' 

/V^o. 

Theorem  2.2 


J  0(f  n  )  is  nonincreasing  in  n  and  in  fact 


/0(/"+i)^/o(/")-  I  I  A, 2  p 

<*.*  )€I 

where 

An  =  f  n  +1  __  f  n 

ik  “  /  ik  J  ik 

The  sequence  /  n  has  at  least  one  limit  point.  Any  limit  point  /  *  of  the  sequence  f  n  is  a  solu¬ 
tion  to  problem  P0  and 

lim/oC/"  )  =  /(>(/■)• 

n  -*oo 


Proof 

The  theorem  is  a  special  case  of  a  general  theorem  on  projected-gradient  descent  methods 
[Goldstein.  1964],  [Luenberger.  1973].  [Knowles.  1981].  See  Segall  (1979)  for  more  details  and 
an  application  in  a  setting  similar  to  ours.  ■ 

Remarks 

(1)  If  x  is  an  initial  state  such  that  =0  for  some  node  i  other  than  d  .  then  the  constraint 
Bf ,  +r,  ^0  can  become  active.  This  constraint  can  be  removed  if  the  term  dt(f  ),  defined  by 

+oo,  Bf  i  +r,  >0. 
dSf  }  -  o.  otherwise, 

is  added  to  the  cost  functional.  For  small  €  the  function  d,  is  well  approximated  by  the  func¬ 
tion 


i€2/(-5/,-rf )+. 


This  is  just  the  term  corresponding  to  node  i  in  the  sum  defining  J  §  (/  ),  where  J  $  is  defined 
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in  the  same  way  as  J0  but  with  the  modified  initial  state  x *  defined  by 

x  *  =  maxU.X/ ). 

For  each  e>0  the  recursive  procedure  given  above  can  be  used  to  compute  a  solution  f  *  to  the 
problem 

min  USif  ):0 </  £f  +r  <0}. 

If  J S  if  *)  is  finite  for  some  €  >  0  then  any  limit  point  f*  of  /  *  as  €  tends  to  zero  is  a  solution 
to  problem  P0  and  the  optimal  costs  /<>(/*)  decrease  to  J 0i  f  ' )  as  €  decreases  to  zero.  Unfor¬ 
tunately.  if  €  is  small  then  solving  problem  P*  by  the  iterative  method  above  will  require  a 
very  small  step  size  multiplier  tj  and  so  convergence  may  be  very  slow.  Alternative  decentral¬ 
ized  methods  involving  Lagrange  multipliers  which  directly  account  for  zero  initial  states  are 
being  investigated. 

(2)  The  recursive  step  is  readily  implemented  in  a  decentralized  fashion:  The  numbers  r/* 
are  determined  solely  by  the  initial  state  of  node  i  and  by  the  flows  /  *  for  links  leading  either 
into  or  out  of  node  i.  In  turn,  the  updated  flow  /5+1  depends  only  on  f  &  ,r".  and  T*. 

(3)  A  choice  for  M  in  the  initialization  step  is  not  difficult  to  find.  However,  a  fundamen¬ 
tal  limitation  of  the  method  is  that  the  multiplier  7}  required  to  guarantee  descent  will  be  quite 
small.  Presumably,  larger  values  of  tj  will  work  in  practice  and  achieve  a  satisfactory  linear 
convergence  rate,  but  this  requires  numerical  experimentation. 

A  summary  of  some  variable  step  size  rules  is  given  by  Bertsekas  (1982).  Such  rules  can  be 
used  to  identify  those  constraints  from  0^  f  which  are  active  in  a  finite  number  of  steps. 
More  sophisticated  choices  of  descent  directions  based  on  second  derivatives  of  J0  can  be  used 
to  obtain  superlinear  convergence  [Bertsekas].  Typically  the  Hessian  matrix  of  J0  is  not  diago- 


nal.  so  it  is  not  clear  how  the  more  sophisticated  algorithms  can  be  implemented  in  a  decentral¬ 
ized  fashion. 

(4)  The  recursion  step  of  the  method  described  above  provides  an  appealing  dynamic 
decentralized  routing  policy  for  time-varying  and  possibly  stochastic  inputs.  The  recursion 
step  would  be  performed  periodically  in  real  time,  and  the  constants  x,  would  be  replaced  by 
the  amount  of  traffic  at  node  i  at  the  time  of  the  recursion  plus  perhaps  a  correction  term  to 
reflect  expectations  of  future  traffic  input  fluctuations.  It  would  be  interesting  to  establish 
optimal  results  along  these  lines.  See  oegall  (1977)  for  a  discussion  of  the  routing  problem  with 
stochastic  inputs. 

2.6.  An  O  ( I N  1 4)  Algorithm  for  Computing  an  Optimal  Control 

An  efficient  algorithm  is  given  in  this  section  for  finding  a  solution  /  to  problem  P0  for  a 
given  network  N  =  {N  X  .C  .d  )  and  demand  (x  r  ).  According  to  Theorem  2.1  this  produces  an 
optimal  control  for  problem  P  by  flow  relaxation. 

The  algorithm  is  suggested  by  Lemma  2.3  of  Section  2.4.  If  for  some  <r>  0  the  set 
A  CiV  —d  maximized  z  (A  .<r)  then  condition  (a)  of  Lemma  2.3  specifies  the  optimal  flow  on 
links  connecting  A  and  N  —  A  .  Then  the  original  problem  can  be  decomposed  into  two  indepen¬ 
dent  problems  —  finding  /  for  links  connecting  nodes  in  A  and  finding  /  for  links  connecting 
nodes  in  N  —A  .  To  find  the  set  A  .  we  will  make  use  of  the  famous  max-flow  min-cut  theorem. 

Given  a  network  N  =  (N  ,L  ,C  j  4  )  with  source  node  s  and  destination  node  d  ,  a 
max-flow  is  a  vector  g  =  (g, :  l  €Z, )  which  is  a  solution  to 

maxig, y_s  -g*_, _,:<><$  ^C.giJV_,  -fcy-i.,  =0,  itN-s-d). 

We  then  write  g  -  MAXFLCXN).  An  s  —d  cut  for  N  is  a  set  A  of  nodes  with  s  €A  and  d  $A  , 
and  the  capacity  of  the  cut  is  CA  iV_^  .  An  s  —  d  cut  is  called  a  min  s  —d  cut  if  its  capacity  is 
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minimum  among  the  capacities  of  all  s  —d  cuts.  By  the  max-flow  min-cut  theorem  [Ford  and 
Fulkerson.  1962],  the  value  gs  jj-s  —  gjv-sj  of  a  max-flow  g  is  equal  to  the  capacity  of  a  min 
s  —d  cut  A  .  Moreover, 

gA  ,.V  -.4  =  Ca  „V  -A  .  g.V  -A  A  ~ 

Many  efficient  algorithms  exist  for  computing  max-flow  and  an  associated  min-cut  for  a 
network.  A  particularly  efficient  yet  simple  algorithm  is  given  in  Malhotra  et  al.  (1978).  This 
algorithm  is  a  simplified  version  of  successively  improved  algorithms  by  Ford  and  Fulkerson. 
Edmonds  and  Karp,  Kinits,  Karzanov,  and  others.  For  an  excellent  introduction  and  references, 
see  Papadimitriou  and  Steiglitz  (1982).  The  algorithm  of  Malhotra  et  al.  requires  0(\N  I3) 
computations. 

A  max-flow  algorithm  can  be  used  to  find  a  set  A  maximizing  z  (A  ,<r)  for  a  given  <r>0, 
network  N  =  (N  ,L  ,C  .d  )  and  demand  (x  r  )  in  the  following  way.  Define  a  single-source. 

a 

single-destination  network  N  by 

N  =  (iU  NM.Cj.d)  (2.15) 


>  v 


^7 
>»  , 
*  Ik 


where  s  is  a  source  node  with  s  %N  , 


L  =L  U{(s.j):  j£N-d). 

and 


Cn 


C.J, 

Xj  /<T  +  rj  , 


(i  ,j  )€£  . 
i  -s  ,j  €iV—  d. 


(2.16) 


(2.17) 
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Lemma  2.4 

For  fixed  cr>0,  a  subset  A  of  N  —  d  maximizes  z  ( A  ,cr)  over  all  subsets  A  of  N  —  d  if  and 

a 

only  if  j  UA  is  a  min  s  —d  cut  for  N. 

Proof 

The  capacity  of  any  s  —d  cuts  s  U<4  for  N  (where  A  C.N  —d  )  is 


's  im  w 


-A  =  C$W-A  +  C. 


A  W -A 


=  I 

jiN-A 


+  r, 


+  Caw -a 


XN-A 


+  rN-A  +  C 


AW -A 


XN  , 


z  (A . <r ) 


which  is  a  strictly  decreasing  function  of  z  ( A  ,c r). 


The  algorithm  will  now  be  presented. 


FUNCTION  OPTFLO 


Input :  A  single-destination  network  N  =  (N  X  ,C  4  )  with  demand  (x  r  ). 
Output :  An  optimal  flow  for  problem  P0  if  one  exists  with  finite  cost. 


Step  1.  if  xN  9*0  then  go  to  Step  2: 

else  define  N°  -  ( N°JL0.C°.s  ,d  )  by 

N°  =  s UN , 

L°=L  U|(j  ):  jtN-d). 


r  .?= 


c,j  •  (» .j  )€Z,  , 

,  i  =  s  .j  GN  —d 
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Step  2. 
Step  3. 

Step  4. 

Step  5. 


if  {s  }  is  a  min  5  -d  cut  for  N°  then  let  OPTFLO,  =  MAXFLO,  for  (N°)  for  l  £L 

and  return; 

else  flag  that  J‘0  is  infinite  and  return. 

if  Cfit-d  4  ~fN-4  ^  0  then  flag  that  J'0  is  infinite  and  return. 


Let 

<J  =  xN /(Cv_rf  4  —tn-4  ) 

and  define  the  network  N  by  equations  (2.15)  -  (2.17).  Let  g  «  MAXFLO  (N)  and 
let  £  CA/  — d  be  such  that  sUi  is  a  min  s  —d  cut  for  N. 

if  {j  }  is  a  min  s  —d  cut  for  N 

then  let  OPTFLO,  =  gt  for  l  €L  and  return. 


Let  N'  =  (N  '.L ' ,C ' .d  ’)  with  demand  ( x  'r ').  where 
N'  =  R  U d\ 

L '  =  {(*  ,j  )€£  :  i  ,j  6 R  }  U{(i  .d '):  iCR). 


<v- 


i  .j  €R . 
Cijf-R’  i£R.j=d'. 


ri  -  ri .  x‘  =  Xi .  i  6 R. 


Let  N"  =  (AT  ,Lm  .C"  xt )  with  demand  (x*  r"  ),  where 


AT  =  N-R ,  <T  =  d. 

LT  =  {(i  ./)€/.:  i ./  SN-R). 

Cij  =Cjj  for  (i  .j  )  in  L"  . 

r"  =  r,  +  Cj? , .  x,"  =  Xj ,  i£N—R. 


r.  1 


Step  6.  Let  /  *  =  OPTFLO(N>  '.x ').  / "  =OPTFLO(N"  r  a  ).  and 


OPTFLO,,  =  /,, 


•-'ij-  i  €R .  j  €N—R  , 

0,  i€N-R.j€R. 
i.j€N-R: 

return. 

Theorem  13  (OPTFLO  works) 

Given  a  network  N  =  (N  L  .C  ,d  )  and  demand  (x  x  ).  if  the  minimum  cost  J’0  for  problem 
P0  is  Unite,  then  /  =  OFT  FLO  (Njc  r  )  is  an  optimal  flow.  Otherwise  OPTFLO  flags  that  J‘0  is 
infinite. 

Proof 

The  proof  is  by  induction  on  the  number  of  nodes  in  the  network.  If  there  is  only  one 
node  in  the  network  (namely  d  )  then  OPTFLO  works  as  claimed;  that  is.  return  is  hit  in  Step 
1.  Now  suppose  that  N  is  a  given  network  with  demand  Cx  r  )  and  that  the  theorem  has  been 
verified  for  all  networks  with  strictly  fewer  nodes  than  are  in  N. 

We  now  assume  that  Jq  is  finite  for  N  and  (x  x  );  we  prove  that  /  =  OPTFLO  is  an 
optimal  flow  for  problem  P0.  Verification  that  OPTFLO  flags  that  J’0  is  infinite  when  it  is  sup¬ 
posed  to  can  then  be  readily  carried  out  by  the  reader. 

First,  consider  the  case  in  which  xN  =0.  Then  by  definition,  a  flow  /  is  optimal  for  prob¬ 
lem  P0  if  and  only  if  it  is  feasible  for  problem  P0  (and  the  minimum  value  is  zero).  Since  a 
feasible  flow  is  returned  in  Step  1  if  such  a  flow  exists,  the  induction  step  is  proved  in  the  case 
x  v  =0. 


Thus,  we  now  suppose  that  xN  >0.  Then  the  condition  of  Step  2  implies  that  there  is  no 
feasible  flow  for  problem  P0.  We  have  assumed,  however,  that  /q  is  finite  so  that  the  algo¬ 
rithm  will  not  return  in  Step  2.  Then  Step  3  is  executed  and  defines  <r>0.  g  ,  and  R . 


Now  suppose  that  {j  }  is  a  min  s  —d  cut  for  N.  Then,  using  the  fact  that  the  net  flow  for  g 
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into  any  node  t  €N  —d  is  zero,  we  have 

g,  ji-i  ~  g.v-i =  gs,  =  4  =  x,  /O'  +  r,  for  iZN-d  . 

so  for  each  node  t  €  A 1  —d  . 

x,  /(g,  ..v  ~  g.v  — i -  r,  )  =  cr.  x,  >  0. 

Now  by  the  definitions  of  C  and  <x. 

n  _  <r'  "'•i  ,  _  Xfj i 

V'S  .X-d  ~  JL  —  +  ri - — -  +  riV  —  '-'TV  -t/  ^  ■ 

it  X -d  V  a 

Therefore  j  U  A/  —  d  is  also  a  min  s  —d  cut  for  N,  and  so 


(2.18) 


(2.19) 


gid  =  Cid  •  gd,  =  0.  for  i  €N  —d. 

This  fact  and  equation  (2.18)  show  that  if  T  is  defined  by  T;  =  <r  for  all  i  £N  —d  then  (g  .t) 
satisfies  the  optimality  conditions  for  problem  P0.  Therefore  OPTFLO  returns  and  optimal  flow 
if  1  s  }  is  a  min  s  —d  cut. 

Now  suppose  that  {5  }  is  not  a  min  s  —d  cut  for  N.  Since  by  equation  (2.19)  the  s  —d  cuts 
{j }  and  s  UN— d  have  the  same  capacity,  we  conclude  that  j  UA^-rf  is  also  not  a  min  s—d 
cut.  The  set  R  computed  in  Step  2  is  therefore  neither  empty  nor  equal  to  N  —  d  .  and  so  the 
networks  N'  and  N"  each  have  strictly  fewer  nodes  than  the  original  network  N.  Our  induc¬ 
tion  hypothesis  implies  that  the  theorem  is  valid  for  OPTFLO  applied  to  (N'.x'.rO  and  to 
(N"  .r"  ). 

Let  / o(/')  and  Jq  (/" )  be  defined  for  these  networks  and  demands  in  the  same  way  that 
J0(f  )  is  defined  for  N  and  (x  ,r  ).  and  let  Pq  and  Pq  denote  the  corresponding  minimization 
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problems.  We  claim  that  the  minimum  costs  for  problems  Po  and  Po  are  each  finite.  To  prove 
this,  let  /  be  any  solution  to  problem  P0.  Then  by  the  definition  of  R  ,  Lemma  2.3,  and  Lemma 
2.4. 


/  R  _V  -R  ~  CR  j 

;-R  •  f  .V  -R  ,R  —  0. 

(2.20) 

Thus  if  /  ’  and  /  " 

are  defined  by 

T‘J- 

f  l>  C,  ,.V  -R 

i  .j  €R . 

(2.21) 

,  i  €R  .  j  —d. ' . 

1^ 

II 

* 

1^ 

i.j€N-R. 

(2.22) 

then  /  '  and  /  "  are  feasible  for  problems  Pq  and  P£  .  respectively,  and 

Jo  =  /«(/  )  =  J'iT ')  +  /o(/"  ).  (.2.23) 

and  so.  since  J"o  by  assumption  is  finite,  problems  Pq  and  Po  have  finite  minimum  costs  as  well. 

Since  the  theorem  is  valid  for  OPTFLO  applied  to  (N'jr '  x ')  and  to  (N*  jc*  x"  ).  and  since 
the  minimum  costs  for  Pq  and  Po  are  finite,  when  OPTFLO  is  executed  in  Step  6  it  will  return 
optimal  solutions  /  ’  and  /*  for  problems  Po  and  Pq  .  respectively.  Since  J  <ff')  ^  J  o(f  ’)  and 
J"o(f)  ^  Joif  )•  equation  (2.23)  yields  that 

Jo  >  J&f)  +  JZ(ri  (2.24) 

On  the  other  hand  it  is  immediate  from  the  definitions  of  /  (in  Step  6).  /  '.  and  /"  that  /  is 
feasible  for  problem  P0  and  that 


J(lf)  +  Jo  if"  )  ^  Joif  )■ 
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In  view  of  inequality  (2.24),  J'0  =  /<>(/  ).  and  so  /  is  indeed  optimal  for  problem  P0.  ■ 

Computational  Complexity  of  OPTFLO 

The  algorithm  OPTFLO  is  recursive.  When  it  is  first  called,  either  execution  terminates  in 
Steps  1.  2.  or  4.  after  at  most  one  execution  of  MAXFLO.  or  else  MAXFLO  is  executed  once  and 
then  OPTFLO  is  applied  to  the  two  networks  N’  and  N"  in  Step  6.  In  the  latter  case,  the  set  of 
nondestination  nodes  N  —d  is  partitioned  into  two  nonempty  subsets  which  become  the  sets  of 
nondestination  nodes  in  N '  and  V  .  Therefore,  taking  all  recursions  into  account.  OPTFLO 
(and  hence  MAXFLO)  is  called  at  most  2INI-1  times. 

We  assume  that  an  algorithm  MAXFLO  is  used  which  requires  O  ( I  N  I  3)  computations 
[Malhotra  et  al.].  Therefore,  the  number  of  computations  needed  for  all  executions  of  MAXFLO 
when  OPTFLO  is  called  is  at  most  O  ( I  N  1 4).  Since  this  clearly  dominates  the  number  of  all 
other  computations  needed  to  execute  OPTFLO.  we  conclude  that  OPTFLO  requires  at  most 
O  (  I N  I  4)  computations.  Similarly,  if  other  existing  MAXFLO  algorithms  which  require 
O  ( I  N  I  2  I  L  I  *)  computations  or  O  ( I  N  I  I  L  I  log  I L  I  )  computations  were  used  (see  [Papadim- 
itriou  and  Steiglitz]).  then  OPTFLO  would  require  0(\N  l3IZ,  I1*)  or  0(liV  l2l  L  I  log  I L  I) 
computations,  respectively. 

2.7.  Solution  By  Search  For  Successive  Bottlenecks 

An  algorithm  is  given  in  this  section  for  solving  problem  P0  which  is  based  on  a  search  for 
bottlenecks.  By  a  bottleneck  we  mean  a  subset  of  nodes  R  such  that,  based  on  the  initial  traffic 
xR  (0)  in  R  .  on  the  input  rate  rR  ,  and  on  the  total  capacity  CR  ,  R  will  take  as  long  to  com¬ 
pletely  empty  as  any  other  subset  of  nodes.  The  key  motivation  of  the  algorithm  is  that  the 
optimal  initial  net  flow  out  of  such  a  bottleneck  set  R  must  be  at  capacity.  Thus,  once  R  is 
found,  the  network  decouples:  flows  must  be  found  within  R  so  that  all  nodes  within  R  empty 
simultaneously  and,  independently,  flows  must  be  found  within  N—  R  is  the  same  as  the 


33 


original  problem  for  a  new  network  and  demand  with  node  set  N —R  .  Hence  a  bottleneck 
within  the  new  network  must  be  found.  The  process  continues  until  the  whole  network  is  par¬ 
titioned  into  subsets  of  nodes  which  are  successive  bottlenecks.  The  precise  algorithm  is 
described  below. 

Unfortunately  the  bottleneck  approach  does  not.  as  yet.  appear  to  be  computationally  com¬ 
petitive  with  OPTFLO.  The  main  reasons  for  presenting  the  algorithm  are  twofold.  First,  for  a 
given  network  it  may  be  possible  to  spot  successive  bottlenecks  by  inspection,  and  second,  the 
algorithm  is  similar  to  the  algorithm  OPTFLO  but  is  conceptually  simpler. 


FUNCTION  BOTFLO. 


Input  and  Output :  Identical  to  those  of  OPTFLO. 
Step  1.  Identical  to  Step  1  of  OPTFLO. 

Step  2  For  each  A  C N —d  let 


q(A  )  = 


XA  /(Ca  -V  -A  ~  r  A  ).  C-A  JV  -A  ~  f  A  >  0 
0.  XA  =  CA  v  -a  ~  r A  —  0, 

+00  otherwise. 


Let  R  CN  —  d  be  a  set  which  maximizes  q  {A  )  over  all  A  CN  —  d  . 


( Comment :  R  ^  0  and  q  (R  )  >  0.) 


Step  3  If  q  (R  )  =  +  oo  flag  that  J'o  is  infinite  and  return. 

Step  4  Let  N  =  (iV  ,L  .C  .s  ,d  ')  be  defined  by 

N  =  R  Ur  U d\ 

L  =  l(i  .j  )€L:  i  eR.  j  6/?}U{(r  ,j  ):  j  €/?  }U{U  ,d‘):  i€R  ). 
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c>j-  iGR.jGR. 

C,j  -  Xj  lq  (R  )  +  r,  i  =  s . j  GR . 

Cijr-R-  i  GR .  j  —  d  , 

i 

and  let  f  =  MAXFLO(N). 

Step 5.  Let  (N*  ,x“  r  )  be  defined  as  in  Step  5  of  OPTFLO.  let  -  BOTFLOCN*  ,r"  .x" ). 

and  let 


BOTFLO  = 


f.j. 

i£R.j€R. 

C,J. 

i  GR.  J€N-R. 

0. 

i  GJV-R,  j  GR  , 

f 

J  IJ  • 

i  GN-R.  j  GN-R 

return. 


Theorem  2.4  (BOTFLO  works.) 

Given  a  network  N  =  (A  X  ,C  .d  )  and  demand  (x  x  ),  if  the  minimum  cost  J'0  for  problem 
P0  is  finite  then  /  -  OPTFLO(Njc  x  )  is  an  optimal  flow.  Otherwise  OPTFLO  flags  that  J‘0  is 
infinite. 

Proof 

The  proof  will  only  be  sketched  since  it  closely  parallels  the  proof  of  Theorem  2.3.  We 
argue  by  induction  on  I  N  I  and  assume  that  J‘0  is  finite.  If  there  is  only  one  node  in  N  of  if 
x  v  =  0  then  BOTFLO  returns  an  optimal  flow.  We  thus  suppose  that  xN  >  0.  Then  q  {R  )  >  0 
and  R  is  not  empty.  If  q  (/?)>+  oo  then  no  feasible  solution  exists  to  problem  P0  so  that 
return  will  not  be  executed  in  Step  3  and  Step  4  will  be  executed.  Since  N  —R  is  strictly 
smaller  than  N ,  our  induction  hypothesis  implies  that  BOTFLO  is  valid  when  applied  to  the 
network  and  demand  (N"  jc"  x" ). 

Let  <r  =  q  (R  ).  Now  q  (A  )  ^  cr  for  A  C N  —d  with  equality  if  A  =  R  .  Therefore  both 
A  =  R  and  A  =  0  maximize  z  ( A  .<7)  over  subsets  A  of  N  —d  .  By  Lemma  2.3  s  U  A  and  (j  } 
are  each  min  s  —d  '  cuts  for  N.  so 
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f  s.i  —  r,  +  *i/o\  i  £R  . 


I 


and 


(2.26) 


/  i.d'  ~  Cijt-R  •  * 


Since  the  net  flow  into  each  node  i  in  R  is  zero  for  /  . 


(2.27) 


f  s  ,i  +  f  R  — i .(  f  I  JP  Ud  i  —  0. 


(2.28) 


t 


Now  let  (N‘.x  'r ')  and  problem  P0’  be  defined  as  in  the  algorithm  OPTFLO  and  let  /  ‘  be  defined 
by  f  -  f  i  for  l  €L'.  Then  equations  (2.26)  and  (2.27)  imply  that 

*,/(/  ~  i^R.x^O.  (2.29) 

Conditions  (2.27)  and  (2.29)  show  that  /  '  and  the  vector  r  defined  by  r,  =  tr  for  i  6 R  satisfy 
the  optimality  conditions,  implying  that  /  *  is  optimal  for  problem  P0’. 

Another  consequence  of  the  fact  that  R  maximizes  z  (A  ,<r)  is  that,  by  Lemma  2.3.  if  7"  is 
any  solution  to  problem  P<>  then  equation  (2.20)  holds.  We  conclude  as  in  the  proof  of  Theorem 
2.3  that  the  minimum  delay  for  problem  P0"  is  finite.  The  final  step  of  the  proof  is  the  same  as 
that  of  Theorem  2.3.  ■ 


t 

► 

►  * 

►  ■ 

r 


The  dominant  factor  in  the  computational  complexity  of  BOTFLO  is  finding  the  bottleneck 
set  R  in  Step  2.  We  shall  show  next  that  R  can  be  found  by  solving  the  linear  programming 
problem  (LP)  defined  by 
min  cr 


-h  <0. 
h  -Car  ^  0. 
x  +  rcr  +  Bh  <  0. 

— cr  <  0. 


(LP) 


f 


».-  ».•  <■.  v'.O  ^  •  • — ■- 
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If  (<r.h  )  is  a  solution  to  LP  then  cr  is  the  minimum  amount  of  time  required  to  evacuate  the 
network  and,  if  <r  >  0,  uf  .  where  /  =  h  la,  is  an  admissible  control  which  evacuates  the  net¬ 
work  by  time  <r.  Problem  LP  plays  an  important  role  in  the  work  of  Shats  and  Segall  (1980). 

Proposition  2A 

Suppose  (crJi  )  is  a  solution  to  problem  LP  for  some  network  N  and  demand  ( x  r )  with 
x  v  >  0.  Define  /  —  h  la  and  define 

R  =  {i  €N-d  : 

no  path  from  i  to  d  |J  {k  :xk  +  a(.Bf  k  +  rk  )  <  0}  is  unsaturated  for  /  } 


Then 


xR  >  0, 

(2.30) 

z  (R  .o’)  =  max) z  ( A  .cr):  A  CN—d  }  =  0. 

(2.31) 

fsjv-g  =  Cr  jv  -r  .  f  n-r  jt  =0. 

(2.32) 

x*  +  a(Bf  j  +  rj  )  =  0.  i  €R. 

(2.33) 

Proof 

The  Kuhn-Tucker  optimality  conditions  for  problem 

LP  imply  that  there  exist  vectors 

a.X.fi  with  -  0  such  that 

'-T 

\TC  -nTr  +1=0. 

(2.34) 

-<*ij  +  My  ~  Mi  =0.  (i )€Z, , 

(2.35) 

and 


aij 

>  0, 

<*ij  =0. 

unless  f  ij  =0, 

ifc  * 

Kj 

£  0. 

Kj  =  o. 

unless  f  i,  -  Cij , 

(2.36)  - 

M, 

>  0. 

o' 

ll 

€ 

unless*;  +  c r(r,  +  Bf , )  =  0. 

(2.36) 
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r 


w". 


c 

* 


Let  /?'  =  {*:  /x.i  >0}.  If  xg-  —  0  then 

0=  Z  ftWi  +r.)=  £  f  +  nrr  Z  -  krC  +  fSr  . 

/€iV-rf  (ij)€l 

which  contradicts  conditions  (2.34),  so  we  have  proved  that  >  0.  Conditions  (2.35)  and 
(2.36)  imply  that  R '  C  R  .  so  (2.30)  is  established.  Equations  (2.32)  and  (2.33)  are  immediate 
from  the  definition  of  ^  .  Finally,  for  A  CN—d  . 

Z  u  ,cr)  =  X.4  -  cKQ  -  rA ) 

=  AJi-A  ~  Ca  y-A  ~  f  N-A  y  )  +  [*A  +  °i-f  N-A  A  ~  /  A  y  -A  +  rA  )!■ 

which  is  nonpositive  and  equal  to  zero  if  A  =  R  (by  (2.32)  and  (2.33)).  which  proves  equation 
(2.31).  ■ 

Corollary  2.4 

Theorem  2.4  remains  true  if  Steps  3  and  4  of  BOTFLO  are  replaced  by  the  following: 

Solve  LP.  If  no  feasible  vector  for  LP  exists  then  flag  that  J‘0  is  infinite  and  return. 
Else  let  /  —  h  Ar  and  compute  R  defined  in  Proposition  2.4  (with  /  =  /  ). 


Proof 

In  view  of  Proposition  2.4.  the  proof  of  the  corollary  is  a  straightforward  modification  of 
the  proof  of  Theorem  2.4.  ■ 

Remarks 

(1)  In  Section  VI  we  solved  the  problem  of  maximizing  z(A  ,<r)  over  A  C N —d  by  recog¬ 
nizing  that  the  solution  R  is  min-cut  for  a  max-flow  problem  (which  is  linear).  It  is  well 
known  that  a  min-cut  for  a  max-flow  problem  corresponds  to  an  optimal  choice  of  dual  vari¬ 
ables  (see  [Ford  and  Fulkerson,  p.29]).  Therefore,  the  problem  of  maximizing  z(A  ,<r)  over 


•  •  »  m  ■  p  ■  a  •  a  *  ■«*»  *  k  *  >  *  •  *  •  "  •  "  « 

*  ,  «  *  .  '  v'  ’  .  '  .  *  v  i  *  -  V’  *•  *  * 


A, 


A  C.N  —d  is  dual  to  a  max-flow  problem.  Similarly,  as  indicated  in  Proposition  2.4.  the  prob¬ 
lem  of  maximizing  q  (A  )  over  A  C.N  —d  corresponds  to  the  dual  of  the  linear  problem  LP. 

(2)  By  arguments  similar  to  those  at  the  end  of  the  preceding  section  it  can  be  shown  that 
if  R  in  Step  2  of  BOTFLO  is  chosen  maximally  (for  example,  if  R  is  computed  by  the  method 
of  Corollary  2.4)  then  BOTFLO  will  be  called  at  most  K  +  1  times  during  complete  execution, 
where  K  is  defined  in  the  previous  section  and  satisfies  K  ^  \  N  —d  I .  Thus,  if  special  proper¬ 
ties  of  the  problem  LP  can  be  exploited  to  generate  an  algorithm  with  about  the  same  complex¬ 
ity  as  efficient  max-flow  algorithms,  then  BOTFLO  would  be  computationally  competitive  with 
OPTFLO. 

(3)  We  will  describe  a  possibly  more  efficient  method  to  find  the  bottleneck  set  R  needed  in 
Step  2  of  BOTFLO.  Modify  OPTFLO  so  that  N"  in  Step  5  and  /"in  Step  6  are  no  longer  com¬ 
puted.  Then  when  OPTFLO  is  executed  (assume  that  /o  is  finite)  it  either  returns  after  one 
execution  of  MAXFLO  or  else  OPTFLO  is  applied  again  to  the  network  N',  which  has  strictly 
fewer  nodes  than  N.  OPTFLO  will  continue  to  generate  successively  smaller  networks  until  at 
some  level  of  recursion  a  network  is  found  for  which  OPTFLO  terminates  in  Step  1  or  Step  4. 
The  set  of  nondestination  nodes  in  this  network  is  a  bottleneck  set  R  needed  for  Step  2  of  BOT¬ 
FLO.  While  this  may  be  a  relatively  efficient  method  for  computing  a  bottleneck  set  R  (it 
requires  at  most  O  ( I N  1 4)  computations),  if  it  is  used  in  Step  2  of  BOTFLO  then  BOTFLO  will 
clearly  require  more  computations  (although  at  most  O  ( IN  1 5))  than  are  required  by  OPTFLO. 


3.  MINIMUM-DELAY  DYNAMIC  ROUTING:  PIECEWISE-CONSTANT  CAPACITIES 


3.1.  Introduction 

In  this  chapter  we  consider  a  single-source  single-sink  network  where  the  link  flows  are 
real-valued  measurable  functions  defined  on  a  time  interval  (0,  T  ]  and  where  storage  is  allowed 
at  the  nodes.  The  flow  on  each  link  (i  ,j  )  at  each  time  t  represents  the  net  rate  at  which  some 
commodity  moves  from  node  i  to  node  j  at  time  t .  and  thus  the  storage  at  each  node  changes 
at  a  rate  equal  to  the  sum  of  the  flows  on  all  links  entering  that  node.  The  flow  on  each  link  is 
constrained  to  be  dominated  by  a  piecewise-constant  capacity  assigned  to  that  link,  and  the 
storage  at  each  node  is  constrained  to  be  dominated  by  a  piecewise-constant  storage  capacity 
assigned  to  that  node. 

By  scaling  time  on  each  interval  of  constant  capacities,  we  assume  without  loss  of  general¬ 
ity  that  T  is  an  integer  and  that  the  link  and  storage  capacities  change  exactly  at  the  times 

1 . T— 1.  For  example,  if  the  capacities  are  constant,  then  we  assume  T  =  l.  Note,  however. 

that  the  flows  may  change  at  any  time.  Thus,  the  model  we  consider  is  not  equivalent  to  a 
discrete-time  model. 

By  a  r-maximum  flow  we  mean  a  flow  which  maximizes  the  total  amount  of  commodity 
reaching  the  sink  before  time  r.  The  main  problem  considered  in  this  chapter  is  that  of  com¬ 
puting  a  flow  which  is  simultaneously  a  r-maximum  flow  for  each  r€(0,  T).  As  shown  in  Sec¬ 
tion  3.6.  such  a  flow  solves  a  minimum-delay  problem  of  the  type  considered  in  Chapter  2. 
However,  the  model  in  Chapter  2  is  not  as  general  since  the  capacities  there  were  assumed  to  be 
constant  in  time. 

The  main  result  of  this  chapter  is  an  algorithm  which  computes  the  desired  flow  in 
O  ( I N  1 4T*)  computations,  where  I N  I  denotes  the  number  of  nodes.  This  algorithm  first  finds 
the  breakpoints  of  the  optimal  flow  (i.e..  the  points  of  time  at  which  it  changes)  by  solving  a 
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3.  MINIMUM-DELAY  DYNAMIC  ROUTING:  PEECEWISE-CONSTANT  CAPACITIES 

3.1.  Introduction 

In  this  chapter  we  consider  a  single-source  single-smk  network  where  the  link  flows  are 
real- valued  measurable  functions  defined  on  a  time  interval  (0,  T  ]  and  w  here  storage  is  allowed 
at  the  nodes.  The  flow  on  each  link  ( i  .j  )  at  each  time  t  represents  the  net  rate  at  which  some 
commodity  moves  from  node  i  to  node  j  at  time  t .  and  thus  the  storage  at  each  node  changes 
at  a  rate  equal  to  the  sum  of  the  flows  on  all  links  entering  that  node.  The  flow  on  each  link  is 
constrained  to  be  dominated  by  a  piecewise-constant  capacity  assigned  to  that  link,  and  the 
storage  at  each  node  is  constrained  to  be  dominated  by  a  piecewise-constant  storage  capacity 
assigned  to  that  node. 

By  scaling  time  on  each  interval  of  constant  capacities,  we  assume  without  loss  of  general¬ 
ity  that  T  is  an  integer  and  that  the  link  and  storage  capacities  change  exactly  at  the  times 

1 . T  —  1.  For  example,  if  the  capacities  are  constant,  then  we  assume  T  —  1.  Note,  however. 

that  the  flows  may  change  at  any  time.  Thus,  the  model  we  consider  is  not  equivalent  to  a 
discrete-time  model. 

By  a  r-maximum  flow  we  mean  a  flow  which  maximizes  the  total  amount  of  commodity 
reaching  the  sink  before  time  r.  The  main  problem  considered  in  this  chapter  is  that  of  com¬ 
puting  a  flow  which  is  simultaneously  a  r-maximum  flow  for  each  r€(0.  T  ].  As  shown  in  Sec¬ 
tion  3.6.  such  a  flow  solves  a  minimum-delay  problem  of  the  type  considered  in  Chapter  2. 
However,  the  model  in  Chapter  2  is  not  as  general  since  the  capacities  there  were  assumed  to  be 
constant  in  time. 

The  main  result  of  this  chapter  is  an  algorithm  which  computes  the  desired  flow  in 
O  ( I N  I  *T 4)  computations,  where  I N  I  denotes  the  number  of  nodes.  This  algorithm  first  finds 
the  breakpoints  of  the  optimal  flow  (i.e..  the  points  of  time  at  which  it  changes)  by  solving  a 


series  of  static  maximum  flow  problems.  It  then  constructs  the  desired  flow  by  piecing  together 
flows  which  are  each  r-maximum  for  some  breakpoint  r.  In  the  special  case  where  the  capaci¬ 
ties  are  constant  in  time,  the  minimum-delay  problem  of  Section  3.6  is  equivalent  to  the  one 
considered  in  Chapter  2,  and  the  algorithm  computes  an  optimal  flow  in  0(  I  N  I4)  computa¬ 
tions.  as  did  the  algorithm  OPTFLO  of  Chapter  2.  Other  researchers  (see  Section  2.1)  have  pro¬ 
duced  less  efficient  algorithms  for  the  the  case  of  constant  capacities,  but  this  chapter  presents 
the  first  polynomial  algorithm  for  the  case  of  piecewise-constant  capacities. 

This  chapter  is  organized  as  follows.  In  Section  3.2  we  formulate  the  problem,  and  in  Sec¬ 
tion  3.3  we  review  the  well-known  static  max-flow  min-cut  theorem  in  the  required  form  for 
our  application.  In  Section  3.4  we  present  a  dynamic  version  of  the  max-flow  min-cut  theorem 
and  show  how  a  r-maximum  flow  can  be  computed  by  solving  a  static  maximum  flow  problem. 
We  also  derive  necessary  and  sufficient  conditions  for  a  flow  to  be  r-maximum.  In  Section  3.5 
an  algorithm  is  derived  for  computing  a  flow  which  is  r-maximum  for  all  r  6(0.  T  ].  and  in  Sec¬ 
tion  3.6  we  show  that  such  a  flow  also  solves  a  minimum-delay  problem.  Section  3.5  also  con¬ 
tains  an  interesting  characterization  (Lemma  3.6)  of  a  flow  which  is  r-maximum  for  all  r.  and 
therefore  of  a  minimum-delay  flow. 

3.2.  Problem  Formulation 

Let  N  be  a  finite  set  of  nodes,  including  a  source  s  and  a  sink  d  ,  and  fix  T  >0.  By  a  link 
we  mean  an  element  of  N2.  For  each  link  (i  ,j )  we  assign  a  nonnegative  left-continuous 
piecewise-constant  function  btj  on  (0.  T  ]  giving  its  capacity  at  each  time,  and  for  each  node 
j  €N  we  assign  a  nonnegative  lower-semicontinuous  piecewise-constant  function  aj  on  (0.  T  ] 
giving  its  storage  capacity  at  each  time.  Letting  b  —  ( bi} :  i  .j  )  and  a  -  (aj :  j  6  N  ).  we 
define  the  continuous-time  dynamic  network  D  =  (N ,  b  .  a  .  s  .  d  ). 

We  define  a  feasible  flow  for  D  to  be  an  assignment  u  =  (u,  :  i  .j  6/V  )  such  that  each  utJ  is 


a  Lebesgue-measurable  function  on  (0.  T]  and  the  following  conditions  are  satisfied  for  all 
i  €N  and  all  t  €(0.  T  ]: 

U,j  (t  )  =  —Uj,  ( t  ) 

u,j  (f )  ^  b,j  (t ) 

0  <  Xj  (t  )  £  a}  (t  ) 
where  x  =  (x; :  j  €N  ).  called  the  storage  function  due  to  u  .  is  defined  by 

t 

Xj(t)  =  f  ^Uijit^dt'.  j€N—{s.d } 

0  i  €  -V 

xs  (t  )  =  xd  (t  )  =  0. 

The  set  of  all  feasible  flows  for  D  will  be  denoted  U .  We  interpret  utJ  as  the  net  rate  at  which 
commodity  is  moved  from  node  i  to  node  j ,  and  Xj(t  )  as  the  amount  of  commodity  stored  at 
node  j  at  time  t.  Note  that  (3.1)  implies  uu  (t  )=0  and  that  (3.1),  (3.2)  imply 
—bj,  (t  (r  )^b,j  (t  ).  Also  note  that  we  are  taking  the  initial  storage  to  be  zero.  However, 
we  will  see  in  Section  3.6  that  this  restriction  does  not  prevent  us  from  allowing  nonzero  initial 
storage  in  the  minimum-delay  problem. 

As  explained  in  the  Introduction,  we  assume  without  loss  of  generality  that  T  is  a  positive 

integer  and  that  b tj  and  aj  are  constant  on  each  interval  (m— l^n  ),  m  =1 . T .  Thus,  since 

bjj  is  left-continuous,  the  value  of  b^  on  (m  —  l,m  ]  may  be  represented  by  b,j  (m  ).  and  since  aj 
is  lower-semicontinuous,  a^(m)  is  no  greater  than  the  value  of  aj  on  (m— l^n)  or  on 
(m  ,m  +1).  Note  that  since  we  are  not  requiring  a  flow  u  €£/  to  be  constant  on  each  interval 
(m  — l.m  ),  the  above  model  is  not  equivalent  to  a  discrete-time  model. 

When  sets  of  nodes  are  used  as  subscripts,  the  summation  convention  applies.  Thus 


(3.1) 

(3.2) 

(3.3) 
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uAJ(.t)=  £  £  UijU ). 

I $A  i € B 

For  each  r  €  (0,  T  ]  we  define  v  T  on  U  by 

T 

vT(u  )  =  f  uN  d  Cr )  dr. 
o 

This  represents  the  total  amount  of  commodity  reaching  the  sink  before  time  r.  By  a 
r— maximum  flow  we  mean  a  flow  u  €U  such  that  vT(u  )^vT(u ')  for  all  u'£U.  The  main 
problem  considered  in  this  chapter  is  that  of  computing  a  flow  which  is  r-maiimum  for  all 
t€(0,  T  ].  At  this  point  it  is  not  even  clear  that  such  a  flow  exists. 

33.  Static  Max-Flow  Min-Cut  Theorem 

The  next  section  will  make  use  of  the  well-known  max-flow  min-cut  theorem  for  static 
networks,  and  so  we  now  briefly  review  that  theorem  in  the  required  form.  By  a  static  net¬ 
work  we  mean  a  quadruple  N  =  (IV .  c ,  S  ,D ),  where  N  is  a  finite  set  of  nodes, 
c  =  (0^ j i  ,j  €A/ ).  Cjj  being  a  nonnegative  number  giving  the  capacity  of  link  (i  ,j ).  and  S ,  D 
are  disjoint  subsets  of  N  called  the  source  and  sink,  respectively.  By  a  link  we  mean  any  ele¬ 
ment  of  N2.  Some  authors  prefer  to  define  the  set  L  of  links  to  be  an  arbitrary  subset  on  N2 
and  to  define  c  on  L .  Their  model  can  be  converted  to  our  model  by  letting  =0  for  links 
( i  ,j  )  not  in  L  . 

A  feasible  flow  for  N  is  defined  to  be  an  assignment  /  *(/,;:  i  ,j  €N)  such  that 
—f  ji  —  f  ij  ^  c^j  and  f  x  j  =  0  for  all  j  €N—  (S  U D).  Note  that  if  /  is  feasible,  then  /„  =0 
and  — Cjj  <  f  ^  < ctJ .  Some  authors  omit  the  condition  f  tJ  =  —f  yl  and  instead  require  that  f  tj 
be  nonnegative.  Either  model  is  easily  converted  to  the  other.  The  value  of  a  flow  /  is  defined 
to  be  f  5  D  .  where  D  =  N  —D .  Since  /  jv  ,  =0  for  all  j  not  in  S  or  D .  it  is  easy  to  show  that 
/ 5  .d  ~  f  s3  ■  A  separating  cut  for  N  is  defined  to  be  a  cut  K  such  that  D  CK  and  S  C.K  . 


The  capacity  of  a  cut  K  is  defined  to  be  c%  jc  ■ 

The  max-flow  min-cut  theorem  states  that 

max  i  f  5  j>'-  f  is  a  feasible  flow}  =  min  {cjp  ^  :  K  is  a  separating  cut}. 

A  feasible  flow  which  achieves  the  above  maximum  is  called  a  maximum  flow,  and  a  separating 
cut  which  achieves  the  above  minimum  is  called  a  minimum  cut.  A  corollary  to  the  max-flow 
min-cut  theorem  is  that  if  /  is  a  maximum  flow  and  A"  is  a  minimum  cut,  then 

/  ,j  =  Cjj  for  all  i  €K  and  j  €K . 

Many  efficient  algorithms  exist  for  computing  a  maximum  flow  and  an  associated 
minimum  cut.  An  algorithm  which  is  fastest  for  dense  networks  is  the  one  in  [Malhotra  et  al.], 
which  requires  at  most  O  (  I  N  I  3)  computations.  Faster  algorithms  exist  when  the  network  is 
sparse.  For  an  excellent  introduction  and  references,  see  [Papadimitriou  and  Steiglitz,  1982]. 

3.4.  Dynamic  Max-Flow  Min-Cut  Theorem 

In  this  section  we  present  a  dynamic  version  of  the  max-flow  min-cut  theorem  and  show 
that  a  r-maximum  flow  can  be  computed  by  solving  a  static  maximum  flow  problem.  We  also 
derive  necessary  and  sufficient  conditions  for  a  flow  in  U  to  be  a  r-maximum  flow. 

We  define  a  dynamic  separating  cut  for  the  network  D  to  be  a  set-valued  function 
G :  (0,  T]-*2N  which  is  constant  on  each  interval  (m—  l.m],  m=l,...,7’.  and  satisfies 
d  €G  (r )  and  s  €G  (t )  for  all  t .  The  set  of  all  dynamic  separating  cuts  will  be  denoted  G. 
Note  that  G  is  a  finite  set.  For  each  r€(0,  T]  we  define  wT  on  G  by 
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w  r(G  )  -  f  bc  (/  )G  )(r ) 

o  — 

|rl— 1 

+  £  X  djim). 

m  j  €  G  (m  +  1)0(7  (m  ) 

where  fr]  denotes  the  smallest  integer  not  less  than  r.  We  define  a  t— minimum  cut  to  be  a  cut 
G  €G  such  that  wT(G  )  ^  wT(G  ')  for  all  G  '€G. 

Theorem  3.1 

For  each  r€(0.  7], 

max  { vT(u  ):  u  }  =  min  {wr(G  ):  G  6G}. 

Furthermore,  if  G  is  a  r-minimum  cut.  then  a  flow  u  €U  is  a  r-maximum  flow  if  and  only  if 

the  following  conditions  hold:  "" 


Ujj(t)  =  bij(t).  (i  ,j  )€G  (t  )xG  (t  ),  a.e.  t  6(0.r],  (3.4) 

Xj  {m  )  =  a;  (m  ).  j  €G(m  +l)HG  (m  ),  m  =  1,  .  .  .  ,  [t| —  1 .  (3.5) 

Xj  (m  )  =  0,  j  €G  (m  )DG  (m  +1).  m  =  1,  .  .  .  ,  [t|— 1.  (3.6) 

xj(r)  =  0.  j  €G  (r).  (3.7) 


Moreover,  a  r-maximum  flow  u  such  that  x  (r)  =  0,  and  a  r-minimum  cut.  can  be  computed  in 
O  (  I  N  I  3  I  T  I  3)  computations. 


The  proof  of  Theorem  3.1  will  be  based  on  the  static  max-flow  min-cut  theorem.  We  first 


present  three  lemmas. 
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Lemma  3.1 

Let  u  €U  .  G  €G.  and  r€( 0,  T ].  Then 

T 

vT(u)  =  f  uSU)0U)(t )  dt  -  £  Xj(t) 

0  y€G(r) 

|rl— 1 

+  £  £  Xj(m)  —  £  Xj{m)  .  (3.8) 

m—  1  j  €  G  (m  +l)f)G  (m  )  j  €  G  (m  )  PIG  (/n  +1 ) 

It  follows  that  vT(u)  ^  wr(G  ).  and  that  vT(u  )  =  wr(G  )  if  and  only  if  conditions  (3.4)  -  (3.7) 
of  Theorem  3.1  hold. 

Proof 

If  G{t)  —  [d  |  for  all  t  €(0.  T ].  then  (3.8)  holds  trivially.  Now  let  G  be  a  cut  in  G  not 

identically  equal  to  jV—  s  .  let  j  be  a  node  in  G  (m  )—s  for  some  m  £{  1 . T).  and  define  the 

cut  H  in  G  by  H  (m  )  =  G  (m  )  +  j  and  H(k)  =G(k)  for  all  k  ^ m  .  Then  letting  Q  (G  ) 
denote  the  right  side  of  (3.8)  as  a  function  of  the  cut  G  .  it  is  easily  verified  that 

m  -r 

Q(H)  —  Q(G)=  f  ujy  y  (r  )  df  +  ((m —l)«r)  —  X;  (/n-r). 

(m  -1)-t 

which  is  zero  by  the  definition  of  Xj .  Therefore,  since  any  cut  in  G  can  be  obtained  by  starting 
with  the  cut  identically  equal  to  \d  }.  and  repeating  the  above  procedure.  (3.8)  holds  for  all  cuts 
G  €G.  The  last  statement  of  the  lemma  now  follows  from  the  definition  of  wT(G  )  and  the  con¬ 
straints  (3.2)  and  (3.3).  ■ 

For  each  r€(0,  T],  we  define  Nr  to  be  the  static  network  (N T,  c  .  S .  D )  with  node  set 
N T  =  N  x[l . |t|}.  source  S  =  s  x{  1 . |t]},  sink  D  =  d  x[l . |r|},  and  with  the  capa¬ 

city  c  defined  as  follows,  where  m  «■  min  [m  ,  r}: 
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C(ifn).(j^n)  =  (m~r  ~m  +  1  )6,y(m).  m  =  1, .  .  .  ,  [r|. 

C(//»).0/*+l>  =  «y  )■  m=1 . M-1- 

all  other  links  having  zero  capacity. 

UK  is  a  separating  cut  for  NT,  we  define  Km  =  { j  €7V :  (j  ,m  )€  K }.  Note  that,  since 
D  CK  and  S  CK  .  we  have  d  £Km  and  s  £Km  for  all  m. 

Lemma  3.2 

If  K  is  a  separating  cut  for  NT  and  G  is  any  cut  in  G  such  that  G  (m  )  =  Km  for 
m  =  1,  ,  [t],  then  w  r(G  )  =  cy  x  • 

Proof 

This  is  easily  verified  from  the  definitions  of  w.(G  )  and  .  ■ 

Lemma  33 

Let  /  be  a  feasible  flow  for  Nr.  Then  the  flow  u  defined  by 
“ijit  )  -  t  (i/n)Xj  ~m  +  1).  t  €(m-l.m-T]. 

Uij(t)  =  0.  t  €(rj], 

is  in  (/ ,  and  the  storage  function  x  due  to  u  satisfies  Xj  (r)  =  0  for  all  y€iV.  Moreover. 
vr(«  )  =  f  d  ,d  ■ 

Proof 

Since  / =—fj,  and  / iy  ,  the  flow  u  clearly  satisfies  (3.1)  and  (3.2).  Since  /.vj  =0 

for  all  j  €/Vt— (S  UZ?  ).  we  have  for  m  =1,  ...  ,  |r|, 

f  (j  /n  ).(j  m  +1)  /  (j  m  -l).(j  jm  )  ~  f  (N  jn  ),(j  /n  ) 

=  uN  j  (m  «r)(m  —  m  +  l). 


where  we  take  f  a  m  ).(j  m  +i)=0  if  m  =0  or  |r|.  Summing  both  sides  of  the  above  equation  over 
m  .  we  get 


m 

/  (j  /n  ),(j  +1)  —  J“  ,)  j  ( m  -T). 

0 

Letting  m  =|r|  in  the  above  equation  gives  0  —  x}  (r).  Condition  (3.3)  for  r  =  1,  . . .  ,  (r|— 1  now 
follows  from 

0  —  — C  o  /n  +D.0  jn  )  ^  f  (j  /n  ).(]  /n  +1)  ^  c  0  /n  ).<y  ,m  +1)  ~  aj  (m  )■ 

Condition  (3.3)  for  all  t  €(0.r]  now  follows  from  the  facts  that  xy  is  affine  on  each  interval 
[m— lm-r]  and  that  a  j  is  lower-semicontinuous.  Now  since  Xj  (t)=0  and  uNj  (t  )=0  for  all 
f  >r,  we  have  Xj(t  )=0  for  all  t  >r.  and  so  condition  (3.3)  is  satisfied  for  all  t  €(0,  T  ].  There¬ 
fore  u  €U .  Finally, 

vr(u)  =  f  uN  4  (t  )  dt  =  ^  uN  —  m  +1) 

0  m  =1 

=  /  (Ns>).(dsn)  =  f  5  £>■ 

m  =1 


Proof  of  Theorem  3.1 

The  static  max-flow  min-cut  theorem  (Section  3.3)  implies  there  exists  a  feasible  flow  / 
and  a  separating  cut  K  for  NT  such  that  f  &  x>  ~  CK  jc  ■  Therefore,  Lemmas  3.2  and  3.3  imply 
there  exists  a  cut  G  €G  and  a  flow  u  €f/  such  that  vT(u  )  =  f  b  j>  =  cjf  #  —  wr(G  ).  Therefore, 
by  Lemma  3.1.  u  is  a  r-maximum  flow  and  G  is  a  r-minimum  cut.  The  theorem  is  now 
immediate  from  Lemma  3.1.  Since  the  network  Nr  has  no  more  than  \N  \  T  nodes,  the  flow  / 
and  the  cut  K  can  be  computed  in  0(IAM3r3)  computations  by  the  algorithm  in 


[Malhotra  et  al..  1978].  ■ 


3.5.  Computing  a  Flow  which  is  r -Maximum  for  all  r 

In  this  section  we  are  concerned  with  the  problem  of  finding  a  flow  in  U  which  is  a  T- 
maximum  flow  for  all  r€(0,  T  ].  At  this  point  it  is  not  even  clear  that  such  a  flow  exists.  We 
begin  by  singling  out  a  special  T-minimum  cut  for  each  time  t.  If  G  ,H  6G,  we  define  G  U H . 
G  C\H  €G  by  (G  UH)(t  )  =  G  (r  )UH  (t  )  and  (G  DH  )(t)  =  G  (t  )Off(c  ).  Similarly,  we  write 
G  C.H  if  G  (t  )<ZH  (t )  for  all  t . 

Lemma  3.4 

For  each  r,  wr  is  submodular,  i.e.,  if  G  .H  €G  then 

wr(G  GH)  +  wr(G  C\H)  <  wrCG)  +  wrCff  ).  (3.9)' 

It  follows  that  for  each  r  there  exists  a  r-minimum  cut  Rr  such  that  RTC.G  for  all  T-minimum 
cuts  G . 

Proof 

It  is  easily  verified  (and  well  known)  that  ifN  =  UV,c,S.Z))isa  static  network  as  in 
Section  3.3.  then  c*  jf  is  a  submodular  function  of  K .  It  therefore  follows  from  Lemma  3.2 
that  wr  is  submodular.  Now  suppose  G  and  H  are  both  T-minimum  cuts,  and  let 
a  =  wr(G )  =  wt(H  ).  Then  wT(G  and  wT(GC)/f)^a.  and  so  (3.9)  implies 

wr(G  U H )  —  wt(G  C\H)  —  a.  implying  that  G  U H  and  G  C\H  are  both  T-minimum  cuts. 
Therefore,  since  the  set  of  all  separating  cuts  is  finite.  Rr  is  given  by  the  intersection  of  all  t- 
minimum  cuts.  ■ 

The  cut  R  f  of  Lemma  3.4  will  be  called  the  canonical  T-minimum  cut.  We  will  show  that 
if  is  known  for  all  t,  then  a  flow  which  is  T-maximum  for  all  t€(0,  T]  can  be  obtained  by 


piecing  together  flows  which  are  each  r-maximum  for  a  single  r.  We  first  need  some  properties 
of  Rt. 

Lemma  3.5 

(a)  If  cr<r.  then  RaCRr.  That  is.  RT  is  increasing  in  r. 

(b)  Rr  is  a  left-continuous  function  of  r. 

(c)  Let  m  €  { 1 . f }  and  let  m -1  <<r<r^m  .  If  )  =  RT(m  ).  then  Ra  =  RT. 

(d)  Define  W  =  {r€(0.r]:  H/iJUll,  . .  .  ,T\.  Then  I  W  I  <  \N\T. 

(e)  Let  m€{l,...,r}.  Then  Rm+  is  an  m  -minimum  cut  and  a;(m)= 0  for  all 

j  €Rm+(m  +  l)n^m+(m ). 

(f )  For  each  r  €  (0.  T ].  R  T+  is  a  r-minimum  cut. 

(g)  Let  K  —  Rt(t)  for  some  r.  Then  for  all  /  such  that  d  €/  C.K  .  bg  #  (r)  <  bj  j(t). 

(h)  Let  m  € { 1 ,  .  .  .  ,T\  and  let  m  —1  <or<0  ,  and  suppose  G  is  a  cut  which  is  both  a- 

minimum  and  0-minimum.  Then  G  is  r-minimum  for  all  t  €[a,0].  Furthermore,  if 
G  =  Rt  for  some  r€[o,0],  then  G  =  R,  for  all  r  6(a,0]. 

(i)  Suppose  the  capacities  b i}  and  aj  are  rational  valued  with  common  denominator  p.  Let 

m  €|0 . r-li.  Z  =  £  \bij  (m  +1):  i  ,j  SN  },  and  a  =  m  +(pZ  )-1.  Then  for  all 

r  €(m  .a],  R,  =  i.e..  W  fl(m  .a)  =  0. 

Proof 

(a)  It  is  easily  verified  from  the  definition  of  w  T  that 


so 


r)  -  w<T(^(TU^T)  <  wt(Rt)  -  w(T(^T). 

Since  £T  is  a  r-minimum  cut  we  have  wT(RT)  =$  w  T{.R  R  T) ,  and  so  the  above  inequality 

implies  w  ^(RT)  ^  Furthermore,  since  w  a  is  submodular,  we  have 

w(7(.R(rU.Rr)  +  •w0.(j?(Tni?T)  <  w^R^)  +  waQ?T). 

and  so  w<r(.£{rn.£T)  ^  w  ^(R^).  Therefore,  since  R„  is  the  canonical  tr-minimum  cut. 
RaCR0nRr,  i.e.,  R„CRt. 

(b)  Since  the  set  G  of  all  separating  cuts  is  finite,  part  (a)  implies  that  Rr  can  change  only 
finitely  many  times  as  r  increases  from  0  to  T .  Therefore,  given  t€( 0.  T  ].  there  exists  a  posi¬ 
tive  number  €  such  that  R,  —  G  for  some  G  €G  and  all  t  €(r— €,r).  Then  G  is  a  r-minimum 
cut.  hence  w,  (G  )  ^  w(  (RT)  for  all  t  €(r— €.r).  Therefore,  since  w,  (G  )  is  clearly  left- 
continuous  in  t.  wT(G)  ^  wTCRr).  and  so  by  the  definition  of  Rr,  RrCG  —  R,  for  all 
t  €(t— €.t).  Part  (a)  now  implies  RT  =  R,  for  all  t  €(r— e.r),  and  so  Rr  is  left-continuous  in  t. 

(c)  Since  R^m  )  =  RT(m  ).  the  definition  of  w,  implies 

—  w^iRa)  =  w  t(R  r)  —  w  a(if  T). 

Now  by  the  definition  of  R&  we  have  w  ,/.£<,.)  <  RT),  and  so  the  above  equation  implies 
wr(Rv)  <  wt(Rt).  But  by  part  (a).  R„<ZRT.  Therefore  by  the  definition  of  R T,  Ra  =  Rr. 

(d)  If  r€(m  — l.m  )  for  some  m€{l . T\,  and  RT  ^  RT+,  then  parts  (a)  and  (c)  imply 

RT(m  )CZRT+(m  )  and  RT(m  )  ^  /fr+(m  )  respectively,  and  so  Rr+(m  )  has  at  least  one  more  node 
than  RT(m ).  Therefore,  since  for  each  r.  Rr(m )  must  contain  d  but  not  s ,  there  can  be  no 
more  than  I I  — 2  points  r€(m— l,m)  such  that  RT^RT+. 

I  W  I  <  (IN  1-2 )T  +T  ^  IN  IT. 


Therefore, 
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(e)  Let  G  be  an  m  -minimum  cut  such  that  G  (m  +1)  =  G  (m  ).  (Since  wm  (G  )  does  not 
depend  on  G  (m  +1),  such  a  cut  exists.)  Then 

wm  +(G  )  =  w,„(G  )  ^  wm(jRm+)  ^  wm  +(Rm+). 

But  since  Rm  +  is  an  (m  +€)-minimum  cut  for  sufficiently  small  e>0,  wm+(Rm+)  ^  wm+(G). 
Therefore.  wm(G  )  =  wm  (Rm  +).  implying  Rm  +  is  an  m  -minimum  cut.  and 
wm  (Rm  +)  =  wm  +{Rm  +).  which  from  the  definition  of  w  T  implies  a;  (m  )  =  0  for  all 
j  £Rm  +(m  +1  )D/?m+(m  ). 

(f)  If  t  is  an  integer,  then  (f)  follows  from  (e).  Now  suppose  r€(m—lsn)  for  some 
integer  m.  Then  (f)  follows  from  the  fact  that  w,  (G  )  is  continuous  in  t  on  (m—  l.m).  for 
wT(i?T)  ^  w  t(Rt+)  and  w  r+f(R  T)  ^  w  r+f(R  r+)  for  sufficiently  small  €>0,  and  so 
wr(£r)  =  wT(^T+). 

(g)  Let  r€(m  — l,m  ]  for  some  integer  m  ,  let  K  =  Rt{t),  and  let  J  be  such  that  d  6/  C.K . 
Define  the  cut  G  by  G(m)  =  J  and  G  (k  )  —  R  T(k  )  for  all  k  ?*m  .  Now  if  bj  j  (t)  <  2>jp  x  (t), 
then  it  is  clear  from  the  definition  of  wT  that  wr(G  )  <  wr(Rr),  contradicting  the  definition  of 

Rr- 

(h)  Let  H  be  any  cut  in  G.  From  the  definition  of  wT  it  is  clear  that  w,  (G  )  and  w,  ( H  ) 

are  affine  in  t  for  t  €[a,0].  Therefore,  since  w  a(G  )  ^  w„(/f )  and  w^G  )  ^  ),  we  have 

w,  (G  )  ^  w,  ( H )  for  all  t  €[a.j3].  Therefore  G  is  a  t  -minimum  cut  for  all  t  €[a,0].  Now  sup¬ 
pose  G  -  Rt  for  some  T€[a.j8]  and  let  H  be  a  cut  in  G  such  that  H  CG  and  H  ^ G  .  Then 
wr(G  )  <  wT(tf)  and  w,  (G  )  ^  w,  ( H  )  for  all  t  €[a./3).  which  implies  w,  (G  )  <  w,  ( H  )  for  all 
t  €(a./3).  Therefore  G  =  Rr  for  all  t  €(a./3).  By  part  (b)  we  also  have  G  =  R$. 

(i)  By  part  (a),  it  suffices  to  show  Ra  =  Rm+.  By  the  definitions  of  Ra  and  Rm+. 
w0( Ra)  <  wa(£m+)  and  wm+(^0)  ^  wm+(Rm+).  We  distinguish  two  cases.  First  suppose 


wm+(Ra)  =  +(.£„,+).  Then  the  affineness  of  w,  on  (m  ,m  +1]  implies  w,  (^a)  <  w,  (7?m+)  for 

all  t  €(m .«].  Therefore,  Ra  is  an  (m  +€)-minimum  cut  for  sufficiently  small  e>0,  an  so  part 
(h)  implies  Ra  =  R,  =  Rm+  for  all  t  €(m  .a]. 

Now  suppose  wm+(Ra)  >  wm+(Rm+).  Since  the  capacities  btJ  ( t  )  and  ay  (t  )  are  assumed  to 
be  rational  with  common  denominator  p  .  it  is  clear  that  w„,  +(G  )  also  has  this  property  for  any 
cut  G  .  Therefore,  wm+(Ra )  ^  wm+(Rm+)  +  p~l.  Now  from  the  definition  of  wT  we  have 

yva(Rm+)-wm+(Rm+)  =  (a  -m)bxjc(m+l)  =  ( pZ  )_16jf  jr  (m +l)  <  p~x. 

where  AT  =  +(/n  + 1 ).  Therefore. 

wa<X,+)  <  wm+(RmJ  +  p~l  <  yvm+iRa)  ^  w 

and  so  Rm  +  is  an  a-minimum  cut.  Part  (h)  therefore  implies  Rm+  =  R,  -  Ra  for  all  t  €(m  ,a]M 

We  can  now  prove  a  simple  necessary  and  sufficient  condition  for  a  flow  to  be  r-maximum 
for  all  t. 

Lemma  3.6 

A  flow  u  €f/  is  r-maximum  for  all  r€(0,  T]  if  and  only  if  u  is  T  -maximum  and  Xj  (t  )=0 
for  all  j  €R,  ( t )  and  all  t  €(0,  T ]. 

Proof 

That  the  conditions  are  necessary  is  immediate  from  condition  (3.7)  of  Theorem  3.1.  Now 
let  (of./3]  be  an  interval  such  that  R,  is  the  same  cut  G  for  all  t  €(ar,0].  By  virtue  of  Lemma 
3.5(d)  it  suffices  to  show  that  if  u  is  a  0-maximum  flow  such  that  Xj  (r  )=0  for  all  j  €G  (t  )  and 
all  t  €(a,0],  then  u  is  a  t  -maximum  flow  for  all  t  €[<x.0],  for  then  we  may  start  with  0=7*  and 
propagate  backward  in  time.  So  let  u  be  such  a  flow  and  let  t  €(a,0].  Since  u  is  0-maximum, 
conditions  (3.4)  -  (3.6)  of  Theorem  3.1  hold  with  r=0.  and  hence  with  T=t .  Now  by 
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assumption,  condition  (3.7)  also  holds  with  r=t .  Therefore,  since  G  is  at  -minimum  cut.  u  is 
a  t  -maximum  flow  by  Theorem  3.1.  Therefore,  u  is  t  -maximum  for  all  t  €(a.0],  and  since 
v,  ( u  )  is  continuous  in  t .  u  is  t  -maximum  for  all  t  ■ 

Remark 

Note  that  Lemma  3.6  does  not  imply  there  exists  a  flow  which  is  T-maximum  for  all  r;  this 
will  be  implied  by  the  next  theorem.  But  Lemma  3.6  does  suggest  a  method  for  computing 
such  a  flow  once  R,  (r  )  is  known  for  all  t :  set  bj  ( t  )=0  for  all  j  £Rt  (.t )  and  all  t .  and  compute 
a  T -maximum  flow  using  the  method  in  Section  3.4.  However,  since  the  corresponding  static 
network  may  have  as  many  as  I  N  I  2T  nodes,  this  is  not  the  best  approach. 

For  each  j  €N  we  define  the  function  qj  on  (0,  T  ]  by 
qjU)  =  max  {r:  r^r  and  j  €R  T(r  )). 

Since  RT(t  )  is  increasing  in  r.  the  above  set  is  empty  if  and  only  if  j  €R,  (t  ),  in  which  case  we 
set  qj  (c  )=0.  Note  that  if  0<t<T  .  then  qj(t  )=r  if  and  only  if  j  €^T+(t  )  —  Rr(.t  ).  and  that 
qj(,t)—T  if  and  only  if  j£Rr(t).  Thus,  if  q3{t  )=r  >0.  then  r  is  in  the  set  W  defined  in 
Lemma  3.5(d). 

For  each  r€W .  let  uT  be  a  T-maximum  flow  such  that  xT(r)=0,  where  xT  is  the  storage 
function  due  to  uT.  The  next  theorem  constructs  a  flow  which  is  T-maximum  for  all  t  by  piec¬ 
ing  together  the  flows  uT.  t€W\  and  the  following  flow.  Let  u°  be  a  (not  necessarily  feasible) 
flow  which  changes  only  at  points  in  W  and  is  such  that  for  each  t  €(0,  T]  the  static  flow 
defined  by  /(J  =  (f  )  is  a  maximum  flow  for  the  static  network  (N  ,  b  (r  ),  R,  (t  ).  d  ).  Since 
W  has  no  more  than  \N  \  T  points,  the  flow  u°  can  be  computed  by  solving  no  more  than 
IN  IT  static  maximum  flow  problems,  each  on  a  network  with  IN)  nodes.  Thus  (using  the 
algorithm  in  [Malbotra  et  al.]),  u°  can  be  computed  in  O  ( I N  I  *T)  computations. 


Theorem  3.2 


The  flow  u ' ,  defined  by 


uljU) 


Uj(t  ).  qj  (t  )  =  T. 

(*  ).  9,  (t  )  >  9;  (*  ). 

).  9i(0  <  qj  {t  ). 


is  a  r-maximum  flow  for  all  r€(0,  T  ]. 

Proof 

Note  that  if  ^(r  )=r  and  0 <t<T  .  then  Rr^RrJr.  but  by  Lemma  3.5(f).  RT  and  RT+  are 
both  r-minimum  cuts.  This  fact  will  be  used  repeatedly  in  the  proof. 

We  first  show  that  u'  is  in  U .  Clearly  conditions  (3.1)  and  (3.2)  are  satisfied.  To  verify 
condition  (3.3).  we  first  show  that  for  almost  all  t .  u’j  (t  )=uj(t  ).  where  r=qj  ( t  ).  This  is  true 
by  definition  if  q,  it  )=qj  (t ).  We  distinguish  three  more  cases. 

First  note  that  if  (r  U  )=T  then  j  £Rr(t  )  and  i  £RT(t  ).  Therefore,  by  Theorem 
3.1.  UjTi(t  )—bji  (r  ).  i.e.  u,J(r  )=-&;;  (r  ),  for  almost  all  t  such  that  0 (r  )<?;  (f  )=t. 

Next  note  that  if  qt  (f  )>?,  (t  )=r>0,  then  i  €Rr+(t  )  and  j  6 Rr+(t ).  Therefore,  since  Rr+ 
is  a  r-minimum  cut.  Theorem  3.1  implies  that  Ujj(t  )=b^  (f  )=«,’;(£  )  for  almost  all  t  such  that 
qt  (t  )>qj  U  )=r>0. 

Finally,  suppose  that  if  ?,  (f  )>?;  (t  )=0.  Then  i  €R,  (e  )  and  j  €/?,  (f ).  Therefore,  since 
Lemma  3.5(g)  implies  that  R,  (t )  is  a  minimum  cut  for  the  static  network  (N  ,b(t ).  R,  (r  ).  d  ). 
the  corollary  to  the  static  max-flow  min-cut  theorem  discussed  in  Section  3.3  implies  that 
Uijit  )-bu  (r  )-Ujj  (r ).  This  completes  the  proof  that  for  almost  all  t .  u‘j(t  )=uj(t  ).  where 

T=qj  ( t  ). 


Let  x  *  and  x T  denote  the  storage  functions  due  to  u  *  and  u  T,  respectively.  We  next  show 
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that  if  qj  (t  )=r  then  x’U  )=x/(t ),  where  we  take  xfic  )=0  by  convention.  This  will  establish 
|  condition  (3.3).  To  show  this,  we  first  note  that  for  almost  all  t.  u^,  j(t  )-urf j(t ),  where 

T-qs  ( t ).  Therefore,  since  xT(0)=0  for  all  r.  it  suffices  to  show  that  if  <r=qi  (r  )^q2  (t  +)~r. 
then  x/(t  )=x/(f  ).  We  distinguish  three  cases. 


« 


Let  <r=q1  it )  and  r=^;  (r  +).  We  first  suppose  cr>r>0.  Then  the  definitions  of  r  and  cr 
imply  j^R^it)  and  y€i?(J(r+)  respectively.  Thus  t  equals  some  integer  m.  and 
j  €Ra(m  +1)C\R  ).  Therefore,  by  Theorem  3.1,  xfit  )=aj  (t ).  Now  j  €R  ff(t  )CR  T+(t ). 

and  the  definition  of  r  implies  i  €Rr+(t  +).  Thus  j  €Rr+(m  +l)H^  r+(m  ).  Therefore,  since  Rr+ 
is  a  r-minimum  cut.  Theorem  3.1  implies  xfU  )=<ij  ( t  )=x  “(t  ). 

Next  suppose  0^cr<r.  Then  the  definitions  of  <r  and  r  imply  j  €Rr(t )  and  j  €Rr(t  +), 
respectively.  Thus  t  equals  some  integer  m  .  and  j  €Rr(m  )C\R  T(m  +1).  Therefore,  Theorem 
3.1  implies  x/(£  )=0.  Now  if  cr>0.  then  j  €R  r(t  +)Cif  „+(.t  +),  and  the  definition  of  cr  implies 
jZRv+it).  Thus  j  e^^+Cm  )fl£  <,+(.m  +l).  Therefore,  since  R a+  is  a  cr-minimum  cut. 
Theorem  3.1  implies  xfit  )=0 =X;r(£  ).  If  cr=0,  then  xf^t  )=0  by  convention. 

Finally,  suppose  <r>r=0.  Then  x/(t  )=0  and  we  need  to  show  x°(t  )=0.  By  definition  of 
cr.  (T^t .  If  cr—t ,  then  xfit  )=Xy<r(cr)=0  by  our  choice  of  uv.  Now  suppose  <r>t .  Then  the 
definitions  of  < r  and  r  imply  ;  and  y€i?r+(r+)  respectively.  Thus  t  equals  some 

integer  m.  and  j  €Rm+(m  +l)DRm+(m  ).  Therefore,  by  Lemma  3.5(e),  Oj{.t)~  0,  and  so 
x,  (t  )=0. 

We  have  shown  that  u'  is  feasible  and  that  x’(t  )=0  when  qj  ( t  )=0.  i.e.,  when  j  €R,  (t ). 
Therefore,  by  Lemma  3.6,  to  show  u '  is  a  r-maximum  flow  for  all  r€(0.  T  ],  it  suffices  to  show 
that  u'  is  a  7 -maximum  flow.  We  will  verify  conditions  (3.4)  -  (3.7)  of  Theorem  3.1  with 
r-T  and  G  -Rr  ■  Condition  (3.4)  holds  by  the  definition  of  u‘ .  for  if  i  SRr(t )  and  j  €Rr(t  ). 
then  qj  (t  )<qt  (t  )~T .  Conditions  (3.5)  and  (3.6)  hold  because  x](t  )=xf(t  )  when  ?;(r)=7. 
i.e.,  when  j  £Rr{t  ).  Finally,  since  x’(t  )=0  for  j  €/?,(£  ),  condition  (3.7)  holds.  ■ 


Corollary  3.1 

If  W  is  given,  then  a  flow  which  is  r-maximum  for  all  r€(0.  T ]  can  be  computed  in 
O  ( I N  1 4r4)  computations. 

Proof 

For  each  r€  W ,  Theorem  3.1  implies  that  uT  and  RT  can  be  computed  in  O  ( I N  I  37*3)  com¬ 
putations.  Therefore,  since  IWI^IAMT.  ur  and  RT  can  be  computed  for  all  r€W  in 
O  ( I  N  1 4r4)  computations.  Also  recall  that  u°  can  be  computed  in  O  ( I N  I  *T )  computations. 
Therefore,  since  the  functions  q ,  j  6  N  ,  are  determined  by  (R  T:  r  €  W  )  the  corollary  follows 
from  Theorem  3.2.  ■ 

It  remains  to  find  a  method  for  computing  W .  The  proof  of  the  following  theorem  is  sum¬ 
marized  by  the  algorithm  FINDW  of  Figure  3.1  (page  58).  which  computes  W  in  0(  I N  I  4T4) 
computations.  FINDW  makes  T  calls  to  the  function  PARTWCa.jS.^,,,.^^).  which  computes 
W  n[a./3)  in  Oi  I  N  l4r3)  computations.  FINDW  together  with  Theorem  3.2  provide  a  method 
for  computing  in  O  (  I N  1 4T4)  computations  a  flow  which  is  r-maximum  for  all  r€(0,  T  ]. 

Theorem  3.3 

Suppose  the  capacities  6,  j  and  a ,  are  rational  valued.  For  each  m€(l . T ). 

W  fl(m— lyn  ]  can  be  computed  in  0(  1 2V  l4T3)  computations.  It  follows  that  W  can  be  com¬ 
puted  in  O  (  I N  I  4r4)  computations. 

Proof 

Let  m€{l,  .  .  .  ,T\,  letm— l<a</3^m,  and  let  z  =  I  R  p(m  )  —  Ra(m  )  I .  We  will  assume 
that  Ra  and  R $  are  given  and  show  that  W  n[a,/3)  can  be  found  by  computing  R,  for  at  most 
max{  2—1,0  }  values  of  t .  First  suppose  z  =0.  Then  Ra(m  )=Rp(m  ),  and  so  by  Lemma  3.5(c), 
Ra~R Therefore,  by  Lemma  3.5(a).  R,  -Ra  for  all  t  €[a,0],  and  so  we  know  W  n[ar,0)=0. 

Now  suppose  r^l.  We  will  show  by  induction  on  z  that  W  n[a.j3)  can  be  found  by  com- 
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puting  R,  for  at  most  z—l  values  of  t .  Since  z  ^1,  Ra(m)  is  contained  in  but  not  equal  to 
R  p(m  ).  Therefore,  by  the  definitions  of  Ra  and  R  p.  w  a(R  a)^w  a(R  p)  and  w p(R „)>w p(R p). 
Now  since  w,  ( R a )  and  w,  (R  p)  are  both  affine  in  t  on  [or,0].  there  is  a  unique  point  f  €[a.0)  such 
that  wT(/Jcl)=w'T(.ff  p).  We  distinguish  three  cases. 

First  suppose  z  =1.  Then  it  is  clear  from  Lemma  3.5(a.b)  that  R,  =Ra  for  all  t  G[ck.t]  and 
R,  =Rp  for  all  t  €(r./3].  Therefore  we  know  that  W  n[a./3)  =  {r}  without  computing  R,  for 
any  t  . 

Now  suppose  T=a.  Then  R  p  is  an  a-minimum  cut.  and  so  by  Lemma  3.5(h).  R,  —Rp  for 
all  t  €(ck./3],  implying  W  n[a./3)  *  {a}. 

Next  suppose  z  >  1  and  TJ*a.  and  compute  RT.  We  distinguish  two  cases.  First  suppose 
wT(Rj=wr(Ra)=wT(Rp).  Then  Ra  and  Rp  are  both  r-minimum  cuts,  and  so  by  Lemma 
3.5(h).  R,  =Ra  for  all  t  €[a.r]  and  R,  =R  p  for  all  t  6(r./3].  implying  W  ft [a.0)=T.  Next  sup¬ 
pose  wT(RT)<wT(Ra)=wT(Rp).  Then  Ra7±RTi*Rp.  and  so  by  Lemma  3.5(c), 
R  a(m)^R  T(m)^R  p(m).  Therefore.  defining  Zj  =  I  RT{m  )  —  Ra(m  )  I  and 
z  2  =  I  R  p(m  )  —  R  T(m  )  I .  z  j  and  z  2  are  both  less  than  z  .  Therefore,  by  the  inductive  assump¬ 
tion.  W  n[<*,T)  (respectively  W  n[r,/3))  can  be  found  by  computing  R,  for  at  most  z  t— 1 
(respectively  Z2— 1)  values  of  t .  Therefore,  since  Rr  also  had  to  be  computed.  W  n[a.0)  can  be 
found  by  computing  R,  for  at  most  (z  j— l)+(z2— 1)+1  =  2—1  values  of  t . 

We  can  now  compute  W  fl(m—  l,m]  as  follows.  Let  a  be  as  in  Lemma  3.5(0,  with  m 
replaced  by  m— 1.  so  that  we  know  W  fl(m— l,a)=0.  Now  compute  Ra  and  Rm  and  apply  the 
above  method  to  find  Wn[a,m)  by  computing  R,  for  at  most  z—l  values  of  f.  where 
z  =  I  Rm  (m  )  —  Ra(m  )  I .  Since  by  definition,  this  gives  us  WfKm— l,m].  Since 

z^\N  1—2  and  since  Theorem  3.1  implies  that  Rt  can  be  computed  for  a  single  t  in 
0(INI37’3)  computations,  this  method  computes  WCHm—  l,m]  in  0(IAH4r3)  computa¬ 


tions.! 


procedure  FINDW 
begin 

let  p  be  a  common  denominator  of  ( btJ  (m  ).  ay  (m  ):  i  ,j  €jV ,  1  <  T ) 

for  m  =1,  .  .  .  ,T  do 
begin 

Z  :=  £  {bij  ( m  ):  i  ,j  €N  }  : 
a  :=  m  —  1  +  ( pZ  )-1  ; 
compute  Ra.Rm  : 

Wm  :=  PARTWCa.  m  .  Ra.  Rm  )U{m  }  ; 

end 

W  :=  IVj  U  •  •  •  U  WT  : 

end 

function  PARTWCa,  /3.  Ra.  R  p) 

(comment:  returns  IV  n[a,0)  ) 
begin 

let  m  be  such  that  m—  1  <  a  <  0  ^  m  ; 
if  Ra  -  Rp  then  return  0  : 
else  begin 

compute  t€[ a.0)  such  that  wT(^0)  =  w r(R p)  ; 

if  I R p(m  )  —  Ra(m  )  I  =1  or  r  =  a  then  return  {r}  ; 

else  begin 

compute  R  T  ; 

if  RT  =  Ra  then  return  {r }  : 
else  begin 

Wa  :=  PARTWCa,  r,Ra,Rr)  : 

Wp  :=  PARTWCt,  0 .Rr.Rp): 

return  VV^UW^  ; 

end 

end 

end 

end 


Figure  3.1.  Algorithm  FINDW  for  computing  W. 
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Remarks 

(1)  When  the  algorithm  FINDW  is  executed,  Rr  and  a  r-maximum  flow  are  computed  for 
each  r€W.  These  can  be  saved  so  that  they  need  not  be  recomputed  when  Theorem  3.2  is  used 
to  construct  the  optimal  flow  u' . 

(2)  If  b,j  and  are  constant  functions  of  time,  then  by  scaling  time  we  may  assume 
7  =  1.  In  this  case,  it  is  easy  to  see  that  Rr  is  constant  in  r  and  that  any  constant  flow  which  is 
7  -maximum  is  also  r-maximum  for  all  r€(0,  7  ].  Therefore,  the  case  of  constant  capacities  is 
not  very  interesting.  However,  if  we  allow  nonzero  initial  storage,  we  show  in  the  next  section 
that  the  case  of  constant  capacities  is  equivalent  to  the  case  with  zero  initial  storage  and  7=2. 

(3)  We  show  that  the  property  of  being  a  r-maximum  flow  for  all  t  characterizes  the 
solutions  of  the  linear  cost  problem, 

maximize  { J  t(u  ):  u  €£/  }  (Pi) 

where 

r 

J1(u)s  J  (7  —t)  uN  d(t)  dt. 

o 

Integrating  by  parts,  we  obtain 

T  T  T 

J i(u  )  =  f  f  uy 4  (t  )  dt  dr  =  J  vr(u  )  dr. 

0  0  0 

Therefore,  since  Theorem  3.2  implies  there  exists  a  flow  which  is  r-maximum  for  all  r€(0,  T  ], 
and  since  vr(u  )  is  continuous  in  r,  a  flow  u  €t/  solves  problem  (Pi)  if  and  only  if  u  is  a  t- 
maximum  flow  for  all  r€(0.  T  ]. 


3.6.  A  Minimum-Delay  Problem 

Suppose  that  instead  of  letting  bsj  ( t  )  represent  the  maximum  rate  at  which  commodity  can 
enter  node  i  from  the  source,  we  let  bS)  (t  )  represent  the  rate  at  which  commodity  must  enter 
node  i  from  the  source.  (In  practice,  bsj  (.t  )  may  represent  the  rate  at  which  commodity  is  gen¬ 
erated  at  node  j  or  enters  node  j  from  outside  the  network).  Thus  we  let  U 0  denote  the  set  of 
flows  u  €U  such  that  for  each  j  €iV  .  the  following  conditions  hold: 

uSj  (r )  =  bsj  ( t ),  a.e.  t  €(0,  T  ], 

Xj(T)  =  0. 

It  follows  from  the  definition  of  Xj  (r)  that  if  u  €f/0  then 

r 

xj  ( r )  =  f  [uv_s  j  (t )  +  bsj  (r )]  dt. 

o 

Assumption 

We  assume  that  U 0  is  nonempty.  Furthermore,  since  usd  does  not  contribute  to  the  storage 
x  ,  we  assume  for  convenience  that  usd  —  bsd  =  bds  =  0. 

Consider  the  cost  functional 

r 

J  (u  )  =  f  x. v  (t )  dt 

0 

which  represents  the  total  delay  due  to  the  flow  u  .  We  are  interested  in  the  problem: 

minimize  {/  (u  ):  u  €{/0} .  (P) 

The  following  proposition  characterizes  the  solutions  to  problem  (P). 


Proposition  3.1 

For  each  r€(0.  T ],  a  flow  u  €U 0  minimizes  xiV(r)  on  U0  if  and  only  if  u  is  a  r-maximum 
flow.  It  follows  that  a  flow  u  is  a  solution  to  problem  (P)  if  and  only  if  u  is  a  r-maximum 
flow  for  all  t€(0,  T ]. 

Proof 

For  any  u  €t/  we  have,  using  usd  ( t  )=0. 

T 

xn(t)~  f  uSrN(.t)dt  —  vT(u). 
o 

Therefore. 

T  T 

vr(u)  =  J  us  jy(t )  dt  —  xN(T )  ^  J  bs  )  dt  , 
o  o 

where  equality  holds  if  and  only  if'u  €{/ 0.  Therefore,  since  U0  is  nonempty,  a  flow  u  6 U  is  in 
Ido  if  and  only  if  u  is  a  T  -maximum  flow.  In  particular,  by  Theorem  3.2,  there  exists  a  flow 
u  €f/ o  which  is  r-maximum  for  all  r€(0,  T ].  This  implies  that  for  each  r, 

max  {vT(u  ):  u.  6t/0}  -  max  {v.(u  ):  u  €(/ }. 

Therefore,  since  for  u  €(J 0  we  have 

r 

xn(t)  =  f  bs  N(t)dt  —  vT(u). 
o 

a  flow  u  €6/ o  minimizes  x^ix)  on  U0  if  and  only  if  u  is  a  r-maximum  flow.  The  last  statement 
of  the  proposition  now  follows  from  the  continuity  of  x  .  ■ 


Remarks 


(1)  In  view  of  Proposition  3.1,  a  solution  to  problem  (P)  can  be  computed  in  0(  t  N  1 4T4) 
computations  using  Theorem  3.2  and  the  algorithm  FINDW. 

(2)  Recall  that  we  are  requiring  that  the  initial  storage  be  zero.  No  generality  is  lost  by 
this  assumption  since  if  we  are  given  an  initial  storage,  we  may  let  b:j  ( t  )=0  for  i ,)  —s  and 
t  €(0.l],  and  choose  bst  ( t )  for  t  €(0,l]  such  that  it  generates  the  given  storage  at  time  1. 

(3)  Suppose  that  the  capacities  btJ  and  a ;  are  constant  functions  of  time,  but  that  a 
nonzero  initial  storage  is  given.  By  scaling  time  and  using  the  method  of  remark  (2).  we  see 
that  this  case  is  equivalent  to  an  instance  with  zero  initial  storage  and  T  =2.  Therefore,  for  the 
case  of  constant  capacities  and  nonzero  initial  storage,  a  solution  to  problem  (P)  can  be  obtained 
in  O  ( I  N  I  4)  computations.  This  is  exactly  the  case  considered  in  Chapter  2,  where  the  algo¬ 
rithm  OPTFLO  computes  a  minimum-delay  flow  also  in  O  ( I N  I  4)  computations. 
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4.  MAXIMUM  FLOWS  IN  CONTINUOUS  NETWORKS 

4.1.  Introduction 

In  this  chapter,  a  generalized  network  is  defined  for  which  the  set  ft  of  nodes  can  be  any 
set  (including  a  continuum),  and  a  flow  is  an  additive  set  function  defined  on  subsets  of  ft2.  A 
theory  of  flows  is  developed  which  generalizes  some  known  results  for  classical  networks. 

The  chapter  is  organized  as  follows.  Section  4.2  discusses  the  extension  of  an  additive  set 
function  dominated  by  a  submodular  function  (Theorem  4.1).  The  continuous  network  (called 
a  product-algebra  network)  is  defined  in  Section  4.3.  and  a  feasible  flow  theorem  (Theorem  4.2) 
which  generalizes  the  classical  feasible  flow  theorem  is  presented  in  Section  4.4.  Source-sink 
networks  are  considered  in  Section  4.5,  where  a  continuous  version  of  the  max-flow  min-cut 
theorem  (Theorem  4.3)  is  proved.  Also  in  Section  4.5.  the  existence  is  proved  of  a  flow  which 
simultaneously  maximizes  the  flow  into  each  set  in  a  given  increasing  family  of  sets  (Theorem 
4.4).  Two  methods  are  then  given  for  decomposing  the  problem  of  finding  such  a  flow  into  two 
subproblems.  Some  of  these  results  will  be  applied  in  Chapter  5  to  dynamic  networks  with 
continuously  varying  capacities. 

We  remark  that  a  version  of  the  max-flow  min-cut  theorem  similar  to  the  one  in  Section 
4.5  has  been  developed  independently  by  Neumann  (1984.  1985).  However,  the  theory  of  con¬ 
tinuous  networks  given  in  Chapter  6  is  more  general  than  that  of  Neumann. 

4.2.  Submodular  Set  Functions  and  Sublinear  Functionals 

Let  ft  be  any  set  and  let  A  be  an  algebra  of  subsets  of  ft,  i.e.,  A  is  closed  under  finite 
unions,  finite  intersections,  and  complementation.  A  set  function  fi :  A—K  is  said  to  be  addi¬ 
tive  if  fi(A  UZ? )  =  fj.(A  )  +  fi(B )  whenever  A  and  B  are  disjoint.  A  set  function 
v  :  A-*E  U{+oo}  is  said  to  be  submodular  if  for  all  A  ,B  6  A, 


u(A  US  )  +  v(A  C\B  )  ^  v(A  )  +  v(B  ). 

We  let  add(A)  denote  the  set  of  all  additive  set  functions  v  on  A  and  we  let  sub(A)  denote  the 
set  of  all  submodular  functions  v  on  A  such  that  v(0)= 0.  If  (i 6  add(A).  then 
/jl(A  U  B  )+fi(A  C\B )  =  fi(A  )+yu(5  ).  and  so  add(A)Csub(A). 

Let  simp( fi .  A)  denote  the  linear  space  of  all  simple  functions  on  (fl.A),  i.e..  the  set  of 
finite  linear  combinations  of  indicators  of  sets  in  A.  We  let  IA  denote  the  indicator  function  of 
the  set  A  .  If  /x€add(A)  and  /  €simp(fl.  A)  with 

/  =  T.°‘iIAl 

■  =i 

we  define  * 

M/  )  =  JfdH  =  (4.1) 

i=  l 

* 

Then  m  is  a  linear  functional  on  simp(  fl .  A).  Conversely,  if  l  is  a  linear  functional  on 
simpCn.  A),  then  l  =  /a,  where  fi  is  the  additive  set  function  defined  by  /i(A  )  =  lUA  ).  There¬ 
fore,  we  may  identify  add(A)  with  the  space  of  linear  functionals  on  simp(0.  A). 

A  functional  <t>  on  simp(O.A)  is  said  to  be  sublinear  if  the  following  conditions  are 
satisfied  for  all  simple  functions  /  ,  g  : 

<f>(af  )  =  ),  a  ^  0,  (nonnegative  homogeneity) 

<£(/  +g  )  <  4>{f  )  +  <f>(g  ).  (subadditivity) 

If  p€sub(A)  we  define  v  on  simpCH.  A)  by 


KK*. 


Kf  )  =  f  f  (o»)  ^  a)  —  v(fi)/(a<0)]  d  a. 


/  £  0. 


(4.2) 


where  /  denotes  the  indicator  function.  Note  that 

Kf  +  a)  =  £(/)  +  aKCl).  f  *S0.  a 6iR.  (4.3) 

Note  that  if  /z€ add(A),  then  f*  -  jjl.  Thus,  the  definition  of  ?  generalizes  the  concept  of  integra¬ 
tion. 

Lemma  4.1 

Let  v€sub(A).  Then  v  is  sublinear. 

Proof 

By  (4.3)  it  suffices  to  show  v  is  sublinear  on  the  nonnegative  simple  functions.  A  proof  of 
this  can  be  found  in  Choquet  (1953).  ■ 

We  now  present  some  properties  of  submodular  set  functions  which  will  be  needed  in  later 
sections.  The  following  theorem  regarding  the  extension  of  an  additive  set  function  dominated 
by  a  submodular  set  function  parallels  the  extension  form  of  the  Hahn-Banach  Theorem  [Roy- 
den]. 

Theorem  4.1 

Let  A.  B  be  algebras  of  subsets  of  Cl  such  that  BCA.  Let  v€sub(A).  and  suppose 
/x€add(A)  is  such  that  ^t^vonB.  Then  there  exists  n  €add(A)  such  that  /z  is  an  extension  of 
fi  and  £i  ^v. 

Proof 

Since  decreasing  KCl )  does  not  affect  the  submodularity  of  v.  we  may  assume 
fi(Cl)  =  KCl).  Then  it  is  clear  that  fi  <  v  on  simp(fl,  B).  The  Hahn-Banach  Theorem  therefore 
implies  there  exists  an  extension  l  of  n  such  that  l  ^  v  on  simp(  Cl .  A).  The  extension  of  fi 
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that  we  seek  is  now  given  by  fi(A)  =  l  (IA  ).  ■ 

Corollary  4.1 

Let  A  be  an  algebra  of  subsets  of  Cl.  and  let  v€sub(A)  with  v(0 )  =  0.  Let  (A,  :  t  ^0)  be  a 
family  of  sets  in  A  such  that  As  C A(  for  s  ,  and  such  that  v(A,  )  is  finite  for  all  t.  Then 
there  exists  /n€add(A)  such  that  fi^v  and  /jl(A,  )  =  v(A, )  for  t  ^  0. 

Proof 

Let  B  be  the  algebra  generated  by  the  sets  A, .  t  ^0.  and  let  n  be  the  additive  set  function 
on  B  such  that  fi(At  )  «  v(A,  )  for  t  ^0.  Let  B  €B:  then  B  has  the  form 

B  =  0  [Atl-ASl] 

i  =1 

where  <t,  <r,+1.  Let 
Bk  =  Lj  [A,~ASj] 

1=1 

so  that  Bn  =  B.  Then 

M**)=  ZMA'^piA,')]. 

i  =1 

Letting  B0  =  0.  we  will  show  by  induction  on  k  that  /j.(Bk  )^v(Bk  ),  O^k  .  The  claim  is 
trivial  if  k  =  0.  Now  suppose  fJ.(Bk  )^v(Bk  )  for  some  k  ^  0.  Since  Bk  is  the  the  intersection  of 
the  sets  Bk  +1  and  A,k .  and  A,k  is  the  union  of  these  sets,  the  submodularity  of  v  implies 

v(Bk  )+v(.Atk^l)<v(.Bk+1)+v(Atk+l). 


Therefore, 
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M**  «)=  Ms,  )+[^,t+lH(A,l+1)] 
^j/(st  H»Ult+l)-i<illt+l) 
^v(S, +1). 


and  the  induction  is  complete.  Therefore  MS  )  for  all  B  €B. 


Now  by  Theorem  4.1  there  exists  fx  €add(A)  such  that  p-  extends  n  and  p  ^v.  Since 
fi  {A,  )  =  v{A,  )  for  t  ^0.  the  proof  is  complete.  ■ 


4.3.  Product-Algebra  Net-works 

Let  ft  be  any  set  and  define  r  :  ft2-*  ft2  by  *  (cuo.^).  If  A  Cft2,  we  call  r(A  ), 

also  denoted  tA  .  the  transpose  of  A  .  A  is  said  to  be  symmetric  if  tA  —A  ,  and  antisymmetric 
if  A  Dr A  =0.  For  any  A  Cft2.  the  symmetric  part  of  A  is  defined  to  be  A  HrA  and  denoted 
cr(A  ).  and  the  antisymmetric  part  of  A  is  defined  to  be  A  OrA  and  denoted  cr  (A  ).  where 
A  =  ft2— A  .  Thus  A  =cr(A  )U<r  (A  ),  and  A  is  symmetric  if  and  only  if  A  =  oiA  ).  and 
antisymmetric  if  and  only  if  A  =  cr  (A  ). 

It  is  easy  to  see  that  the  class  of  all  symmetric  sets  is  an  algebra,  but  the  class  of  all 
antisymmetric  sets,  although  closed  under  finite  intersections,  is  not  closed  under  complementa¬ 
tion.  The  following  properties  of  <r  will  be  needed  later.  Their  proof  is  by  straightforward 
verification. 


Lemma  4.2 


Let  A  ,  B  €A. 


(a)  <r  (A )  =  cr  (B )  if  and  only  if  A  —B  and  B  —A  are  both  symmetric,  where 


AAB  =(A-B)U(B~A). 


v  s' 


A.' 

r.  *  -  -  .i..-. 
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(b)  cr  (A  U5  )  =  <7  (cr  (A  )U«r  (5)). 

A  class  A  of  subsets  of  Cl2  is  said  to  be  symmetric  if  tA  €A  whenever  A  €A.  By  a  prenet¬ 
work  we  mean  a  triple  (Q,VA),  where  V  is  an  algebra  of  subsets  of  SI.  and  A  is  a  symmetric 
algebra  of  subsets  of  Cl2  such  that  V  1xV26A  whenever  V  1,V26V. 

By  a  flow  on  A  we  mean  an  additive  set  function  X  :  A-*.R  such  that  X(A  )=0  whenever 
A  is  symmetric.  By  a  capacity  on  A  we  mean  an  additive  set  function  c  :  A -*R  such  that 
c(A)£0  whenever  A  is  symmetric.  Note  that  every  flow  is  also  a  capacity,  and  that  if  X  is  a 
flow  and  c  is  a  capacity,  then  c  —  X  is  also  a  capacity.  The  following  two  lemmas  give  some 
immediate  characterizations  of  flows  and  capacities. 

Lemma  43 

Let  X  :  A -*1R  be  additive.  The  following  are  equivalent. 

(a)  X  is  a  flow 

(b)  X(A  )+X(rA  )  =  0  for  all  A  €A 

(c)  X«cr  =  X 

Lemma  4.4 

Let  c  :  A-*J?  U{+oo}  be  additive.  The  following  are  equivalent. 

(a)  c  is  a  capacity 

(b)  c  (A  )+c  (tA  )^0  for  all  A  €A 

(c)  c- cr 

If  X  is  a  flow  on  A  and  Vj.V^V,  then  by  Lemma  4.3,  X(V2XV2)  =  —  X(V2xV1).  The 
quantity  X(V'1xV2)  represents  the  net  amount  of  some  commodity  moved  from  Vi  to  V2. 


We  define  a  product  -algebra  network  to  be  a  quadruple  (fl,V.A.A),  where  (fl.V,A)  is  a 


prenetwork  and  A  is  some  set  of  flows  on  A.  called  feasible  flows.  For  example.  A  may  be 
taken  to  be  the  set  of  all  flows  \  on  A  such  that  A^c  .  where  c  is  some  capacity  on  A. 

The  following  example  shows  how  a  classical  network  may  be  represented  by  a  product- 
algebra  network. 

Example  (Classical  Network) 

Let  A  be  a  finite  set.  We  call  the  elements  of  N  nodes  and  the  elements  of  A2  arcs.  Let 
b  =  ( b,j  :  (i )€N2)  where  b,j  is  real  and  btj  +  bji  ^  0.  Let  F  denote  the  set  of  vectors 
/  =  C/i>  :  (i  .  j  )€N2)  such  that  f  tj  +  f  -  0  and  f  u  ^bu  .  We  call  f  u  the  flow  in  arc  (i.j). 
bij  the  capacity  of  arc  (t  ,j ).  and  F  the  set  of  feasible  flows.  Note  that  —bjj  ^b^  for 
f  €F .  and  that  —bjt  can  be  strictly  positive. 

The  above  classical  network  gives  rise  to  a  product-algebra  network  (ft,  V.  A.  A)  where 
ft  =  N  .  V  and  A  are  the  power  sets  of  N  and  N2.  respectively,  and  A  is  the  set  of  flows  on  A 
such  that  A  ,  where  c  is  the  capacity  on  A  defined  by 

c  (A  )  =  X  b‘j  • 

U. JXA 

The  (ft.  V.  A.  A)  represents  the  classical  network  in  the  sense  that  a  flow  /  for  the  classical 
network  is  feasible  if  and  only  if  the  flow  A  on  A  defined  by 

AU)=  Z  fij 

V. jXA 

is  in  A. 

4.4.  Continuous  Feasible  Flow  Theorem 

Define  S  :  V-*A  by  8(V )  =  V  xV .  If  A  is  a  flow  on  A.  then  A.8  is  an  additive  set  function 
on  V.  since  A»8(V )  =  A(V  xV )  =  A (ftxV  ).  Furthermore.  A .8(ft)  =  0.  The  following  theorem 


generalizes  the  feasible  flow  theorem  oe  Gale  and  Hoffman  for  classical  networks  [Ford  and 
Fulkerson,  1962], 

Theorem  4.2 

Let  (ft,  V.  A)  be  a  prenetwork,  let  c  be  a  capacity  on  A,  and  let  y.  be  an  additive  set  func¬ 
tion  on  V  such  that  /t(ft)  =  0.  Then  there  exists  a  flow  X  on  A  such  that  A^c  and  \»8  —(j.  if 
and  only  if  /i^c.5. 

Remark 

The  additive  set  function  (i  is  often  called  the  demand  for  some  commodity.  Since  ) 
represents  the  net  amount  of  commodity  moved  into  V.  the  equation  \.8  =  n  implies  that  the 
demand  is  met. 

The  following  lemma  is  a  key  step  in  the  proof  of  Theorem  4.2. 

Lemma  4-5 

If  c  is  a  capacity  on  A,  then  c°<T  is  a  submodular  function  on  A.  It  follows  that  c°8  is  a 
submodular  function  on  V. 

Proof 

Let  A  ,  B  .  €A.  and  let  L  =  or  (A  U2? )  and  M  =  cr  (A  DB  ).  It  is  easily  verified  that 
L  C\M  =  <rG4  )n<r  (5)and£  UM  C<7  (A  )UcF(5  ). 

Now,  by  Lemma  4.2.  part  (b). 

(?  (L  UA/)=  <7  ((A  U3)U(A  HB )) 

=  <F  (A  UB ) 

=  <7  (<F  (A  )U<7  (5  )). 


Therefore,  by  Lemma  4.2.  part  (a). 


Z  =  <r  (A  )U<7  (B)  —  (L  UM) 


is  symmetric,  and  so  by  the  definition  of  a  capacity,  c  (Z  )^0,  implying 


c(L  U M)^c(a  U  )Uo-  (5)). 


Therefore,  using  the  additivity  of  c  , 

c(L)  +  c(M)  =  cU  U  M  )  +  c  U  r\M  ) 

<  c  (cr  (A  )UcF(B  ))  +  c(cr  {A  )Cta  (5  )) 

=  c(cr  (A  ))  +  c(<r  (B  )). 

This  shows  that  c.tr  js  submodular.  The  submodularity  of  c»8  now  follows  from  the  identity 
CoS(V)  =  c(VxV)»c.ar(nxV).  ■ 


Proof  of  Theorem  4.2 

Let  B  =  I  ClxV  :  V  €V}.  and  define  the  additive  set  function  fi'  on  B  by  £t'(  ftxV  )  =  fj.(V  ). 
Then  n'^c°<r  on  B  since  fi'iSlxV)  =  fi(V)^c°8(V)  =  c°<r  ({IxV).  By  Lemma  4.5.  c°(r  is 
submodular.  and  so  by  Theorem  4.1  there  exists  A€add(A))  which  extends  fi'  such  that 
A^c.<r  .  By  Lemma  4.4,  part  (c).  we  also  have  A^c.  Let  A  be  a  symmetric  set  in  A.  Then 
A(A  Xc.cr  (A  )  =  c(0)  =  0.  Since  A  is  also  symmetric,  A(A  )<0.  But  0  =  A((12)  - 
A(A  )+A(A  ),  and  so  A(A  )  and  A(A  )  are  both  zero.  Therefore  A  is  a  flow.  Finally,  if  V  €V. 
then  A°8(V )  =  A(V  xV )  =  A(ftxV )  =  fiXClxV  )  =  fi(V  ),  completing  the  proof.  ■ 


4.5.  Maximizing  Flow  into  a  Family  of  Sets 

Let  (O,  V,  A)  be  a  prenetwork,  let  c  be  a  capacity  on  A,  and  let  S  and  D  be  sets  in  V 
called  the  source  and  sink  .  respectively.  We  shall  consider  the  set  A  =  A(c  .  S .  D  )  of  flows  A 


on  A  satisfying  the  following  conditions: 


(Al)  X  ^  c 

(A2)  X.8(V)^0  ifVCS 

(A3)  X.S(V)<$0  if  V  CD 

Conditions  (A2).  (A3)  represent  the  restriction  that  commodity  can  be  generated  only  in  S  and 
absorbed  only  in  D.  If  S  =  D  =  d.  then  (A2),  (A3)  impose  no  restriction.  If  S  —  D  —  0, 
then  (A2).  (A3)  imply  X»8(V  )  =  0  for  all  V.  in  which  case  X  is  called  a  circulation .  Note  that 
if  c»8(V  )<0  for  some  V  C.S  .  then  A  is  empty.  We  shall  assume  A  is  nonempty. 

Let  B  €V.  We  are  concerned  with  the  problem 

maximize  {X°S(f?  ):  X€A|.  (PI) 

By  Theorem  4.2.  problem  (PI)  is  reduced  to 

maximize  {ft(B ):  }.  (P2) 

where  M  is  the  set  of  additive  set  functions  fi  on  V  such  that 

(Ml)  fi  <  c«8 

(M2)  ifVCS 

(M3)  m(VO  <0  if  VCD 

(M4)  fji(n)  =  0 

Lemma  4.6 

A  set  function  fi  is  an  element  of  M  if  and  only  if  /n  =  X.S  for  some  X€A.  It  follows  that 
if  fi  is  a  solution  to  problem  (P2).  then  any  X  such  that  ft  —  X«5  is  a  solution  to  problem  (Pi). 
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Proof 

If  A€A,  then  it  is  clear  that  p.  =  Conversely,  let  .  Since  /ui^c>8  and 

=  0,  Theorem  4.2  implies  there  exists  a  flow  A  such  that  A^c  and  A =8  =  p.  Now  condi¬ 
tions  (M2),  (M3)  imply  (A2).  (A3),  and  so  A6A.  ■ 

We  next  show  that  the  constraints  (Ml)  -  (M3)  may  be  replaced  by  a  single  constraint  of 
the  form  where  v  is  submodular.  If  V  lt  V2  are  disjoint  sets  in  V,  we  define  K(V  1(  V2)  to 
be  the  set  of  all  K  €V  such  that  V  lGK  and  V 2CK .  We  call  an  element  of  K(V j.  V 2)  a 
V i.  V 2  -  separating  cut .  For  V  € V  we  define 


v(V  )  =  inf  {c.S(A- )  :  K  €K(S  fl  V  .  D  n  V  )> 


Lemma  4.7 

The  set  function  v  is  submodular.  and  M  =  {/a€add(V)  :  fi^v.  //(ft)  =  0|. 


Proof 


Let  Vj,  V2€V,  and  let  €>0.  Then  by  definition  of  v  there  exists  Kj  €K(S  DV, .  D  flV, ) 
for  i  «  1.  2  such  that  c°8(if,  )<v(Vl  )+€.  Then 


a-iUa-2  €  K(sn(F7m777.£)  n(vluv2)) 


k jDAj  €  k(s  n 07771V 2 ;.  d  n(v,nv2)) 


and  so  v(V,iUV2)^c.8(Jir1U^2)  and  v(V1nV2)^Co8(A'inA'2).  Recall  from  Lemma  4.5  that 
c  »8  is  submodular.  Therefore, 


KV,UV2)  +  *»(v,nv2) 


\  -  v-v-i  •  v-y-y-v-y-I-'-  --y-y  ’•*  r 
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<  c£{KlUK2)  +  c  >8(^,0  jr2) 

<  c  »5(AT x)  +  c3{K2) 

<  i/(Vj)  +  v(V2)+2e. 

Since  the  inequality  holds  for  each  €>0,  v  is  submodular. 

Now  let  /u€add(V)  with  and  /u(fl)  =  0;  we  will  show  fj.£M .  Since 

V  €K(S  flV  .  D  DV ),  we  have  v(V  )^co8(V  ),  and  so  establishing  (Ml).  If  V  C.D  , 

then  0€K(S  C\V  .  D  C\V),  which  implies  v(V)^C.  Since  this  establishes  (M3).  If  V  CS  . 
then  0€K(5nV'.DnV),  which  implies  v(V  )^0.  Therefore  /x(V )  =  —  /jl(V  )^~  v(V  )^0, 
establishing  (M2).  Therefore.  fi£M  . 

Finally,  let  :  we  will  show  fj.^v.  Let  V  €V  and  let  K  €K(S  flV  .  D  flV').  Then 

K—  V  CS  and  V~ K  CD  .  Therefore  yu( K—  V  )^0  and  —K  )^0.  and  so 


)+M (K  -V  )-/*(  V  -K  )  =  nUC  Kc  .8U  ). 

Therefore  /x(  V  )^Co8(/C )  for  all  K  €K(S  (IV.DnV),  implying  /u(V  )^v(V  ).  ■ 

The  following  theorem  generalizes  the  max-flow  min-cut  theorem  for  classical  networks. 

Theorem  43  (Max-Flow  Inf-Cut) 

Let  B  €V  and  suppose  v(B  )  is  finite.  Then 


max  !\.$(5  )  :  A€A}  =  v(B  )  =  inf  (cSU  )  :  K  €K(S  DB  .DHB  )}. 


Proof 

Since  c.8(0)  =  0.  we  have  v(0)^O.  But  if  v(0)<O.  then  no  additive  set  function  is  dom¬ 


inated  by  v,  and  so  by  Lemma  4.7.  M  is  empty,  which  by  Lemma  4.6  implies  that  A  is  empty. 
But  since  we  assume  A  is  nonempty.  v(0)  =  0.  Similarly.  v(Cl)  =  0. 
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Therefore,  since  v  is  submodular,  Corollary  4.1  implies  there  exists  /x€add(V)  such  that 
)  =  v(B  ),  and  -  v(Cl)  =  0.  By  Lemma  4.7,  fi€M  .  Therefore  v(B  )  =  max 

\n(B  ) :  fx€M  >.  which  by  Lemma  4.6  is  equal  to  max  (X»5(5  ) :  A€A}.  ■ 

We  state  as  a  corollary  the  special  case  where  S  and  D  are  disjoint  nonempty  sets  and 
BCD.  This  case  more  closely  resembles  the  classical  max-flow  min-cut  theorem. 

Corollary  4.2 

Suppose  S  and  D  are  disjoint  nonempty  sets  and  let  B  €V  with  B  CD  .  If  v(B  )  is  finite, 

then 


max  |AoS(£  )  :  A€A|  =  v(B  )  =  inf  {08U: )  :  K  €K(S  .  B  )}. 

Remark 

Note  that  as  long  as  B  CD  ,  the  infimum  in  Corollary  4.2  does  not  depend  on  D  . 

The  following  theorem  shows  that,  given  an  increasing  family  of  sets  in  V.  there  exists  a 
single  A€A  which  simultaneously  maximizes  the  flow  into  each  set  in  in  the  family. 

Theorem  4.4 

Let  ( B , ,  t  £0)  be  a  family  of  sets  in  V  such  that  Bs  CB,  for  s  ,  and  suppose  v(B,  )  is 
finite  for  all  t .  Then  there  exists  a  flow  X*  6  A  such  that  for  all  t , 

\*°8(B, )  =  max  {X»S (Br ) :  X€A}. 

Proof 

By  Corollary  4.1,  there  exists  //€add(A)(V)  such  that  fi^v,  fi{£l)  =  KH)  =  0,  and 
fi(B,  )  =  v(B, )  for  all  t .  By  Lemma  4.7,  fi€M .  and  so  by  Lemma  4.6  there  exists  X*  €A  such 
that  A*  .8  =  /x.  Then  X *°8(B, )  =  (j.(Bc  )  =  v(B, )  =  max{A.8(5( ) :  A€A}.  where  the  last  equality 


is  from  Theorem  4.3.  ■ 


The  following  corollary  shows  that,  if  B,  CD  for  all  t.  then  the  problem  of  finding  the 
flow  X*  of  Theorem  4.4  can  be  solved  by  first  finding  a  flow  which  maximizes  k>8(Bt  )  for  all 
t  and  then  finding  a  flow  on  a  modified  network  which  maximizes  \-8(B,  )  for  all  t  ^r. 

Corollary  4.3 

Suppose  S  and  D  are  disjoint  nonempty  sets,  and  let  ( B ,  :  t  ^0)  be  as  in  Theorem  4.4 
except  that  B,  CD  for  all  t.  Fix  r> 0,  let  Aj  =  A(c  ,  S .  BT).  and  let  XxGAx  maximize 
{Xo8(5f ) :  X€ Aj)  for  all  t<r.  Now  let  A2  =  A(c S .  D  ),  where  c '  *  c  —  Xj,  and  let  X2€A2  max¬ 
imize  {X°S(Z?,  )  :  X€A2}  for  all  t^r.  Then,  letting  X*  =  X!  +  X2.  X*  €A  =  A(c  .  S  ,  D)  and  X* 
maximizes  {X«8(21,  )  :  X€A}  for  all  t  ^  0. 

Proof 

It  is  easily  verified  that  X*  €A.  First  let  t  ^r.  Then  B,  CBrCD .  Therefore,  since  the 
infimum  in  Corollary  4.2  does  not  depend  on  D  . 

KV8(B,  )  =  max  {X°8(2?r  )  :  X€Axl  =  max{X.8(5,  )  :  XGA|. 

Now  since  X2°8(5,  )>0.  we  have  \*°8(B,  )^\v8(Bt  ),  and  so  X*  maximizes  {X°8(2?(  )  :  X6A}  for 
all  t  <r. 

Since  XxGAx,  Xj-SCV )  =  0  whenever  V  is  disjoint  from  both  S  and  B T.  Therefore,  by  the 
additivity  of  Xr8,  we  have  Xx.8(A )  =  Xx»8(5r)  for  K  €K(S  ,  BT).  So  for  t  ^r. 


\*°8(Bt  )  =  Xj.SfB,  )  +  X2"8(5r ) 

=  XrS(flT)  +  inf  {c  '°8(K  )  :  K  €K(S  ,  Bt  )| 
=  inf  {(Xx+c ') ,8(A  )  :  K  6K(S  .  B,  )( 

=  max  (X«8(5, ) :  X€A} 


where  the  last  equality  follows  from  Corollary  4.2  since  Xj+c  '  =  c  .  ■ 
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Remark 

If  B,  .  t  ^0  takes  on  only  finitely  many  values,  then  Corollary  4.3  implies  that  X*  can  be 
obtained  by  solving  a  finite  number  of  maximum  flow  problems. 

Although  Corollary  4.3  decomposes  the  problem  of  finding  k*  into  two  subproblems,  these 
subproblems  are  dependent  in  the  sense  that  one  must  be  solved  before  the  other  can  be  solved. 
The  next  theorem  shows  that  if  there  exists  a  cut  which  minimizes  {c°8(K )  :  K  €K(S.5X)}, 
then  the  problem  of  finding  X*  decomposes  into  two  independent  subproblems.  We  first  need 
the  following  lemma. 

Lemma  4.8 

Let  S  .  D  .  and  B  be  as  in  Corollary  4.2.  Suppose  there  exists  Km  €K(S  .  B  )  such  that 

c&(Km  )  =  v(B )  =  inf  |c.S(Ji: ) :  K  €K(S  .  B  )}. 

Then  any  flow  X*  €A  maximizing  {X.8(i?  )  :  X€A)  must  satisfy 
X*,8Um)  =  c.8Um). 

Proof 

Since  Km  €K(S  M ).  we  have  Km—BCS,  and  so  X* °8(Km  —B  )^0.  Therefore. 
X*-8Cfi  )^X*»8(Arm  ).  But  Corollary  4.2  implies  X*.8(5  )  =  v(B  )  =  c>8(K„  ).  and  so 
k*>8(JCm)^c°8(Km).  The  lemma  now  follows  from  X*  .  ■ 

Theorem  4.5 

Let  S.D.  and  (B, .  t  ^0)  be  as  in  Corollary  4.3.  Fix  r> 0,  and  suppose  there  exists 

Km  €K(S .  Br)  which  minimizes  {c.8(A" )  :  K  €K(S ,  5T)}.  Let  Xj€A  maximize  {X»8(5, ) :  X€A} 

for  all  t  <r,  and  let  X2€A  maximize  {X«8(5, )  :  X€A)  for  all  t  ^r.  Now  define  the  flow  X*  by 
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X*  (A  )  =  X,(A  C\K }  )+X2(A  C\K})  +  c{A  n8XTT)-cU  DSKt). 

Then  X*  €A,  and  A*  maximizes  iX»8(.fi, ) :  X6Al  for  all  t  ^0. 

Proof 

It  is  easily  verified  that  X*  €A.  Since  Xt  and  X2  both  maximize  {X=8(fiT)  :  X€A},  Lemma  4.8 
gives 

AV8{KT)  =  X2*S(if  T)  =  c°8(1Ct). 

Therefore,  since  X^c  and  Xj.  X2.  c  are  additive,  we  have  Xt(A  HS/f,.)  =  X2U  flSXf)  * 
c(Afl8iifT)  for  all  A  €A.  Since  X;  (rA  )  =  — X;  (A  ),  i— 1,2.  we  also  have  Xj(A  H8Ar)  « 
X2(A  HSif  T)  =  — c(tA  n8^Tr).  Substituting  these  equations  into  the  definition  of  X*  .  we  see 
that 

X**8(V  )  =  X1»8(V  )  if  VdKr 

and 

A*  *8(  V  )  =  X2«8(V  )  if  V  CK  T  . 

Now  if  t  ^t,  then  B,C.Kr.  and  so  X*"8(5,  )  =  X^Cff,  ).  Therefore,  X*  maximizes 
{X*8(B,  )  :  X€A)  for  all  t  <t. 

Now  suppose  t  ^r.  Then 

A*  MB,  )  =  AV8(B,  r\KT)+A2MB,  HK  T). 

Let  i  «  1  or  2.  Since  Xj*8(5T)  -  Xj^(ifT)  and  X,-.S(V  )^0  for  V  CS  ,  we  have  Xj*8(V )  =  0  for 
VC Kt—Bt,  implying  X,*8(5T)  =  X,-*8(.B(  C\K T).  Therefore,  since  Xr8(5r)  ■  c*8(KT)  « 
X 2*8(5 T),  we  have  AV8(B,  C\K T)  -  X2*8(5,  C\Kr).  Therefore. 
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5.  MAXIMUM  FLOWS  IN  DYNAMIC  NETWORKS  WITH  CONTINUOUSLY  VARYING 
CAPACITIES 

5.1.  Introduction 


In  this  chapter,  we  consider  a  dynamic  network  of  the  type  in  Chapter  3,  but  where  the 
link  and  storage  capacities  are  no  longer  required  to  be  piecewise  constant.  Recall  that  in 
Chapter  3,  the  dynamic  max-flow  m in-cut  theorem  followed  easily  from  the  static  max-flow 
min-cut  theorem  for  finite  networks.  This  is  not  the  case  for  the  more  general  model  con¬ 
sidered  in  this  chapter.  Instead,  we  will  use  the  theory  of  product-algebra  networks  presented 
in  Chapter  4  to  prove  a  dynamic  max-flow  min-cut  theorem  which  generalizes  the  one  in 
Chapter  3.  This  theorem  is  similar  to  the  one  given  in  Anderson,  et  al.  (1982),  but  their 
theorem,  as  stated  is  incorrect.  However,  after  I  informed  these  authors  of  their  error,  they 
have  proved  (but  not  yet  published)  a  modified  version  of  that  theorem.  Their  method  of  proof 
involves  a  dynamic  version  of  the  flow  augmentation  technique  of  Ford  and  Fulkerson  (1962) 
and  is  therefore  very  different  from  the  method  used  in  this  chapter. 


Recall  that  an  important  discovery  in  Chapter  3  was  the  existence  of  a  flow  which  is  r- 
maximum  for  all  r.  We  also  extend  that  result  to  the  model  of  this  chapter.  It  is  hoped  that 
these  results  will  contribute  to  the  understanding  of  flow  optimization  in  dynamic  networks 
and  may  lead  to  new  methods  for  computing  optimal  dynamic  flows. 


5.2.  Problem  Formulation 

Let  N  be  a  finite  set  of  nodes,  including  a  source  s  and  a  sink  d  ,  and  fix  T  >0.  By  a  link 
we  mean  an  element  of  N2.  For  each  link  (t )  we  assign  a  nonnegative  bounded  Lebesgue- 
measurable  function  bi}  on  (0.  T  J  giving  its  capacity  at  each  time,  and  for  each  node  j  €N  we 
assign  a  nonnegative  lower-semicontinuous  function  a;  on  (0,  T  ]  giving  its  storage  capacity  at 
each  time.  Letting  b  -  {btj:i  ,j  &N )  and  a  —  (a, :  j€N).  we  define  the  continuous-time 


dynamic  network  D  =  (N  .  b  .  a  ,  s  .  d  ). 

We  define  a  feasible  flow  for  D  to  be  an  assignment  u  =  (u:J :  i  ,j  €N  )  such  that  each  u:J  is 
a  Lebesgue-measurable  function  on  (0,  T  ]  and  the  following  conditions  are  satisfied  for  all 
i  .j  €jV  and  all  t  €(0.  T  ]: 


Utj  {t  )  =  -Uj,  ( t  ) 

(5.1) 

ut j  (f  )  <  b,j  ( t  ) 

(5.2) 

0  <  Xj  (r )  <  a}  (r ) 

(5.3) 

where  x  =  (xy :  j  €N  ),  called  the  storage  function  due  to  u  .  is  defined  by 
x,  (f  )  =  f  £  u,j  (t  *)  dt \  i€N—  {r  .d  } 

0  i  €  ,V 

xs  (t )  =  xd  (r )  =  0 

The  set  of  all  feasible  flows  for  D  will  be  denoted  U .  We  interpret  as  the  net  rate  at  which 
commodity  is  moved  from  node  t  to  node  j  .  and  xy  (t  )  as  the  amount  of  commodity  stored  at 
node  j  at  time  t.  Note  that  (5.1)  implies  uyy(t)= 0  and  that  (5.1),  (5.2)  imply 
— bj,  ( t  )^uiy  ( t  )^bij  ( t ).  Also  note  that  we  are  taking  the  initial  storage  to  be  zero.  However, 
as  seen  in  Chapter  3.  this  restriction  does  not  prevent  us  from  allowing  nonzero  initial  storage 
in  the  minimum-delay  problem  considered  in  that  chapter. 

When  sets  of  nodes  are  used  as  subscripts,  the  summation  convention  applies.  Thus 
“.-t  s  (t  )  =  £  L  uu  )• 

,iA  iSB 


For  each  r€(0,  T]  we  define  vT  on  U  by 
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l 


vT(u  )  =  f  uN  d  (l  )  dt. 
o 

This  represents  the  total  amount  of  commodity  reaching  the  sink  before  time  T.  By  a 
r— maximum  flow  we  mean  a  flow  u  61/  such  that  vT(u  )^vT(u  ')  for  all  u  ‘€U  . 

5.3.  Dynamic  Max-Flow  Min-Cut  Theorem 

In  this  section  we  present  a  dynamic  version  of  the  max-flow  min-cut  theorem  and  derive 
necessary  and  sufficient  conditions  for  a  flow  in  U  to  be  a  r-maximum  flow.  We  also  prove  the 
existence  of  a  flow  u  €£/  which  is  r-maximum  for  all  r€( 0.  T ]. 

We  define  a  dynamic  separating  cut  for  the  network  D  to  be  a  set-valued  function 
G  :  (0.  T  ]-*2‘v  which  is  left-continuous  with  at  most  a  finite  number  of  jumps,  and  satisfies 
d  €G  (r  )  and  s  €G  (t )  for  all  t .  Note  that  this  definition  is  different  from  the  one  in  Chapter  3. 
The  set  of  all  dynamic  separating  cuts  will  be  denoted  G. 

If  G  is  a  dynamic  cut  and  _/  €  N  .  we  define 
Gj  =  \t:j€G(t) ). 

Then  G j  has  the  unique  representation 

Gj  =  ()  (<Xjk.P)k]  (5.4) 

lc  =1 


where 

0  ^  a jk  <  (ijit  <  ctjjt+i  <  Pj  jt  +i  ^  T ,  1  ^  k  ^  Zj  —1. 

Thus  no  interval  in  the  union  is  empty,  and  the  intervals  have  disjoint  closures.  For  each 


r  €  (0,  T  ]  we  define  w  f  on  G  by 


w  t(G  )  —  f  bG  )  G  ,(r )  dt 
o 

+  £  L  ay  (<*/*  )• 

y  €  JV  <*jA  <  r 

We  define  a  t— minimum  cut  to  be  a  cut  G  €G  such  that  wT(G  )  ^  wT(G  ')  for  all  G  ’€G. 

Theorem  5.1 

For  each  r€(0.  T  ], 

max  {vT(u  ):  u  €G  |  =  inf  (wr(G  ):  G  €G}. 

Furthermore,  if  G  is  a  r-minimum  cut.  then  a  flow  u  €U  is  a  r-maximum  flow  if  and  only  if 
the  following  conditions  hold: 


Uij(t  )  =  bjj(t).  (i  ,j  )€G  (f  )xG  (f  ).  a.e.  t  €(0,t],  (5.5) 

Xj(ajk  )  —  aj(ajlt  ).  7  <  T.  (5.6) 

x;(0j*)  =  O.  <  r.  (5.7) 

*;(r)  =  0.  j  €G  (t).  (5.8) 


where  a;i  and  0 ]k  correspond  to  the  representation  (5.4)  for  G, . 

The  proof  of  Theorem  5.1  will  follow  three  lemmas. 

Lemma  5.1 


Let  u  €f/ .  G  €G.  and  r€(0.  T  ].  Then 
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vT(u)  =  f  uoU)fOU)(t)  dt  -  £  Xj(t ) 

0  y€G(r) 


+  £  £  xj  (otjk )  -  £  xj  ifijk ) 

oryt<T  I3jk  <T 


It  follows  that  vT(u  )  ^  wr(G  ).  and  that  vT(u)  =  wr(G  )  if  and  only  if  conditions  (5.5)  -  (5.8) 


of  Theorem  5.1  hold. 


Proof 


Since  this  proof  is  similar  to  the  one  for  Lemma  3.1.  we  omit  it. 


In  view  of  Lemma  5.1.  proving  Theorem  5.1  amounts  to  showing  there  exists  a  flow  u  €£/ 


such  that  vT(u  )  =  inf  |wr(G  ):  G  €G}.  We  will  do  this  by  constructing  a  product-algebra  net¬ 


work  which  represents  the  dynamic  network  and  then  applying  Corollary  4.2. 


Let  I  denote  the  class  of  intervals  of  the  form  (t  j.  t2],  O^C  l^t24T .  Let  Cl  ~  N  x(0.  T ]. 


and  let  V  denote  the  smallest  algebra  of  subsets  of  ft  containing  {i  }x/  for  each  i  £N  and  I  61. 


For  convenience  we  shall  let  /,  denote  (t  \xl .  Let  A0  be  the  class  of  all  sets  of  the  form  I,  X/ j 


with  i  .j  €iV  and  I  J  €1.  and  let  A  be  the  algebra  generated  by  A0.  It  follows  that  A  contains 


VjXV^  for  all  V,.  V2€V. 


Define  c  on  Ao  as  follows: 


c  (It  XJj )  =  f  bij(t)dt.  i^j- 


+00.  inj*0. 

c(IjXJj)=  cij(t).  max/  =inf/  =t 
0,  otherwise. 


where  we  take  at  ( t  )=0  for  j  —s  ,d  .  Thus,  if  I  and  J  are  disjoint  but  are  "touching"  at  t ,  then 
c  (Ij  Xjj  )  is  equal  to  either  <z;  (f  )  or  0  depending  on  whether  J  is  to  the  right  or  to  the  left  of 


r«t- 
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I .  Note  that  c  (/,  Xjj  )=0  if  i^j  and  I  fl  /  =0.  It  is  easy  to  verify  that  c  is  additive  on  A*. 
Since  Ao  is  a  semialgebra  which  generates  A.  it  follows  [Royden.  1968]  that  each  set  in  A  is  a 
finite  union  of  disjoint  sets  in  Ao  and  that  c  has  a  unique  additive  extension  to  A.  given  by 


m 

c  (A  )  =  £  C  (Ai  ) 
i=l 

where  {A  1(  .  .  .  ,  Am  }  is  any  finite  partition  of  A  with  Ai  €  Ao- 

Letting  5  =  {s  }x(0,  T]  and  D  =  {d  |x(0.  T],  we  define  the  product-algebra  network 
N  =  (ft,V.  A.  A),  where  A  =  A (c.S.D)  is  defined  in  Section  4.5.  Recall  that  a  V  v  V  2~ 
separating  cut  for  N  is  a  set  K  in  V  such  that  V  jCAT  —Cl—K  and  V 2CK .  We  let  K(V  ,.V 2) 
denote  the  class  of  all  such  cuts.  If  A  is  such  a  cut.  we  define  K  (r  )  =  { j  €N :  (j  X  )ZK  |.  For 
each  t€(0.  T]  we  define  the  set  Dr  —  \d  }x(0.r].  If  K  €K(S  DT).  then  since  DTC.K .  we  have 
d  €  A  (*  )  for  all  t  <T. 

Lemma  5.2 

If  K  €K(S  ,Dt)  and  G  is  any  cut  in  G  such  that  G(t)~K(t)  for  all  t  then 
wt(G)  ^c(AxA). 

Proof 

This  is  easily  verified  from  the  definitions  of  wr(G  )  and  c  .  ■ 

Lemma  5J 

Let  \€A.  There  exists  a  dynamic  flow  u  €U  such  that  for  i  ,j  €  A  and  /  61. 

uu  (t  )  dt  =  X(/j  XI j  ).  (5.10) 

and  for  j  €N  and  t  €(0,T  J. 


Xj  (f )  =  X((0.r  ]j  x(f  ,T  ] j ). 


(5.11) 


It  follows  that  for  each  tG(0,  T\.  vt(u  )  =  X(25  TXDT). 

Proof 

For  i  ,j  €N  .  define  h>j  on  I  by 

fjL.j  U)  =  \UiXlj). 

Since  X  is  additive  with  XfViXVj)  =  —  X(V2xV,1).  Hij  is  additive  with  h,j  ——  Hlt .  Therefore, 
from  X^c  we  get 

-/  bj,(t  )  dt  </xi;(/)<  f  b,j  (t )  dt. 

1  i 

Using  basic  measure  theory,  this  implies 

HU  (/  )  =  f  Ujj(t)dt,  I  €1. 

7 

for  some  measurable  function  utJ  such  that  —bj,  ^bij .  Since  Hij  =~Hj,  •  we  may  redefine 
(utj:  i  ,j  €N)  on  a  set  of  zero  measure  so  that  utJ(,t  )=— u^(r  )  for  all  r  €(0.  7*].  Then 
u  -(.Uij :  i  ,j  )  satisfies  (5.1)  and  (5.2). 

We  let  I,  denote  O—  /, ,  not  to  be  confused  with  (  /  )>  =((0,  T  ]— /  )* .  Since  X  is  additive. 


X(  /  j  Xlj)=  £  X(7,  XI j  )  +  £  X((  /  );  XI j  )  +  X((  /  );  X/  ^  ). 

I  €  ,V  i  €  V 


Since  c  (/j  X7,  )=0  if  i  and  I  C\J  =0,  the  second  summation  in  the  above  equation  is  zero. 


and  so 
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X(  /  j  X/,  )  =  £  (/  )  -  A(7,  x(  /  ), ).  (5.12) 

UN 

Since  A(V  XV  )=0  for  V  Cft— (S  U D  ),  the  left  side  of  (5.12)  is  zero  for  t  ^ s  ,d  .  Furthermore, 
if  /  =(0 ,t  ].  the  summation  in  (5.12)  is  equal  to  x,  ( t  )  for  j  r*s  ,d  .  Therefore. 

Xj  (t )  =  \((0.r  x(r  J], ),  ;  ,d. 

Using  A^c  and  A(V!XV2)  =  — X(V2xV'1).  the  right  side  of  the  above  equation  is  bounded 
between  0  and  a,(t)  if  0 <t  <T  and  is  equal  to  zero  if  t  =0  or  t  -T .  Therefore  x,  ( t  )  satisfies 
(5.3).  proving  that  u  €U . 

Finally. 

T 

vT(u)  =  f  uN  d(t)dt  -  Z  Vjd  ((O.r])  =  TxDr). 

0  ;€.V 

where  the  last  step  holds  because  Z)T=(0,r]</  and  the  last  term  on  the  right  side  of  (5.12)  is  zero 
for  j—d.  ■ 

Remark 

If  we  visualize  N  x(0.  T  ]  as  a  graph  where  the  horizontal  axis  represents  N ,  then  (5.10) 
implies  the  vertical  part  of  X  represents  /  ,  and  (5.11)  implies  the  horizontal  part  of  X 
represents  the  storage  x  .  It  can  be  proved  that  (5.10)  determines  a  one-to-one  correspondence 
between  A  and  U  .  but  we  do  not  require  this  fact. 

Proof  of  Theorem  5.1 

Corollary  4.2  (the  max-flow  m in-cut  theorem  for  product-algebra  networks)  with  B  =DT 
implies  there  exists  a  flow  A€A  such  that 
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k(D  txDt)  =  inf  {c  (AT  XAT ):  K  €K(S  ,  DT)\  2  inf  (wT(G  ):  G  €G|. 

where  the  last  inequality  follows  from  Lemma  5.2.  Now  by  Lemma  5.3  there  exists  a  flow 
u  €U  such  that 

v T(u  )  =  \{DTxDT )  Z  inf  {wr(G  ):  G  €G}. 

The  theorem  is  now  immediate  in  view  of  Lemma  5.1.  ■ 

Theorem  5.2 

There  exists  a  flow  u'  €(J  which  is  r-maximum  for  all  r€(  0,  T  ],  i.e.,  such  that 
vT(u'  )  =  max  {vr(u  ):  u  €U  }  for  all  r€(0.  T ]. 

Proof 

Theorem  4.4  with  B,  =  Dt  implies  there  exists  a  flow  \*  6  A  such  that  for  each  tG(0,  T  ]. 

( D  TXD  t)  =  max  {X(D  TxD  T):  X  6  A} 

=  inf  {c  ( K  XK  ):  K  6K(S  .  D T)\ 

^  inf  {wr(G  ):  G  €G| 

=  max  {vT(u  ):  u  €G  }, 

where  the  last  three  steps  follow  from  Corollary  4.2.  Lemma  5.2,  and  Theorem  5.1.  respec¬ 
tively.  Now  by  Lemma  5.3  there  exists  a  flow  u  such  that  wT(u‘  )  =  X‘  (D  TxD T)  for  all 
r€(0,  r ].  and  the  theorem  is  proved.  ■ 
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6.  MAXIMUM  FLOWS  IN  VECTOR  SPACE  NETWORKS 

6.1.  Introduction 

In  this  chapter,  we  generalize  the  theory  of  abstract  networks  presented  in  Chapter  4.  We 
still  consider  a  flow  to  be  an  additive  set  function  on  an  algebra  of  subsets  of  Cl2.  but  we  no 
longer  require  the  capacity  to  be  an  additive  set  function.  Instead,  the  flows  are  identified  with 
linear  functionals  and  they  are  constrained  by  a  sublinear  functional.  Using  the  Hahn-Banach 
Theorem,  a  generalization  of  the  feasible  flow  theorem  and  of  the  max-flow  min-cut  theorem 
are  proved  in  Section  6.3  under  a  certain  condition  on  the  constraining  sublinear  functional. 
This  condition  involves  the  concept  of  vertical  additivity,  which  is  presented  in  Section  6.2.  In 
Section  6.3,  we  consider  the  special  case  (polymatroid  networks)  where  the  flows  are  con¬ 
strained  by  a  submodular  set  function.  As  an  application,  we  consider  a  class  of  finite  net¬ 
works  which  are  related  to  the  polymatroid  networks  of  Lawler  and  Martel  (1982).  In  Section 
6.4.  we  generalize  even  further,  taking  flows  to  be  linear  functionals  on  an  arbitrary  vector 
space,  and  we  prove  another  version  of  the  feasible  flow  theorem. 

6.2.  Vertically  Additive  Sublinear  Functionals 

Let  Cl  be  any  set  and  let  A  be  an  algebra  of  subsets  of  Cl.  We  let  add(A)  denote  the  set  of 
all  additive  set  functions  v  on  A.  sub(A)  denote  the  set  of  all  submodular  functions  v  on  A 
such  that  v(0)—O.  and  simp(0.  A)  denote  the  linear  space  of  all  simple  functions  on  (CIA.). 
See  Section  4.2  for  the  definitions  of  additive  and  submodular  set  functions  and  simple  func¬ 
tions.  We  let  IA  denote  the  indicator  function  of  the  set  A  .  If  pi€add(A)  and  /  €simp(n.  A) 
with 

/  = 

t  =1 


we  define 
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M(/  )  =  f  fdn  =  £a,  MUi  )• 

<  =1 

Recall  from  Section  4.2  that  the  mapping  p-+p  defines  a  one-to-one  correspondence  between 
add(A)  and  the  set  of  all  linear  functionals  on  simp( ft.  A). 

Recall  that  a  functional  0  on  simpCfl,  A)  is  said  to  be  sublinear  if  the  following  conditions 
are  satisfied  for  all  simple  functions  /  ,g  : 

)  =  a0(/  ).  a  ^  0.  (nonnegative  homogeneity) 

0(/+g)^0(/)  +  0(g ).  (subadditivity) 

Given  any  subset  M  of  add(A).  we  define  the  support  functional  of  Af  to  be  the  function  <f>  on 
simp(  fi .  A)  given  by 

0(/  )  =  sup  {/*(/  ):  }. 

It  is  clear  that  the  support  functional  of  any  subset  of  add(A)  is  sublinear  if  it  is  finite. 

We  are  interested  in  subsets  M  of  add(A)  that  have  the  form 

M  -  {/i€add(A):  p.  <  <f>)  (6.1) 

for  some  sublinear  functional  <f>  on  simp(  fl ,  A).  It  is  easy  to  show  that  a  set  M  has  this  form 
if  and  only  if  it  is  the  intersection  of  an  arbitrary  number  of  closed  half-spaces 
i/x€add(A):  f  )  <  a}  and  has  a  finite  support  functional.  For  example,  if  v  is  any  (not  neces¬ 
sarily  additive)  set  function  on  A.  then  the  sets  (yu€add(A):  /x(0)=KO)i  and 

|/i€add(A):  each  have  the  form  (6.1)  for  some  sublinear  functional  0. 


Let  0  be  a  sublinear  functional  on  simp(0.  A)  and  let  /  €simp(Q.A).  Then  the  Hahn- 


Banach  Theorem  [Royden.  1968]  implies  there  exists  /*€add(A)  such  that  ft  ^  <f>  and 
fiif  )  =  0(/  ).  Equivalently. 

0(/  )  =  max  {fi(f  ):  /x€ M  > 

where  M  is  given  by  (6.1).  Therefore,  if  M  has  the  form  (6.1)  for  some  sublinear  functional 
0.  then  0  must  be  the  support  functional  of  M . 

We  next  show  how  a  submodular  set  function  naturally  determines  a  sublinear  functional. 
If  i/€sub(A)  we  define  v  on  simp(Q.  A)  recursively  by 

QO 

v(f  )  ~  f  v{(u:  f  (ot)  ^  a)  d  a.  f  ^  0.  (6.2) 

o 

v(f  +«)  =  ?(/ )  +  ai/(0).  /  ^  0.  a€JR.  (6.3) 

Since  (6.2)  implies  v(f  +a)  »  v(f  )+av(ft)  whenever  /  +a  is  nonnegative,  v  is  well  defined. 
Note  that  if  /z€add(A).  then  fi  —  fi.  Thus,  the  definition  of  v  generalizes  the  concept  of  integra¬ 
tion. 

Lemma  6.1 

Let  v€sub(A).  Then  v  is  sublinear. 

Proof 

By  (6.3)  it  suffices  to  show  v  is  sublinear  on  the  nonnegative  simple  functions.  This  is 
implied  by  Theorem  54.1  of  Choquet  (1953).  ■ 

A  sublinear  functional  <f>  on  simp(fl.  A)  will  be  called  vertically  additive  if  0  =  v  for 
some  v€sub(A).  Clearly  any  linear  functional  is  vertically  additive.  An  example  of  a  verti¬ 
cally  additive  functional  which  is  not  linear  is  0(/  )  -  max  {/  (to):  o>€  fl).  Defining  v  by 
v(0)=O  and  v(A  )=1  for  all  nonempty  A  .  it  is  clear  that  ?=0.  The  following  lemma  provides 


two  useful  methods  for  determining  when  a  sublinear  functional  is  vertically  additive. 


Lemma  6.2 

Let  0  be  a  sublinear  functional  on  simp(0.  A).  define  the  set  function  v  on  A  by 
u(A  )  *  4>Ua  )•  and  suppose  that  0(/  +a)  «  0(/  )  +  av(Cl)  for  all  /  €simp( SI .  A)  and  «€^. 
Then  the  following  are  equivalent. 

(a)  0  is  vertically  additive. 

(b)  If  A  . . An  is  an  increasing  collection  of  sets  in  A.  then 


1=1  1=1 

(c)  If  A  ! . An  is  an  increasing  collection  of  sets  in  A.  then  there  exists  /x€add(A)  such 

that  (i  <  0  and  /oi(A,  )  =  v(A,  )  for  all  i . 

Proof 

That  (a)  implies  (b)  follows  easily  from  the  fact  that  0  =  v  if  0  is  vertically  additive. 
Now  suppose  (b)  holds  and  let  A  .  . . ,  An  be  an  increasing  collection  of  sets  in  A.  By  the 
Hahn-Banach  Theorem,  there  exists  /i€add(A)  such  that  n  ^  0  and  *  0(2^, )• 

£/x(Aj  )  -  ).  Since  implies  n^v,  we  have  n(At )  =  v(A, )  for  all  i .  Therefore,  (b) 

implies  (c). 

Finally,  suppose  (c)  holds,  and  let  /  €simp(  fl .  A)  be  nonnegative.  Then  /  can  be 
expressed  as 

/  *  z  «</*,. 

i=l 

where  A  i . An  is  an  increasing  collection  of  sets  in  A.  Thus,  by  (c)  there  exists  M^addlA) 


such  that  (i  ^  4>  and  n(Ai )  =  v(A,  )  for  all  i .  Therefore. 


M/  )  =  Z  “i  V^A‘  )  =  *'(/)• 

i  =1 


and  since  n^<t>.  we  have  K  f  )<<£(/  ).  But  the  sublinearity  of  <f>  implies 


<t>(f  )  ^  Z  a<  p(Ai  )**(/). 

i=! 


and  so  <{>  =  v.  ■ 

The  following  lemma  states  a  useful  consequence  of  a  sublinear  functional  being  vertically 
additive. 

Lemma  63 

Let  v€sub(A).  Then 

{>x€add(A):  /jl  ^  v]  =  |/x€add(A):  fi  ^  p.  =  p(SI)}. 

In  particular,  v  is  the  support  functional  of  the  set  on  the  right. 

Proof 

First  suppose  /i^p.  Then  clearly  Also.  — n(Sl)  •  fiC—l)  ^  ?(—  1)  -  — Kfl).  implying 

=  v(H).  Conversely,  suppose  and  m(O)  =  v(O).  Let  /  be  a  nonnegative  function 
in  simp(0.  A)  and  let  o€l? .  Then 

M/  +«)  =  m(/  )  +  ^  )  +  av(fl)  =  P(/  +a). 

Therefore. ■ 


Remarks 


(1)  Suppose  we  are  only  interested  in  nonnegative  additive  set  functions,  i.e..  we  are 
interested  in  finding  the  support  functional  of  the  set 

A/,  =  i/i€add(A):  0  <  /i  ^  v,  /z(ft)  -  v(ft)} 

For  example,  if  A  is  a  O’ -algebra,  v(  ft)=l.  and  An  -*0  implies  viA„  )-*0.  then  M  i  is  a  set  of 
probability  measures.  First  suppose  v  is  increasing,  i.e..  A  CB  implies  v(A  )£v(B  ).  Then  n'&O 
is  implied  by  fi^v  and  for  if  fi{A  )<0  for  some  set  A  then  n(Cl~-A  )  >  m(SI)  - 

v(Cl)  ^  v((l—A  ).  Thus  Lemma  6.3  implies  that  v  is  the  support  functional  of  M x.  If  v  is  not 
increasing,  then  it  can  be  modified  to  be  increasing  by  defining 

i >'(A  )  =  inf  v(B  ). 

BOA 

Then  since  -  v(fl).  and  since  if  and  only  if  Lemma  6.3  implies  that  v' 

is  the  support  functional  of  M 

(2)  Suppose  further  that  we  do  not  wish  to  require  /*(  O  )  -  v(  ft  ).  i.e.,  we  wish  to  find  the 
support  functional  of  the  set 

M 2  ~  {/u€add(A): 

For  example,  if  ft  is  finite,  then  Af2isa  polymatroid  [Hassin.  1982],  [Schrijver.  1984].  If  v  is 
increasing,  then  the  functional  v  defined  by 

oo 

~v  if  )  =  ?(/+)  =  f  v{a >:  /  M>«|rfo 
0 

is  sublinear  since  ?(( /  +g  )+)  ^  v(f  ++g  +)  ^  v(  f  +)  +  ?(g  +).  We  will  show  that  v  is  the  sup- 


port  functional  of  Af  2.  i.e..  that 


M2  -  {^t€add(A):  p  <  v  ). 

First,  it  is  clear  that  p  ^  v  if  and  only  if  p  <  v  .  Now  suppose  p  ^  v  and  let  A  €A.  Then 
— p(A  )  =  p(—IA  )  ^  v  (~IA  )  -  0.  Therefore  p  ^  0.  ■ 


6.3  Product-Algebra  Networks 

By  a  product  -algebra  network  we  mean  a  quadruple  ( Cl .  V.  A.  A),  where  V  is  an  algebra 
of  subsets  of  Cl.  A  is  an  algebra  of  subsets  of  V  such  that  V'jXV'jCA  whenever  V!.  V2€V,  and 
A  is  a  subset  of  add(A)  of  the  form  (X€add(A):  X  ^  p).  where  p  is  a  sublinear  functional  on 
simp(fl2.  A).  For  example,  as  discussed  in  Section  6.2.  A  may  have  the  form 
|X€add(A):  0  <  X  <  c  }  where  c  is  a  nonnegative  (not  necessarily  additive)  set  function  on  A. 
The  elements  of  A  will  be  called  flows .  If  V  t  and  V2  are  disjoint  sets  in  V,  then  X(V  txV2) 
represents  the  amount  of  some  commodity  moved  from  V ,  to  V2.  For  finite  networks,  some 
authors  require  flows  to  be  nonnegative  or  antisymmetric  (i.e.,  X(V  iXV2)  =  — X(V2xV  j)).  Still 
others  [Rockafellar.  1984]  do  not  require  either  of  these  conditions.  Our  model  is  general 
enough  to  accommodate  all  three  possibilities. 

Let  Z  be  the  set  of  functions  /  €simp(02,  A)  of  the  form 

/  =  y  (cuj)  -  yCwj) 

for  some  y€simp(fl.V).  Then  Z  is  a  subspace  of  simp(02.  A)  and  will  be  called  the 
differential  space.  If  /  and  y  are  related  as  above,  we  write  /  =  Vy  and  call  /  the 
differential  of  the  potential  y .  The  following  theorem  generalizes  the  feasible  flow  theorem 
[Ford  and  Fulkerson],  [Rockafellar]  for  finite  networks. 


Theorem  6.1 


Let  (ft.  V.  A.  A)  be  a  product-algebra  network,  let  p  be  the  support  functional  of  A,  and 
define  v  on  V  by  v(V )  =  p(V/v  )•  Suppose  the  composition  p  «V  is  vertically  additive  (defined 
in  Section  6.2),  and  let  add(V).  Then  there  exists  a  flow  X€A  such  that 

X(V  xV)-\(VxV)  = /x(V).  V  €V  (6.4) 

if  and  only  if 

fi  ^  v  and  /i(fl)  =  0.  (6.5) 

Proof 

Since  V/v  =  h  xv  ~  Iv  xv  ■  (6.4)  is  equivalent  to  X(V/v  )  =  )-  Therefore.  (6.4) 

implies  n(V  Xp(V/y  )=*'(V  )  and  #t(G)=/>(V/ft)=p( 0)*0.  and  so  (6.5)  is  a  necessary  condi¬ 
tion.  Now  define  the  linear  functional  /  on  Z  by  l  (Vy  )  =  /x(y  )-  Since  Vy  =0  implies  y  is 
equal  to  some  constant  a.  which  implies  fi(y  )  »  atfii O)  =  0,  it  follows  that  Vy  j  —  Vy 2  implies 
fi(y ,)  =  fi(y  2).  and  so  l  is  well  defined.  Now  if  (6.5)  holds,  then  since  />  °V  is  vertically  addi¬ 
tive  and  v((l)  =  0.  Lemma  6.3  implies  that  n  ^  p  »V.  i.e.,  that  l  ^  p  on  Z .  Therefore,  by  the 
Hahn-Banach  Theorem,  there  exists  an  extension  X  of  l  to  simp(ft2.  A)  such  that  X  ^  p .  Since 
X  equals  l  on  Z  .  we  have  X(V/y  )  =  )•  ■ 

The  following  corollary  generalizes  the  max-flow  min-cut  theorem  [Ford  and  Fulkerson] 
for  finite  networks. 

Corollary  6.1 

Let  ( Q.  V.  A,  A),  p .  and  v  be  as  in  Theorem  6.1,  and  suppose  p  »V  is  vertically  additive. 
Let  S  and  D  be  disjoint  sets  in  V.  let  T  denote  the  set  of  flows  X€A  such  that 


A(VxV)-A(VxV)  =  0, 


vcn-isuD}. 


where  V  =  SI— V  .  and  define  v  (A)  =  A(Z)  xD  )  —  A (D  XD  ).  Then 

max  {v  (A):  A€T}  =  inf  [v(V  ):  V  €V.  D  C  V  CS  }.  (6.6) 

Proof 

Let  V  €V  with  D  C V  CS  .  Note  that  for  any  V  .  A(V  xV )  —  A(V  XV  )«  A (V/v  )•  There¬ 
fore.  if  A6T  and  D  CV  CS  ,  then 

v  (A)  =  A(V/d  )  =  A(  V/v  )  ~  A(V/v-o  > 

=  A(VA  )  <  p(V/v)  =  v(V). 

Therefore,  the  left  side  of  (6.6)  is  no  greater  than  the  right  side. 

Let  a  denote  the  right  side  of  (6.6).  To  establish  (6.6)  it  now  suffices  to  show  there  exists 
A€T  such  that  v(A)  —  a.  Let  V'  be  the  subalgebra  of  V  consisting  of  all  sets  V  €V  such  that 
V  C\S  is  either  0  or  S  and  V  C\D  is  either  0  or  D  ,  i.e.,  V’  is  the  algebra  obtained  from  V  by 
making  S  and  D  into  atoms.  Then  (O,  V'.  A.  A)  is  a  product-algebra  network.  Let  fi  be  the 
member  of  add(V)  such  that  /x(D  )  -  a  -  — /*(S  )  and  fi(.V  )  =  0  for  V  C  SI— (S  U£) ).  Then  n^v 
and  so  by  Theorem  6.1  there  exists  A€A  such  that  A(V  XV  )  — A(V  xV  )  -  jjl(V  )  for  all  V  €V.  It 
follows  that  A€T  and  v(A)  =  a.  establishing  (6.6).  ■ 

6.4.  Polymatroid  Networks 

Let  (Q.V.  A.  A)  be  a  product-algebra  network  and  suppose  A  has  the  form 
(A€add(A):  0  <  A  ^  c  }  where  c  is  an  increasing  submodular  set  function  on  A.  Such  a  net¬ 
work  will  be  called  a  polymatroid  network  (since  A  is  in  fact  a  polymatroid  when  SI  is  finite 
[Schrijver.  1984]).  and  c  will  be  called  the  capacity  for  the  network.  We  will  say  that  the 


98 


capacity  c  is  progressively  additive  if  for  any  increasing  finite  collection  V  v  .  .  .  ,V„  of  sets  in 
V. 


n  n  — k  +1  n  _ 

£c(U  ^cCV.xVi). 

k  =1  i  =1  i  =1 


(6.6) 


Define  c  on  simp(  H2.  A)  by 


c  (/  )  =  c  (/  +)  =  J  c  {a>:  /  (<u)  ^  a}  d  a  . 


(6.7) 


We  see  from  Remark  2  in  Section  6.2  that  c  is  the  support  functional  of  A.  It  is  easy  to  verify 
that  (6.6)  is  equivalent  to 


c  ®V 


I'v, 

i=l 


=  £  cv(/Vl). 

<=» 


which  by  Lemma  6.2  is  equivalent  to  c  °V  being  vertically  additive.  Therefore,  c  is  progres¬ 
sively  additive  if  and  only  if  c  «V  is  vertically  additive.  We  can  therefore  apply  Theorem  6.1 
with  p  -~c  and  v(V  )  ="c  ( V^v  )  =  c  (V  xV  )  to  obtain  the  following  feasible  flow  theorem  and 
max-flow  min-cut  corollary. 


Theorem  6.2 

Let  (Q.  V.  A.  A)  be  a  polymatroid  network  whose  capacity  c  is  progressively  additive, 
and  let  /x€add(V).  There  exists  X€A  such  that  X(V  xV )  —  X(V  xV  )  =  p(V )  if  and  only  if 
/t(V )  <  c  (V  xV' )  for  all  V  €V. 


Corollary  6.2 

Let  (0,  V.  A,  A)  be  a  polymatroid  network  whose  capacity  c  is  progressively  additive. 
Let  5,  D,  be  disjoint  sets  in  V,  let  T  denote  the  set  of  flows  X€A  such  that 
X( V  xV )  —  X(V  xV  )  =  0  whenever  V  C II — (5  U D).  and  define  v  (X)  =  \(D  xD  )  —  X(Z>  xD  ). 


'.-■A 


.  i 


••1 


•'*  »*’  ,*•  ,%  .*•  ,  • 
a  f-  «  .*  *•».'  a  *  • 


s.V 


\  % 
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Then 

max  {v  (X):  X€T}  =  inf  {c  (V  XV ):  V  €V.  D  CV  CS  }.  (6.8) 

For  the  remainder  of  the  section,  we  consider  the  special  case  of  finite  networks.  Suppose 
that  (fl.  V.  A.  A)  is  a  polymatroid  network  where  fl  is  a  finite  set,  V  and  A  are  the  power  sets 
of  fl  and  fl2,  respectively,  and  A  =  (X€add(A):  0  ^  X  ^  c  }.  In  this  case,  the  elements  of  fl 
and  fl2  are  called  nodes  and  arcs,  respectively.  Let  F  denote  the  set  of  all  functions 
/  :  fl2-»J?  such  that 

0^  £  /(i./)  <c(A).  A  C fl2. 

U  .)  'fcA 

There  is  clearly  a  one-to-one  correspondence  between  F  and  A  which  is  defined  by 

X(A)  =  £  /  (» ,j ). 

U  .j  %A 

We  next  show  how  to  construct  a  large  class  of  finite  polymatroid  networks  whose  capaci¬ 
ties  are  progressively  additive.  These  networks  are  important  because  they  satisfy  the  max- 
flow  min-cut  property  (6.8)  for  any  choice  of  S  and  D .  If  (i ,/  )  is  an  arc  in  fl2,  then  i  is 

called  the  tail  of  the  arc  and  j  is  called  the  head  of  the  arc.  Let  B  . BM  be  a  collection  of 

sets  in  A  such  that  for  each  m  ,  either  all  arcs  in  Bm  have  a  common  head  or  all  arcs  in  Bm  have 
a  common  tail.  For  each  m  let  rm  be  an  increasing  submodular  function  on  A  such  that 
rm(A)  =  0  if  A  C\Bm  =0.  equivalently,  such  that  rm(A  )  =  rm(A  C\Bm ).  Now  define  the 
capacity  c  on  A  by 

c  (A  )  *  £  rm  (A  ).  (6.9) 

m  =1 


and  let  A  =  |X€add(A):  0  ^  X  ^  c  }.  Note  that  the  sets  B  t, .  .  .  ,  BM  need  not  be  disjoint,  nor 


Cl 


need  their  union  be  equal  to  Cl2.  However,  if  they  have  these  two  properties,  then  it  is  clear 
that  a  set  function  \€add(A)  is  in  A  if  and  only  if  for  each  m  . 

0  4  AC4  )  rm  (A  ).  A  CBm  .  (6.10) 

Theorem  6.3 

The  capacity  c  defined  by  (6.9)  is  progressively  additive,  i.e.,  c  »V  is  vertically  additive. 
Proof 

Since  the  sum  of  progressively  additive  capacities  is  clearly  progressively  additive,  it 
suffices  to  show  that  for  each  m  ,rm  is  progressively  additive,  i.e.,  that  rm  »V  is  vertically  addi¬ 
tive.  where  rm  is  defined  as  in  (6.7)  with  rm  replacing  c  .  So  let  m  G  { 1 . M  }  and  suppose 

that  all  arcs  in  Bm  have  the  common  tail  k  .  (The  case  where  all  arcs  in  Bm  have  a  common 
head  is  similar.)  Since  BmdkxCl  and  rm(A)  =  rm  (A  ClBm)  and  rm  is  increasing,  we  have 
rm{A  )  =  rm(A  D (A  Xfl))  =  rm(k  x[j\(k  ,j)€A  )}.  Therefore, 


rm  “VOy  )  =  f  rm  {(i  ,j  ):  y  O' )— y  (i  )  ^  a}  d  a 

I 

ao 

=  /  rm  (k  xlj :  y  (j )— y  (*  )  ^  a})  d  a. 
o 


.Note  that  k  €{/  :  y  (j  )— y  (k  )  ^  a}  if  and  only  if  o^O.  Therefore,  defining  the  submodular 
function  vm  on  V  by 


^m(V)  = 


rm(*xV). 

0. 


k  ev, 

k  €  V  . 


we  have 


rm  »v(y  )  =  /  Vm{  j:y(j  )-y  (k  )  ^  a}  d  a 

— "OO 

oo 

=  /  >  a)  da. 

— oo 

and  so  rm  »V  =  ^  .  proving  that  rm  »V  is  vertically  additive.  ■ 

Remarks 

(1)  Since  the  capacity  defined  in  (6.9)  is  similar  to  the  polymatroid  model  considered  in 
[Lawler  and  Martel],  we  should  compare  the  two.  In  [Lawler  and  Martel],  a  collection  of  sub- 
modular  functions  and  arc  sets,  essentially  the  same  as  rm  and  Bm ,  m  =1. . .  .  ,  M .  are  given. 
But  instead  of  requiring  0  <  X  <  c  ,  where  c  is  defined  in  (6.9).  they  impose  the  constraint 

given  in  (6.10).  which  is  equivalent  to  0  ^  A  ^  c  only  when  the  sets  B  > . BM  are  disjoint. 

However,  the  model  of  Lawler  and  Martel  can  easily  be  converted  to  our  model  by  modifying 
the  network  in  the  following  way.  For  each  arc  (i  ,j  )  which  belongs  to  two  sets  Bm  and  Bn  (in 
their  model  an  arc  will  never  belong  to  more  than  two  such  sets),  add  a  node  *  to  ft.  replace 

(i  .j  )  with  (t  Jc  )  in  Bm  .  and  replace  (i  .j  )  with  (*  ,j  )  in  B„  .  This  makes  the  sets  B  . . BM 

disjoint.  It  is  not  clear  that  our  model  can  be  easily  be  converted  to  Lawler  and  Martel’s 
model. 

(2)  Once  the  model  in  [Lawler  and  Martel]  is  converted  to  our  model  as  suggested  in  the 
first  remark,  their  version  of  the  max-flow  min-cut  theorem  is  implied  by  our  Corollary  6.2. 
However,  their  method  of  proof,  which  uses  a  flow-augmentation  argument,  is  radically 
different  from  our  method. 

6.5.  General  Vector  Space  Networks 

In  Sections  6.3  and  6.4,  flows  were  taken  to  be  elements  of  add(A).  or  equivalently,  linear 
functionals  on  simp(ft2.  A).  In  this  section,  we  generalize  further,  letting  flows  be  linear  func- 


tionals  on  an  arbitrary  vector  space  Z  .  We  now  present  a  very  general  feasible  flow  theorem. 

Theorem  6.4 

Let  Y  and  Z  be  topological  vector  spaces  with  respective  duals  Y*  and  Z*  [Kelley  and 
Namioka.  1976].  Let  H :  Y  -*Z  be  a  linear  mapping,  and  let  L  =  {l  6 Z* :  l  }.  where  p  is  a 
continuous  sublinear  functional  on  Z.  Let  H*  :Z*~*Y*  be  the  adjoint  of  H .  defined  by 
H*(l)  =  l»H.  Let  r  GY* .  Then  there  exists  l  GL  such  that  H*(l)  =  r  if  and  only  if 
r  ^  p  *H  . 

Proof 

That  r  ^  p»H  is  a  necessary  condition  is  trivial.  Conversely,  suppose  r  ^  p»H ,  and  let 
l '  be  the  linear  functional  on  H  (Y  )  defined  by  l  '{H  (y  ))  =  r  (y  ).  The  inequality  r  <  p»H 
implies  that  r  (y  )=0  whenever  H  (y  )=0.  which  in  turn  implies  that  l '  is  well  defined.  The  ine¬ 
quality  r  ^  p“H  also  implies  that  1‘  <  p  on  HiY).  and  so  by  the  Hahn-Banach  Theorem.  I ' 
has  an  extension  l  €Z*  such  that  l  ^  p .  i.e..  I  €L  .  Since  /  -l '  on  H  (.Y  ).  we  have  H*l  (y  ) 
=  Umy))  =  r(my»  =  r(y).  ■ 

In  some  applications,  we  may  wish  to  generalize  the  condition  H*  {l  )-r  to  Rl  =r  .  where 
R :  Z*  -*Y*  is  not  necessarily  the  adjoint  of  any  linear  mapping.  We  then  have  the  following 
(even  more  general)  feasible  flow  theorem. 

Theorem  6.5 

Let  Y .  Z  ,  L  ,  and  p  be  as  in  Theorem  6.4,  and  let  R :  Z*  be  a  linear  mapping.  Define 
q  on  Y  by  q  (y  )  -  sup  {Rl  (y  ):  l  €L  |.  and  let  r  €T* .  Then  there  exists  l  €L  such  that  Rl  —r  if 
and  only  if  r  ^q  . 

Proof 

Let  U :  Y  -*Z**  be  the  restriction  of  the  adjoint  of  R  to  K.  i.e..  for  y  GY  and  l  €Z* . 
U  (y  )(/ )  -  ( Rl  )(y  ).  Now  define  p '  on  Z**  by  p  '(z  )  -  sup  <1  (z  ):  l  GL  }  where  we  identify  Z* 


with  Z***  in  the  natural  way.  Then  clearly  L  -  | /  €Z* :  Z  and  p  ‘=p  on  Z  .  Now 

p‘*U(y  )  =  sup  {Z(.R(y  )):  Z  €2, }  =  sup  {(.ffZ  )(y  ):  l  £L)  -  q(y). 

Now  apply  Theorem  6.4  with  Z**  in  place  of  Z  and  p '  in  place  of  p .  Since  r  ^ q  —p  '•£/ . 

Theorem  6.4  implies  that  there  exists  Z  €2,  such  that  l«U  -r  .  i.e.,  ^Z  =r  .  ■ 

The  next  feasible  flow  theorem  treats  the  special  case  where  Y  =  simp(ft.  V).  so  that  Y* 
can  be  identified  with  add(V).  In  this  case,  if  q  is  vertically  additive,  then  we  can  express  the 
condition  r  ^q  in  terms  of  set  functions,  as  we  did  in  Theorem  6.1. 

Theorem  6.6 

Let  Z  .  L  .  and  p  be  as  in  Theorems  6.4  and  6.5.  Let  ft  be  any  set.  and  let  V  be  an  algebra 
of  subsets  of  ft.  Let  R :  Z*  -*add(V)  be  a  linear  mapping.  For  Z  €Z* ,  we  denote  Rl  by  pt . 

Define  q  on  f  by  q(y)  -  sup  { Ht(y):l^p)  and  define  v  on  V  by  v{V)=q(Iv)  - 

sup  \p  i  (VU^p).  Suppose  q  is  vertically  additive,  i.e.,  q  -v.  and  let  /x€add(V).  Then  there 
exists  Z  €2,  such  that  p,  =fi  if  and  only  if  p^v  and  /x(ft)=Kft). 

Proof 

Since  q  is  vertically  additive.  Lemma  6.3  implies  that  p^q  if  and  only  if  p^v  and 

p((l)~ v(Q).  The  theorem  therefore  follows  immediately  from  Theorem  6.5.  ■ 


104 


REFERENCES 


Anderson,  E.J.,  P.  Nash  and  A.B.  Philpott,  "A  class  of  continuous  network  flow  problems." 
Math.  ofOper.  Res.,  7  (1982).  pp.  501-514. 

Bertsekas.  D.P..  "Projected  Newton  methods  for  optimization  problems  with  simple  constraints." 
SIAM  J.  Control  and  Optimization,  20  (1982).  pp.  221-246. 

Choquet  G..  "Theory  of  Capacities,"  Arm.  Inst.  Fourier,  vol.  5.  pp.  131-292,  1953. 

Ford.  L.R.  and  D.R.  Fulkerson.  Flows  in  Networks,  Princeton  University  Press.  Princeton.  NJ. 
1962. 

Goldstein.  A.A..  "Convex  programming  in  Hilbert-Space,"  Bulletin  AMS.  70  (1964),  pp.  709- 
710. 

Hajek  B.  and  R.  Ogier.  "Optimal  dynamic  routing  in  communication  networks  with  continuous 
traffic,"  Networks,  14  (1984).pp.  457-487. 

Hassin.  R..  "Minimum-Cost  Flow  with  Set  Constraints."  Networks,  vol.  12,  pp.  1-21.  1982. 

Jodorkovsky.  M.  and  A.  Segall.  "A  maximal  flow  approach  to  dynamic  routing  in  communica¬ 
tion  networks."  EE  Publ.  358.  Technion-Israel  Institute  of  Technology.  August  1979. 

Jodorkovsky.  M.  and  A.  Segall.  "An  optimal  control  approach  to  dynamic  routing  in  communi¬ 
cation  networks.  Part  II:  A  maximal  flow  approach,"  submitted  to  IEEE  Trans.  Automatic  Con¬ 
trol,  July  1981. 

Kelley.  J.L.  and  I.  Namioka,  Linear  Topological  Spaces,  Springer-Verlag.  New  York,  1976. 

Knowles.  G.,  An  Introduction  to  Applied  Optimal  Control,  Academic  Press.  New  York,  1981. 

Lawler.  E.L.  and  C.U.  Martel.  "Computing  Maximal  Polymatroidal  Network  Flows,"  Mathemat¬ 
ics  of  Operation  Research,  vol.7,  no.3.  August  1982. 

Luenberger.  D.G..  Introduction  to  Linear  and  Nonlinear  Programming,  Addison-Wesley.  Read¬ 
ing.  MA.  1973. 

Malhotra.  V.M.,  M.P.  Kumar  and  S.N.  Maheshwari,  "An  0(IVI  )  algorithm  for  finding  maximal 
flows  in  networks."  Inf.  Proc.  Letters,  7  (1978),  pp.  277-278. 

McCormick,  G.P.,  "Anti-zig-zagging  by  bending,"  Management  Science  (Theory),  15  (1969). 
pp.  315-320. 

Moss.  F.H..  "The  application  of  optimal  control  theory  to  dynamic  routing  in  communication 
networks."  Ph.D.  dissertation.  Massachusetts  Institute  of  Technology.  Cambridge: 
Elect.  Syst.  Lab.  Rept.  ESL-R-721.  February  1977. 

Moss.  F.H.  and  A.  Segall.  "An  optimal  control  approach  to  dynamic  routing  in  data  communica¬ 
tion  networks.  Part  II:  Geometrical  interpretation,"  EE  Publ.  319.  Technion-Israel  Institute  of 
Technology,  January  1978. 


Moss,  F.H.  and  A.  Segall,  "An  optimal  control  approach  to  dynamic  routing  in  networks."  IEEE 
Trans.  Automatic  Control,  27  (1982),  pp.  329-339. 

Neumann,  M.M.,  "Flows  in  infinite  networks,"  Proc.  of  the  11th  Winter  School.  Circolo 
Matematico  di  Palermo.  Palermo.  Italy  (1984). 

Neumann.  M.M.,  The  theorem  of  Gale  for  infinite  networks  and  applications,  Springer  Lecture 
Notes  in  Economics  and  Math.  Systems.  Springer-Verlag,  New  York.  (1985). 

Rockafellar.  Network  Flows  and  Monotropic  Optimization,  John  Wiley  and  Sons.  New  York. 
1984. 

Royden.  H.L.,  Real  Analysis,  MacMillan.  New  York.  1968. 

Rudin.  W..  Functional  Analysis,  McGraw-Hill.  New  York.  1973. 

Rudin.  W.,  Real  and  Complex  Analysis,  McGraw-Hill,  New  York.  1974. 

Schrijver.  A..  "Submodular  Functions."  unpublished  notes.  University  of  Amsterdam.  Amster¬ 
dam.  Holland.  1984. 

Segall.  A..  "The  modeling  of  adaptive  routing  in  data  communication  networks."  IEEE 
Trans.  Common.,  25  (1977).  pp.  85-95. 

Segall.  A..  "Optimal  distributed  routing  for  virtual  line-switched  data  networks."  IEEE 
Trans.  Common.,  27  (1979).  pp.  201-209. 

Shats.  S.  and  A.  Segall.  "Open-loop  solutions  for  the  dynamic  routing  problem."  Rept.  LIDS-R- 
922.  Laboratory  for  Information  and  Decision  Systems.  M.I.T..  1980. 


FILMED 

t 

t-86 

K 


DTIC 


